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INTRODUCTION: 

When  we  wrote  our  application  we  proposed  to  identify  the  microRNAs  that  are 
expressed  in  three  different  prostate  cancer  cell  lines,  PC3,  DU145  and  LnCAP. 
Comparison  of  the  microRNA  expression  patterns  was  expected  to  identify  microRNAs 
that  are  specific  to  certain  cancer  cell  lines.  As  the  project  progressed,  we  realized  that 
we  would  learn  more  about  microRNAs  in  prostate  cancer  progression  if  we  compared 
the  microRNA  profiles  of  an  early  stage  and  and  an  advanced  stage  cell  line  derived  from 
the  former.  We  therefore  focused  on  the  androgen-dependent,  non  metastatic  LnCAP 
cells  and  its  derivative  C4-2  that  is  androgen-independent  and  metastatic. 

We  also  wanted  to  characterize  microRNAs  that  are  induced  or  repressed  when  an 
androgen  responsive  prostate  cancer  cell  line,  LnCAP  is  grown  in  the  presence  and 
absence  of  androgen.  Comparison  of  the  microRNA  expression  patterns  was  expected  to 
identify  microRNAs  that  may  have  a  role  in  the  cell’s  response  to  androgens. 


BODY: 

The  Statement  of  Work  and  the  status  of  each  task: 

Task  1.1:  Grow  DU  145,  HeLa  and  HCT1 16  cancer  cell  lines,  isolate  size-selected 
microRNAs  from  each  of  them.  (These  cells  were  substituted  by  C4-2  and  LNCaP 
growing  in  serum.)  Grow  LnCAP  cells  in  the  absence  and  presence  of  androgens  and 
isolate  size-selected  microRNAs.  COMPLETED 

Task  1.2:  Ligate  adapters,  make  cDNAs  of  the  microRNAs,  clone  the  cDNAs  to  make 
libraries  representing  microRNA  complements,  and  sequence  the  clones  in  the  libraries. 
COMPLETED 

Task  1.3:  De-convolute  the  sequences  and  do  bioinformatics  on  public  sequence 
databases  to  identify  predicted  and  novel  microRNAs.  COMPLETED 

NEW  TASK  1.3b:  Assay  the  microRNAs  also  on  Locked  Nucleic  Acid  based  microarrays 
of  microRNAs  to  profile  which  microRNAs  are  altered  between  LNCaP  and  C4-2  and  in 
LNCaP  plus  and  minus  androgen.  COMPLETED 

Task  1 .4:  Determine  which  microRNAs  should  be  pursued  in  task  2.  The  two 
independent  methods  for  evaluating  the  microRNAs  (sequencing  and  microarray)  helped 
us  to  prioritize  which  microRNAs  gave  reproducible  results  and  should  be  pursued. 
COMPLETED 

Task  2.1:  Develop  primer  extension  assays,  Rnase  protection  assays  and  Northern  blot 
assays  for  20-30  microRNAs  prioritized  in  Task  1.4  Quantitative  RT-PCR,  RNase 
Protection  assays  and  Northern  blots  were  used.  COMPLETED 
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Task  2.2:  Assay  levels  of  the  microRNAs  to  confirm  differences  in  their  levels  between 
cell  lines  and  in  response  to  androgen.  COMPLETED 

Task  2.3:  Prioritize  5-8  microRNAs  for  follow  up  in  task  3  COMPLETED 

Task  3.1:  Downregulate  microRNAs  in  PC3  (changed  to  LNCaP  and  C4-2)  cells  and 
assay  effects  on  cell  proliferation.  ONGOING 

Task  3.2:  Upregulate  microRNAs  in  PC3  (changed  to  LNCaP  and  C4-2)  cells  and  assay 
effects  on  cell  proliferation.  COMPLETED 

Task  3.3:  Alter  microRNA  levels  in  LnCAP  cells,  treat  cells  with  androgens  and  harvest 
total  RNA  at  specified  time-points.  ELIMINATED.  The  microRNA  alteration  was  not 
done.  Instead  we  went  straight  to  Task  3.4  because  a  short-cut  became  apparent. 

Task  3.4:  Microarray  hybridization  of  the  RNA  from  task  3.3.  Instead,  the  genes 
changed  by  androgens  were  obtained  from  microarray  experiments  in  public  databases 
where  LNCaP  cells  were  grown  in  the  presence  and  absence  of  androgen.  COMPLETED 

Task  3.5:  Analysis  of  results  from  Task  3.4  to  determine  whether  the  tested  microRNAs 
modulate  the  expression  of  genes  affected  by  androgens.  COMPLETED 

Deliverables 

At  12  months:  Annual  report  for  year  1  DELIVERED 
At  24  months:  Annual  report  for  year  2  DELIVERED 
At  36  months:  Final  report  of  project.  CURRENT 

2-3  manuscripts  for  peer  review  reporting  on  the  microRNAome  of  prostate  cancer  cells 
with  an  emphasis  on  microRNAs  that  are  critical  for  cell  proliferation  and/or  androgen- 
driven  gene  expression  program.  COMPLETED 

Summary  of  our  results. 

Papers:  Two  manuscripts  have  been  prepared  (attached  as  appendix): 

1)  “ MicroRNAs  in  the  androgen-independence  and  androgen  regulation  of  prostate 
epithelial  cell  lines”.  Under  revision  after  review.  Copy  included  in  the  Appendix  as 
LEE1.  Figures  and  Tables  in  the  following  section  refer  to  Figures  and_  Tables  in  the 
text  o£ LEE1_  and_  in  the.  Supplemental  material  (or  LEEL 

A  systematic  comparison  of  the  short  RNA  universe  of  two  prostate  cancer  cell 
lines  was  performed  by  cloning  and  deep  sequencing  to  identify  microRNAs  whose 
expression  may  correlate  with  androgen  dependence/independence  or  microRNAs  that 
are  regulated  by  androgens  (Tables  1  and  2  in  attached  paper,  Leel).  The  results  were 
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confirmed  by  two  independent  methods  of  quantitating  microRNAs:  locked  nucleic  acid 
microarrays,  (Supplemental  Tables  SI,  S2  and  S3)  and  microRNA-specific 
qPCR/RNase  protection/Northem  assays  (Figure  1  and  Supplemental  Figure  S5). 
Changes  that  were  consistently  detected  by  all  three  methods  are  the  decrease  of  five 
microRNAs  and  increase  of  three  microRNAs  in  the  androgen- independent  C4-2  cell  line 
relative  to  its  androgen-dependent  predecessor,  LNCaP  (Fig.  1A).  miR-125b,  100,  -99a, 
-99b  and  19b  were  downregulated  in  C4-2  and  miR-196b,  -9  and  -557  were  upregulated 
in  C4-2.  Student’s  t-test  showed  that  all  the  changes  measured  by  Q-RT-PCR  except  that 
of  miR-19b  had  a  p<0.05.  For  the  Q-RT-PCR  measurement  of  the  decrease  of  mir-19b  in 
C4-2,  the  p=0.18.  Nevertheless,  we  have  decided  to  keep  this  microRNA  in  the  list  of 
microRNAs  decreased  in  C4-2  relative  to  LNCaP  because  the  difference  was  noted 
independently  in  the  microarrays  and  in  the  cloning/sequencing  approach. 

Ectopic  expression  of  the  microRNAs  that  were  downregulated  in  C4-2  decreased 
the  growth  of  C4-2  in  the  absence  of  androgen  but  not  in  its  presence  (Fig.  2A  and  B), 
suggesting  that  microRNA  changes  contribute  to  the  androgen-independence.  Except  for 
the  transfection  of  single  microRNAs  miR-99b  and  -125b  in  Fig.  2B,  transfection  of  all 
five  microRNAs  (Fig.  2A)  or  of  single  microRNAs  miR-99a,  -100  and  -19b  (Fig.  2B) 
repressed  C4-2  proliferation  in  androgen-depleted  medium  (p<0.05,  Student’s  t  test).  All 
five  downregulated  microRNAs  have  been  reported  to  be  downregulated  in  clinically 
advanced  prostate  cancer  [1],  suggesting  that  they  constitute  a  microRNA  signature  of 
prostate  cancer  progression. 

Four  microRNAs  were  consistently  shown  to  be  downregulated  in  multiple  assays 
by  androgens:  miR-221,  -196b,  -125b  and  -99a  (Fig.  ID,  Supplemental  Table  S3,  Fig. 
SI  and  S4).  The  difference  in  the  level  of  microRNA  in  the  highest  concentration  of 
R1881  relative  to  no  androgen  was  found  to  be  significant  in  all  four  cases  (p<0.05, 
Student’s  t  test).  We  were  disappointed  by  the  failure  to  detect  any  microRNAs  induced 
by  androgens  and  plan  to  re-visit  the  question  by  additional  microRNA  microarrays  on 
other  androgen-dependent  cell  lines.  This  is  why  we  have  been  unable  to  ascertain 
whether  transfection  of  androgen-induceed  microRNAs  change  the  behavior  of  the  cells 
(the  former  Task  3.3).  We  have,  however,  examined  gene  expression  changes  induced  by 
androgens  (new  Tasks  3.4  and  3.5)  to  ascertain  whether  induction  of  a  gene  by  androgen 
could  be  explained  by  the  repression  of  the  four  microRNAs.  Intersection  of  the  mRNA 
changes  induced  by  androgens  (available  in  public  microarray  datasets)  with  the 
computationally  predicted  targets  of  the  four  microRNAs  (the  TargetScan  and  MiRanda 
databases)  yielded  four  mRNAs  that  could  be  regulated  by  the  microRNAs.  The  genes 
that  will  be  studied  are  a  multi-drug  transporter  involved  in  drug  resistance  (ABCC4),  an 
anti-apoptotic  transcriptional  factor  (Stat3  ),  PTPN 1 8  (a  protein  tyrosine  phosphatase 
involved  in  signal  transduction)  and  CPD  (carboxypeptidase  D,  a  gene  that  regulates 
intracellular  NO  levels).  The  first  three  are  predicted  to  be  repressed  by  miR-125b  and 
the  last  one  by  miR-196b.  All  the  genes  are  induced  by  androgens  and  both  the 
microRNAs  are  repressed  by  androgens. 

Despite  the  broad  concordance  in  changes  in  microRNAs  detected  by 
cloning/sequencing  versus  microarrays  in  nearly  four  hundred  microRNAs  studied  , 
there  was  a  significant  discrepancy  in  the  changes  reported  by  cloning/sequencing  versus 
microarray  hybridization  for  six  microRNAs  (Supplemental  Fig.  S6).  miR-16,  -21,  200c, 
106a,  20a  and  -30b  appeared  to  be  5-50  fold  higher  in  androgen  treated  LNCaP  than  in 
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androgen-depleted  cells  in  the  cloning/sequencing  experiments.  Yet,  on  the  microarrays 
these  microRNAs  were  equally  abundant  in  the  two  types  of  cells.  Q-RT-PCR  confirmed 
the  equal  abundance  (Fig.  S6).  This  suggests  that  there  is  a  systematic  loss  of  these 
microRNAs  when  they  were  cloned/sequenced  from  androgen-depleted  cells.  Using 
miR-2 1  as  an  example,  we  followed  the  yield  of  this  microRNA  at  various  steps  in  the 
cloning  procedure  (Fig.  3),  and  discovered  that  an  as  yet  unknown  3’  end  modification  of 
miR-2 1  in  androgen-depleted  cells  decreases  3’  adaptor  ligation  and  thus  decreases  the 
cloning  frequency. 

Finally,  cloning  and  sequencing  reveals  a  large  number  of  non-micro-small 
RNAs  (nmsRNAs)  of  unknown  function  that  are  expressed  abundantly  in  the  prostate 
epithelial  cells.  Examples  of  these  nmsRNAs  are  given  in  Table  3  of  Leel  and  in 
Supplemental  Fig.  S8.  As  more  of  these  nmsRNAs  are  characterized,  we  expect  that 
novel  functions  will  be  attributed  to  them  and  that  some  of  them  may  be  as  useful  for 
profiling  cancers  as  microRNAs  or  mRNAs 

2)  “A  novel  class  of  small  RNAs,  tRNA-derived  RNA  fragments  (tRFs),  is  highly 
expressed  in  proliferating  cells  and  required  for  cell  proliferation.  ”  Reviewed  at  Genes 
&  Development,  currently  being  revised  for  re-submission.  Copy  included  in  Appendix 
as  LEE2.  Figures  and  Tables  in  the  following  section  refer  to  LEE2  and_  Supplemental 
Material  of  LEE2. 

New  types  of  small  RNAs  distinct  from  microRNAs  (miRNAs)  are  progressively 
being  discovered  in  various  organisms.  In  the  course  of  cloning  and  sequencing  17-26 
base  long  RNAs  from  prostate  cancer  cell  lines  we  discovered  a  novel  class  of  short 
RNAs,  the  tRFs.  A  significant  number  of  the  sequences  are  derived  from  precise 
processing  at  the  5'-  or  3 '-end  of  mature  or  precursor  tRNAs,  to  form  three  series  of  tRFs, 
the  -5,  -3  and  -1  series  (Table  1  and  Fig.  1  of  LEE2;  also  Supplemental  Material, 
Table  SI,  S2  and  Fig.  S1-S3).  tRF-5  molecules  are  generated  by  the  5’-most  21-22 
nucleotides  being  cleaved  from  mature  tRNAs.  tRF-3  molecules  are  generated  by  the 
cleavage  of  the  3’-most  21-22  nucleotides  from  mature  tRNAs.  tRF-1  molecules  are 
generated  from  precursors  of  tRNAs  by  the  cleavage  that  releases  the  3’  end  of  the 
mature  tRNA.  These  sequences,  collectively  named  as  tRFs  (tRNA-derived  RNA 
Eragments),  constitute  a  class  of  short  RNAs  that  are  second-most  abundant  to  miRNAs. 

Northern  hybridization  and  qRT-PCR  assays  independently  measured  the  levels 
of  at  least  16  tRFs  (Fig.  2  and  S4).  Their  ease  of  detection  is  comparable  to  that  of  some 
of  the  best-studied  microRNAs. 

The  best  evidence  that  the  tRFs  are  not  random  byproducts  of  the  turnover  of 
tRNAs  or  their  precursors  would  be  if  knocking  down  a  specific  tRF  produces  a 
phenotype  suggestive  of  a  specific  biological  function.  Towards  that  end,  we  focused  on 
tRF-1001,  derived  from  the  3'-end  of  a  Ser-TGA  tRNA  precursor  transcript  that  is  not 
retained  in  the  mature  tRNA.  tRF-1001  is  expressed  highly  in  a  wide  range  of  cancer  cell 
lines  but  much  less  in  tissues  and  its  expression  in  cell  lines  was  positively  correlated 
with  cell  proliferation  (Fig.  2).  siRNA-mediated  knockdown  of  tRF-1001  impaired  cell 
proliferation  with  the  specific  accumulation  of  cells  in  G2  (Fig.  3C),  phenotypes  that 
were  reversed  by  co-introducing  a  synthetic  2'-0-methyl  tRF-1001  oligoribonucleotide 
resistant  to  the  siRNA  (Fig.  3D).  Although  tRF-1001  is  localized  in  the  cytoplasm  and 
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similar  in  size  to  miRNAs  and  siRNAs,  it  did  not  repress  a  target  mRNA  with  sites  that 
are  complementary  to  the  tRF  molecule  (Fig.  4),  suggesting  that  tRF-1001  does  not  enter 
into  RISC  or  miRNP  complexes.  Our  data  suggest  that  tRFs  are  not  random  by-products 
of  tRNA  degradation  or  biogenesis  but  an  abundant  and  novel  class  of  short  RNAs  with 
precise  sequence  structure  that  have  specific  expression  patterns  and  specific  biological 
roles. 

3)  MicroRNAs  in  Cancer,  Annual  Review  of  Pathology,  2009  4: 199-227  (A  review) 

Within  the  past  few  years,  studies  on  microRNA  (miRNA)  and  cancer  have  burst 
onto  the  scene.  Profiling  of  the  miRNome  (global  miRNA  expression  levels)  has  become 
prevalent,  and  abundant  miRNome  data  are  currently  available  for  various  cancers.  The 
pattern  of  miRNA  expression  can  be  correlated  with  cancer  type,  stage,  and  other  clinical 
variables,  so  miRNA  profiling  can  be  used  as  a  tool  for  cancer  diagnosis  and  prognosis. 
miRNA  expression  analyses  also  suggest  oncogenic  (or  tumor-suppressive)  roles  of 
miRNAs.  miRNAs  play  roles  in  almost  all  aspects  of  cancer  biology,  such  as 
proliferation,  apoptosis,  invasion/metastasis,  and  angiogenesis.  Given  that  many  miRNAs 
are  deregulated  in  cancers  but  have  not  yet  been  further  studied,  it  is  expected  that  more 
miRNAs  will  emerge  as  players  in  the  etiology  and  progression  of  cancer.  Here  we  also 
discuss  miRNAs  as  a  tool  for  cancer  therapy. 

Meeting  Abstracts: 

1)  ImPact,  2007,  Atlanta 

MicroRNA  profile  in  prostate  cancer  cells  and  response  to  androgen  depletion 
Anindya  Dutta;  Yong  Sun  Lee;  Hak  Kyun  Kim 

MicroRNAs  are  short  single-stranded  RNAs  of  18-22  bases  length  that  are  produced  by 
the  processing  of  specific  transcripts  by  cellular  RNAses  Drosha  and  Dicer.  They 
regulate  the  expression  of  a  wide  variety  of  genes  by  annealing  to  the  3  prime 
untranslated  regions  (UTRs)  of  the  mRNAs  and  preventing  the  expression  of  the  genes. 

In  order  to  understand  whether  microRNAs  are  important  for  prostate  cancer  biology  we 
have  begun  examining  the  profile  of  microRNAs  in  prostate  cancer.  We  have  cloned 
short  RNAs  of  20-28  bases  from  three  different  prostate  cancer  cell  lines,  PC3,  LnCAP 
and  DU  145.  The  clones  have  been  sequenced  to  identify  the  microRNAs  that  are 
expressed  in  these  cells.  Over  a  hundred  microRNAs  have  been  identified  to  be 
expressed  in  prostate  cancer  cells,  and  tens  of  new  short  RNAs  discovered  that  could  be 
novel  microRNAs  or  some  other  type  of  cellular  short  RNAs.  A  distinct  microRNA 
expression  profile  separates  the  PC3  and  DU145  cells  from  LnCAP  cells,  consistent  with 
the  different  lineage  of  the  latter  type  of  cells.  A  subset  of  the  microRNAs  are  expressed 
specifically  in  an  androgen-dependent  manner  in  LnCAP  cells,  raising  the  possibility  that 
changes  in  expression  of  microRNAs  may  be  induced  by  androgens  and  that  microRNAs 
may  contribute  to  the  androgen  dependent  cell  proliferation  program  in  prostate  epithelial 
cells.  Finally,  overexpression  of  specific  microRNAs  modulate  the  proliferation  of 
cancer  cells,  so  that  a  systematic  examination  of  the  substrates  of  these  microRNAs  will 
reveal  genes  important  for  cancer  cell  proliferation. 
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2)  Forbeck  Symposium  on  microRNAs  and  cancer,  Hilton  Head,  SC,  2007 
MicroRNAs  in  cell  quiescence  and  cancer 
Anindva  Dutta 


We  are  examining  two  systems  to  discover  microRNAs  that  regulate  cell 
proliferation.  In  the  first  approach,  we  focused  on  miRNAs  induced  during 
differentiation  of  C2C12  myoblasts  into  myotubes  in  vitro.  One  of  the  microRNAs 
induced  during  differentiation,  miR-206,  suppresses  cell  proliferation.  By  surveying 
mRNAs  that  are  downregulated  by  miR-206  for  genes  predicted  to  be  targets  of  the 
microRNA,  we  discovered  the  mRNA  of  DNA  polymerase  alpha  pi 80  is  a  direct  target 
that  is  destabilized  upon  miR-206  transfection.  We  now  report  that  additional 
microRNAs  induced  during  muscle  differentiation  target  other  cell-cycle  regulators. 

After  noting  that  targeted  mRNAs  are  often  destabilized  by  microRNAs,  we  tested 
whether  global  downregulation  of  microRNAs  by  knockdown  of  Dicer  or  Drosha  can 
help  identify  microRNA-repressed  mRNAs.  The  principle  was  proved  by  the  discovery 
that  HMGA2  oncogene  is  repressed  by  the  growth-suppressive  let-7  microRNA. 

We  have  sought  other  examples  of  microRNAs  and  short  RNAs  that  regulate  cell 
proliferation  by  cloning  short  RNAs  from  androgen-dependent  prostate  cancer  cells 
grown  in  the  presence  or  absence  of  androgens.  The  frequency  of  a  number  of  clones 
present  in  the  libraries  is  changed  by  androgens  and  we  are  currently  testing  whether  the 
changes  reflect  changes  in  the  abundance  of  the  microRNAs  upon  androgen  depletion. 
Over  30-40%  of  the  short  RNAs  cloned,  however,  do  not  correspond  to  known 
microRNAs  and  at  least  some  are  produced  by  cleavage  of  known  mRNAs  and 
noncoding  RNAs.  We  are  curious  whether  other  groups  cloning  short  RNAs  from  other 
sources  are  also  observing  this  phenomenon. 

3)  AACR  meeting  on  Molecular  Diagnostics,  Philadelphia,  2008 

MicroRNAs  and  non-micro-short  RNAs  (nmsRNAs)  in  cell  quiescence  and  cancer 

Anindya  Dutta 

We  are  using  three  approaches  to  discover  microRNAs  that  regulate  cell  proliferation.  In 
the  first  approach,  we  focused  on  miRNAs  induced  during  differentiation  of  C2C12 
myoblasts  into  myotubes  in  vitro.  One  of  the  microRNAs  induced  during  differentiation, 
miR-206,  suppresses  cell  proliferation.  By  surveying  mRNAs  that  are  downregulated  by 
miR-206  for  genes  predicted  to  be  targets  of  the  microRNA,  we  discovered  the  mRNA  of 
DNA  polymerase  alpha  pi 80  is  a  direct  target  that  is  destabilized  upon  miR-206 
transfection.  We  now  report  that  additional  microRNAs  induced  during  muscle 
differentiation  target  other  cell-cycle  regulators.  After  noting  that  targeted  mRNAs  are 
often  destabilized  by  microRNAs,  in  the  second  approach,  we  tested  whether  global 
downregulation  of  microRNAs  by  knockdown  of  Dicer  or  Drosha  can  help  identify 
microRNA-repressed  mRNAs.  The  principle  was  proved  by  the  discovery  that  HMGA2 
oncogene  is  repressed  by  the  growth-suppressive  let-7  microRNA.  Indeed,  chromosomal 
translocations  in  lipomas  and  leiomyomas  de-repress  HMGA2  by  deleting  the  3’UTR 
that  is  normally  repressed  by  let-7.  Conversely,  downregulation  of  let-7  has  been  noted  in 
lung  cancers  and  large  leiomyomas  and  is  correlated  with  over-expression  of  the  HMGA2 
oncogene. 
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In  the  third  approach  we  have  cloned  short  RNAs  from  androgen-dependent  prostate 
cancer  cells  grown  in  the  presence  or  absence  of  androgens  and  subjected  them  to  ultra- 
high-throughput  sequencing.  The  frequency  of  a  number  of  clones  present  in  the  libraries 
is  changed  by  androgens  and  we  have  identified  several  microRNAs  that  change  upon 
androgen-depletion.  Over  30-40%  of  the  short  RNAs  cloned,  however,  do  not  correspond 
to  known  microRNAs  and  are  produced  by  cleavage  of  known  mRNAs  and  noncoding 
RNAs.  The  diversity  and  abundance  of  the  non-micro-short  RNAs  (nmsRNAs)  contrast 
with  how  little  is  known  of  their  function  and  suggest  that  much  remains  to  be  done 
before  we  understand  the  biological  functions  of  many  short  RNAs  present  in  the  cell.  In 
an  extension  of  the  project  we  have  identified  several  microRNAs  that  are  reproducibly 
altered  in  expression  level  between  the  androgen-dependent  LnCAP  cells  and  androgen- 
independent  derivative  C4-2  cells.  Interestingly,  the  microRNAs  decreased  in  C4-2  cells 
were  also  reported  by  others  to  be  decreased  as  prostate  cancers  advance  in  Gleason  score 
or  metastasize.  This  raises  the  possibility  that  a  microRNA  signature  may  be  established 
for  following  the  molecular  stage  of  prostate  cancer  and  for  predicting  its  predilection  to 
recur  in  an  androgen-independent  form. 

4)  Human  Genome  Organization.  Annual  meeting,  Hyderabad,  2008 
MicroRNAs  and  non-micro-short  RNAs  (nmsRNAs)  in  cell  quiescence  and  cancer 
Anindya  Dutta 

We  are  using  three  approaches  to  discover  microRNAs  that  regulate  cell  proliferation.  In 
the  first  approach,  we  focused  on  miRNAs  induced  during  differentiation  of  C2C12 
myoblasts  into  myotubes  in  vitro.  One  of  the  microRNAs  induced  during  differentiation, 
miR-206,  suppresses  cell  proliferation.  By  surveying  mRNAs  that  are  downregulated  by 
miR-206  for  genes  predicted  to  be  targets  of  the  microRNA,  we  discovered  the  mRNA  of 
DNA  polymerase  alpha  pi 80  is  a  direct  target  that  is  destabilized  upon  miR-206 
transfection.  We  now  report  that  additional  microRNAs  induced  during  muscle 
differentiation  target  other  cell-cycle  regulators.  After  noting  that  targeted  mRNAs  are 
often  destabilized  by  microRNAs,  in  the  second  approach,  we  tested  whether  global 
downregulation  of  microRNAs  by  knockdown  of  Dicer  or  Drosha  can  help  identify 
microRNA-repressed  mRNAs.  The  principle  was  proved  by  the  discovery  that  HMGA2 
oncogene  is  repressed  by  the  growth-suppressive  let-7  microRNA.  Indeed,  chromosomal 
translocations  in  lipomas  and  leiomyomas  de-repress  HMGA2  by  deleting  the  3’UTR 
that  is  normally  repressed  by  let-7.  Conversely,  downregulation  of  let-7  has  been  noted  in 
lung  cancers  and  large  leiomyomas  and  is  correlated  with  over-expression  of  the  HMGA2 
oncogene. 

In  the  third  approach  we  have  cloned  short  RNAs  from  androgen-dependent  prostate 
cancer  cells  grown  in  the  presence  or  absence  of  androgens  and  subjected  them  to  ultra- 
high-throughput  sequencing.  The  frequency  of  a  number  of  clones  present  in  the  libraries 
is  changed  by  androgens  and  we  have  identified  several  microRNAs  that  change  upon 
androgen  depletion.  Over  30-40%  of  the  short  RNAs  cloned,  however,  do  not  correspond 
to  known  microRNAs  and  are  produced  by  cleavage  of  known  mRNAs  and  noncoding 
RNAs.  The  diversity  and  abundance  of  the  non-micro-short  RNAs  (nmsRNAs)  contrast 
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with  how  little  is  known  of  their  function  and  suggest  that  much  remains  to  be  done 
before  we  understand  the  biological  functions  of  many  short  RNAs  present  in  the  cell. 

Personnel  receiving  pay  from  this  project. 

Anindya  Dutta 
Hak  Kyun  Kim 
Yong  Sun  Lee 
Ankit  Malhotra 
Mirela  Matecic 
Dandan  Sun 

KEY  RESEARCH  ACCOMPLISHMENTS: 

-  Discovering  microRNAs  that  are  altered  when  prostate  cancer  cells  become  androgen 
independent.  Increased  in  C4-2  (more  advanced):  miR-196b,  -557  and  -9.  Decreased  in 
C4-2:  miR  -100.  -125b.  -19b,  -99a  and  -99b 

-  Discovering  three  microRNAs  that  are  repressed  when  androgen  is  added  to  an 
androgen-dependent  prostate  cancer  cells:  miR-221,  -196b,  -125b. 

-  Discovery  that  certain  microRNAs  (miR- 16,  -21,  -200c,  -106a,  -20a  and  miR-30b)  have 
a  3’  end  modification  in  androgen-depleted  cells  that  prevent  their  cloning. 

-  Discovery  of  a  novel  class  of  short  RNAs,  tRNA-derived  fragments  (tRFs)  with 
biological  function. 

-Comparison  of  miRNome  by  two  different  profiling  approaches:  cloning  and  deep 
sequencing  and  microarrays 

REPORTABLE  OUTCOMES: 

-  One  published  paper: 

1)  Lee  YS  and  Dutta  A.  MicroRNAs  in  Cancer,  Annual  Review  of  Pathology,  2009  4: 
199-227 

-Four  Meeting  Abstracts: 

1)  ImPact,  2007,  Atlanta 

MicroRNA  profile  in  prostate  cancer  cells  and  response  to  androgen  depletion.  Anindya 
Dutta;  Yong  Sun  Lee;  Hak  Kyun  Kim 

2)  Forbeck  Symposium  on  microRNAs  and  cancer,  Hilton  Head,  SC,  2007 
MicroRNAs  in  cell  quiescence  and  cancer  .  Anindya  Dutta 

3)  AACR  meeting  on  Molecular  Diagnostics.  Philadelphia.  2008 
MicroRNAs  and  non-micro-short  RNAs  (nmsRNAs)  in  cell  quiescence  and  cancer. 
Anindya  Dutta 

4)  Human  Genome  Organization.  Annual  meeting,  Hyderabad,  2008 
MicroRNAs  and  non-micro-short  RNAs  (nmsRNAs)  in  cell  quiescence  and  cancer. 
Anindya  Dutta 

-  Data  generated  used  as  preliminary  data  for  an  R01  application  to  the  National  Cancer 
Institute. 
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-  Dr.  Yong  Sun  Lee  was  selected  after  a  National  Search  to  be  an  Assistant  Professor  of 
Biochemistry  at  the  University  of  Texas  Medical  Branch  in  Galveston. 

-Dr.  Hak  Kyun  Kim  obtained  his  Ph.D.  degree  from  Seoul  National  University  based  on 
his  work  on  this  project.  He  is  currently  a  postdoctoral  fellow  at  Stanford  University 
School  of  Medicine. 

-Mr.  Ankit  Malhotra  will  obtain  a  Ph.D.  in  Biochemistry  from  the  University  of  Virginia 
at  least  partly  on  the  basis  of  his  work  in  this  project 

-Dr.  Anindya  Dutta  gave  talks  based  on  this  work: 

2007:  The  Forbeck  Symposium  on  microRNAs  in  Hilton  Head,  SC, 

2008:  The  Forbeck  Scholars  Retreat  in  Lake  Geneva,  WI 

2008:  AACR  special  meeting  on  microRNAs  in  cancer  in  Philadelphia,  PA. 

CONCLUSION: 

The  papers  that  are  now  going  through  various  phases  of  peer-review  are  expected 
to  be  published  in  high  impact  journals.  Discovering  that  a  handful  of  microRNAs  are 
altered  as  prostate  cancer  cells  progress  has  allowed  us  to  now  focus  in  on  these  to 
ascertain  whether  they  contribute  to  the  cancer  progression.  We  were  very  pleased  that 
the  microRNAs  we  identified  were  also  identified  as  being  altered  in  a  similar  manner 
when  clinical  prostate  cancer  tissues  are  examined  [1]. 

The  microRNAs  repressed  by  androgens  will  allow  us  to  ascertain  whether  this  is 
a  new  way  by  which  androgens  alter  the  gene  expression  program  in  prostate  cancer  cells. 

The  discovery  of  the  tRFs  was  unexpected.  These  short  RNAs  are  sometimes 
more  abundant  than  bonafide  microRNAs  and  we  have  already  demonstrated  that  some 
of  them  have  biological  functions.  The  biogenesis  and  mechanism  of  action  of  these 
tRFs  are  expected  to  take  up  a  lot  of  future  time  in  my  laboratory  and  in  Dr.  Yong  Sun 
Lee’s  laboratory.  Whether  they  contribute  to  prostate  cancer  progression  will  be  an 
important  topic  of  inquiry. 

We  are  currently  embarking  on  characterizing  which  genes  are  targeted  by  the 
microRNAs  we  have  identified.  We  are  also  beginning  phenotypic  assays  for  cell 
migration  and  invasion  in  C4-2  to  ascertain  whether  overexpression  or  depletion  of  the 
altered  microRNAs  will  alter  their  migration  or  invasion  phenotypes.  If  changing  the 
microRNAs  alter  prostate  cancer  cell  behavior,  we  hope  to  go  to  xenograft  assays  to 
demonstrate  that  microRNA  alterations  can  suppress  metastasis. 

REFERENCES: 

1.  Mattie,  M.D.,  et  al.,  Optimized  high-throughput  microRNA  expression  profiling 
provides  novel  biomarker  assessment  of  clinical  prostate  and  breast  cancer 
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Abstract 

Within  the  past  few  years,  studies  on  microRNA  (miRNA)  and  cancer 
have  burst  onto  the  scene.  Profiling  of  the  miRNome  (global  miRNA 
expression  levels)  has  become  prevalent,  and  abundant  miRNome  data 
are  currently  available  for  various  cancers.  The  pattern  of  miRNA  ex¬ 
pression  can  be  correlated  with  cancer  type,  stage,  and  other  clinical 
variables,  so  miRNA  profiling  can  be  used  as  a  tool  for  cancer  diag¬ 
nosis  and  prognosis.  miRNA  expression  analyses  also  suggest  onco¬ 
genic  (or  tumor-suppressive)  roles  of  miRNAs.  miRNAs  play  roles  in 
almost  all  aspects  of  cancer  biology,  such  as  proliferation,  apoptosis, 
invasion/metastasis,  and  angiogenesis.  Given  that  many  miRNAs  are 
deregulated  in  cancers  but  have  not  yet  been  further  studied,  it  is  ex¬ 
pected  that  more  miRNAs  will  emerge  as  players  in  the  etiology  and 
progression  of  cancer.  Here  we  also  discuss  miRNAs  as  a  tool  for  cancer 
therapy. 
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siRNA:  small 
interfering  RNA 

pre-miRNA: 

precursor  microRNA 

UTR:  untranslated 
region  of  messenger 
RNA 


SYNOPSIS 

Within  the  past  decade,  there  was  a  major  ad¬ 
vance  in  biology:  the  discovery  of  small  RNAs, 
including  microRNA  (miRNA)  and  small  in¬ 
terfering  RNA  (siRNA),  as  highlighted  in  the 
December  2002  issue  of  Science  (1).  Since  the 
discovery  of  RNA  interference  (RNAi)  in  ne¬ 
matodes  (2),  siRNA  has  provided  a  technical 
breakthrough  for  short-term  genetics  in  mam¬ 
malian  systems.  The  big  impact  of  small  RNAs 
was  well  celebrated  by  the  2006  Nobel  Prize 
awarded  to  the  two  scientists  who  discovered 
RNAi. 

Moreover,  miRNAs  shed  new  light  on  the 
posttranscriptional  regulation  of  gene  expres¬ 
sion.  miRNAs  were  first  discovered  in  worms 
(3,  4),  and  they  were  later  observed  in  a  num¬ 
ber  of  animals,  plants,  and  viruses.  During  the 
past  several  years  the  miRNA  field  has  ex¬ 
panded,  with  many  recent  publications  impli¬ 
cating  miRNAs  in  diverse  cellular  processes. 

Cancer  is  a  complex  genetic  disease  caused 
by  the  accumulation  of  mutations,  which  lead 
to  deregulation  of  gene  expression  and  uncon¬ 
trolled  cell  proliferation.  Given  the  wide  im¬ 
pact  of  miRNAs  on  gene  expression,  it  is  not 
surprising  that  a  number  of  miRNAs  have  been 
implicated  in  cancer.  In  this  review,  we  present  a 
comprehensive  discussion  of  the  links  between 
miRNA  and  cancer. 

INTRODUCTION 

MicroRNAs  (miRNAs)  are  small  noncoding 
regulatory  RNAs  ranging  in  size  from  17  to 
25  nucleotides  (see  miRBase,  http://microma. 
sanger.ac.uk/).  The  definition  of  miRNAs  is 
based  on  their  generation  by  the  action  of  Dicer, 
an  RNase  that  processes  hairpin-structured 
precursors  (known  as  pre-miRNA)  into  ma¬ 
ture  miRNAs  (5).  miRNAs  posttranscription- 
ally  repress  gene  expression  by  recognizing 
complementary  target  sites  in  the  V  untrans¬ 
lated  region  (UTR)  of  target  messenger  RNAs 
(mRNAs). 

Eight  years  after  the  first  miRNA  was  re¬ 
ported  by  Ambros’s  and  Ruvkun’s  groups  (3,  4) 
in  1993,  the  miRNA  era  began  in  earnest  when 


three  groups  identified  tens  of  small  RNAs  from 
Caenorhabditis  elegans ,  Drosophila ,  and  human 
(6-8).  During  the  past  seven  years,  the  number 
of  miRNAs  identified,  as  well  as  the  number 
of  publications  describing  them,  has  expanded 
enormously. 

Genes,  Sequences, 

and  Genomic  Organization 

Currently,  more  than  5000  miRNAs  from  over 
50  organisms  are  registered  in  the  miRBase 
database  (8a).  Five  hundred  thirty- three  human 
miRNAs  are  now  known,  but  this  number  may 
increase,  as  bioinformatics  studies  predict  that 
up  to  1000  miRNAs  exist  (9,  10).  miRNAs  are 
named  as  miR-  plus  numbers  (e.g.,  miR-125), 
although  there  are  a  few  exceptions.  miRNAs  of 
similar  sequence  are  usually  distinguished  by  an 
additional  letter  following  the  miRNA  number 
(e.g.,  miR- 12  5b).  A  miRNA  of  identical  mature 
sequence  may  appear  at  several  genomic  loci 
with  different  precursor  sequences.  In  such  a 
case,  the  different  miRNA  genes  would  be  dis¬ 
tinguished  by  the  addition  of  another  number 
at  the  end  of  the  sequence  (e.g.,  miR-125b-l). 

In  humans,  approximately  one-third  of 
miRNAs  are  organized  into  clusters.  A  given 
cluster  is  likely  to  be  a  single  transcriptional 
unit,  suggesting  a  coordinated  regulation  of 
miRNAs  in  the  cluster.  In  silico  analysis  has 
revealed  that  more  than  half  of  the  clus¬ 
ters  contain  two  or  more  miRNAs  of  simi¬ 
lar  sequence  (11).  However,  it  is  very  rare  for 
miRNAs  of  an  identical  mature  sequence  to  be 
duplicated  in  a  cluster.  This  genomic  organi¬ 
zation  confers  simultaneous  expression  of  sim¬ 
ilar  miRNAs,  possibly  leading  to  combinatorial 
diversity  and  synergy  in  the  biological  effects. 
However,  all  the  miRNAs  from  a  single  tran¬ 
scriptional  cluster  are  not  expressed  at  equal 
levels,  suggesting  that  miRNAs  are  also  reg¬ 
ulated  posttranscriptionally.  Also,  a  significant 
portion  of  miRNAs  are  located  in  the  intronic 
region  of  protein-coding  or  -noncoding  tran¬ 
scription  units  (12),  whereas  a  minor  subset  of 
miRNAs  are  mapped  to  repetitive  sequences 
such  as  long  interspersed  nuclear  elements  (13). 
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Expression 

Many  miRNAs  exhibit  characteristic  expression 
patterns.  Some  miRNAs  are  differentially  ex¬ 
pressed  in  developmental  stages,  for  instance 
the  two  founding  members  lin-4  and  let-  7  in 
C.  elegans.  For  this  reason,  they  were  once 
known  as  small  temporal  RNAs  because  they 
are  expressed  in  specific  temporal  phases  of  de¬ 
velopment  and  regulate  developmental  timing. 
Many  miRNAs  are  expressed  in  a  tissue-specific 
manner.  In  the  few  cases  where  a  forward  ge¬ 
netic  approach  in  worms  or  flies  led  to  identifi¬ 
cation  of  a  miRNA,  the  function  of  the  miRNA 
was  inferred  from  the  phenotype  of  the  mu¬ 
tant.  In  mammalian  cells,  however,  a  miRNA 
expression  profile  is  usually  the  first  clue  to  its 
possible  role.  Analogous  to  mRNA  expression, 
miRNA  expression  is  determined  by  intrinsic 
cellular  factors  as  well  as  diverse  environmental 
variables. 

Expression  of  a  miRNA  can  be  measured  by 
molecular  biology  techniques,  such  as  North¬ 
ern  blot,  RNase  protection  assay,  and  primer 
extension  assay.  The  small  size  of  miRNAs 
initially  hampered  polymerase  chain  reaction 
(PCR)-based  methods.  However,  since  the  de¬ 
velopment  of  adaptor-mediated  quantitative 
real-time  PCR  (qRT-PCR;  reviewed  in  Refer¬ 
ence  14),  PCR-based  techniques  have  become 
very  popular  due  to  their  high  sensitivity.  Mi¬ 
croarray  techniques  (reviewed  in  Reference 
15)  are  widely  used  to  comprehensively  as¬ 
say  the  entire  miRNome  (the  global  miRNA 
expression  profile)  in  tissues  or  in  cell  lines 
(Table  1).  In  addition  to  microarray  and  qRT- 
PCR,  both  in  situ  hybridization  (16)  and  se¬ 
rial  analysis  of  gene  expression  (SAGE)  adapted 
for  small  RNAs  (17)  have  been  used  to  obtain 
miRNomes.  Interest  in  the  SAGE  approach  was 
stimulated  by  recent  innovations  in  ultrahigh- 
throughput  sequencing  that  provide  a  powerful 
tool  for  various  genomics  studies.  Overall,  these 
technical  improvements  are  expected  to  greatly 
widen  the  repertoire  of  miRNAs  in  a  variety  of 
biological  systems. 


Biogenesis 

Biogenesis  of  a  miRNA  begins  with  the  syn¬ 
thesis  of  a  long  transcript  known  as  a  pri- 
miRNA  (Figure  1).  In  general,  pri-miRNAs  are 
transcribed  by  RNA  polymerase  II  and  retain 
mRNA  features  such  as  5'  cap  structure  and 
y  poly(A)  tail  (18,  19).  However,  other  path¬ 
ways  generate  in  a  minor  set  of  miRNAs,  espe¬ 
cially  from  genomic  repeats.  For  example,  RNA 
polymerase  III  is  responsible  for  transcription 
of  miRNAs  in  Alu  repeats  (20). 

In  the  nucleus,  pri-miRNA  is  processed  to 
pre-miRNA  by  RNase  III  enzyme  Drosha  and 
its  interacting  partner  DGCR8  (DiGeorge  syn¬ 
drome  critical  region  gene  8)  (2 1-23).  DGCR8 
recognizes  the  stem  and  the  flanking  single- 
stranded  RNA  (ssRNA)  and  serves  as  a  ruler  for 
Drosha  to  cut  the  stem  approximately  1 1  nu¬ 
cleotides  away  from  the  stem-ssRNA  junction 
to  release  the  hairpin-shaped  pre-miRNA  (24). 
A  subset  of  miRNAs  known  as  miRtrons  bypass 
the  Drosha  requirement  by  taking  an  alterna¬ 
tive  pathway,  where  pre-miRNAs  are  derived 
as  a  byproduct  of  a  splicing  event  (25-27).  The 
pre-miRNA  is  exported  by  exportin-5  to  the  cy¬ 
toplasm  (28,  29)  and  is  subsequently  converted 
to  mature  duplex  miRNA  by  another  RNase 
III  enzyme,  Dicer  (30).  The  two  strands  of  the 
duplex  are  separated  by  a  RNA  helicase  (31) 
during  the  specific  assembly  of  a  ribonucleo- 
protein  complex  containing  miRNA  (miRNP). 
As  with  the  strand  selection  of  a  siRNA  (32,33), 
the  strand  whose  5'  end  forms  the  more  unsta¬ 
ble  duplex  with  its  partner  seems  to  preferen¬ 
tially  survive  as  the  miRNA  in  the  miRNP  (24, 
32,33).  However,  the  detailed  molecular  events 
that  take  place  during  the  miRNP  assembly  and 
strand  selection  are  not  yet  clear. 

Mechanisms  of  Action 

miRNAs  posttranscriptionally  suppress  the  tar¬ 
get  mRNA  expression,  mostly  through  interac¬ 
tion  with  the  V  UTR  (Figure  1).  However,  the 
exact  mechanism  for  miRNAs  on  their  targets 


qRT-PCR: 

quantitative  real-time 
polymerase  chain 
reaction 

miRNome:  the  global 
microRNA  expression 
profile 

pri-miRNA:  primary 
transcript  of  a 
microRNA 
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Table  1  Global  microRNA  expression  levels  (miRNome)  analyses  in  cancers 


Tumor  type 

Profiling 

methoda 

miRNAs  deregulated  in  tumors 

Notes 

Reference 

Upregulated 

Downregulated 

Breast  cancer 

Bead 

Classification  of  breast  tumors 
into  five  subtypes 

93 

qRT-PCR, 

array 

miR-2 1 

165 

Array 

Classification  according  to 
ErbB2/ER  status 

158 

miR-2 1 

miR-125b,  -145 

Clinicopathologic  features 

185 

Prostate  cancer 

Array 

let- 7c,  miR-125b, 
miR-145 

186 

158 

Pancreatic 

cancer 

Membrane 

Classification  according  to 
androgen  dependency 

66 

Array 

miR-2 1,  -221,  -222, 
-181a, -181b, -181d, 
-155 

Chronic  pancreatitis  and 
normal  tissues 

187 

miR-196a 

miR-2 17 

Pancreatic  ductal 

adenocarcinoma  versus 
chronic  pancreatitis  and 
normal  tissues 

160 

qRT-PCR  of 
pre-miRNA 

let- 7i,  miR-100,-155, 
-221,-301,-21, 

-181a, -125b,  -212, 
-376a 

miR-375 

Pancreatic  adenocarcinoma, 
paired  benign  tissue,  normal 
pancreas,  chronic  pancreatitis 

188 

Pancreatic 

tumor 

Array 

miR-103,-107, -23a, 
-26b, -342,  -192, 

-204,  -211,-21 

miR-155 

Endocrine  tumors  versus 

acinar  cell  carcinoma 

189 

Lung  cancer 

Array 

miR-2 1,-205 

miR-126* 

Clinicopathological  features 

190 

Ovarian  cancer 

Array 

miR-200a,  -141 

miR-199a,  -140, 
-145,  -125b-l 

Clinicopathological  features 

191 

Cervical  cancer 

Array 

Versus  Drosha  overexpression 

95 

Uterine 

leiomyoma 

Array 

let- 7,  miR-2 1, -23b 

Clinicopathological  features 
(tumor  size,  race,  etc.) 

192 

Hepatocellular 

cancer 

Array 

miR-2 1 

134 

miR-2  21 

let-1  miR-122a 

Versus  liver  cirrhosis 

193 

miR-224,  -18 

miR-199a,  -199a*, 
-200a 

194 

Thyroid  cancer 

Array 

miR-30d,  -125b, 
-26a,  and  -30a-5p 

Anaplastic  thyroid  carcinomas 
versus  normal  tissues 

67 

miR-197,  -346  in 

FTC  relative  to  FA 

FTC  versus  FA 

195 

miR-221,  -222, 

-181b 

Papillary  thyroid  carcinoma 
versus  normal  tissues 

159 

miR-221, -222, 

-146b 

Papillary  thyroid  carcinoma 
versus  normal  tissues 

100 

( Continued ) 
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Table  1  ( Continued ) 


Tumor  type 

Profiling 

method 

miRNAs  deregulated  in  tumors 

Notes 

Reference 

Upregulated 

Downregulated 

Colorectal 

cancer 

Array 

miR-25,  -92  in  MSS 
relative  to  MSI-H 

Microsatellite  stability  versus 
high  microsatellite  instability 

196 

qRT-PCR 

miR-31,  -96,  -135b, 
-183 

miR-133b,  -145 

miR-3 1  according  to  cancer 
stage 

197 

Cloning 

miSAGE 

17 

Pituitary 

adenomas 

Array 

miR-26a,  -149 

miR-21, -141, -144 

Clinicopathological  features 

198 

Neuroblastoma 

qRT-PCR 

Neuroblastoma  subtypes 

69 

Glioblastoma 

Array 

miR-221,  -10b 

miR-128, -181b, 
-181a 

199 

Leukemia 

Bead 

miR-128a,  -128b  in 
ALL  compared  with 
AML 

let-lb ,  miR-223  in 
ALL  compared 
with  AML 

ALL  versus  AML 

200 

Cloning 

miR-21,  -150,  -155 

miR-92, -222 

Various  CLLs 

201 

Array 

miR-29b, -181b  in 
aggressive  CLL 
with  llq  deletion 

Three  subtypes  of  CLL 

202 

miR-16-l,-15a 

CLL  according  to  ZAP-70  and 
IgVH 

203 

Leukemia 

Array 

B  cell  CLL 

65 

Various  tumors 

Array,  cloning 

Various  sarcomas 

204 

Array 

Solid  tumors  (lung,  breast, 
stomach,  prostate,  colon,  and 
pancreatic) 

122 

Bead 

Comprehensive  collection  of 
various  cancers  and  normal 

tissues 

64 

Profiling  methods  are  bead-based  flow  cytometric  profiling  (bead),  quantitative  real-time  polymerase  chain  reaction  (qRT-PCR),  microarray  of  miRNA 
chip  (array),  membrane  spotted  array  (membrane),  cloning,  and  serial  analysis  of  gene  expression  of  miRNA  (miSAGE).  miRNAs  of  notable  change  are 
indicated  by  asterisks. 

Abbreviations:  ALL,  acute  lymphoblastic  leukemia;  AMT,,  acute  myeloid  leukemia;  CLL,  chronic  lymphocytic  leukemia;  FA,  follicular  adenoma;  FTC, 
follicular  thyroid  carcinoma;  MSS,  microsatellite  stability;  MSI-H,  high  microsatellite  instability;  pre-miRNA,  precursor  miRNA. 


(reviewed  in  Reference  34)  is  still  controversial. 
In  contrast  to  the  perfect  sequence  complemen¬ 
tarity  between  siRNA  and  mRNA,  there  are 
mismatches  and  bulges  in  most  miRNA  target 
sites.  Comparison  between  siRNA  and  miRNA 
in  earlier  papers  suggested  that  siRNA  destabi¬ 
lizes  mRNA,  whereas  miRNA  inhibits  mRNA 
translation  without  affecting  the  mRNA  level. 
Therefore,  the  degree  of  complementarity  be¬ 


tween  short  RNA  and  the  target  was  thought  to 
be  a  major  determinant  distinguishing  the  two 
mechanisms. 

Although  the  translational  repression  mech¬ 
anism  still  holds  true  for  many  miRNAs,  it  has 
also  been  demonstrated  (35-38)  that  a  miRNA 
can  decrease  the  level  of  a  target  mRNA  de¬ 
spite  imperfect  sequence  complementarity  be¬ 
tween  the  miRNA  and  the  target.  mRNA 
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Strand  selection  and  RISC  assembly 


i  Dicer,  TRBP 
Ago,  PACT 


Translational  inhibition 

(initiation/elongation  step) 


Cytoplasmic  P-body 

(RNA  degradation) 


Canonical  pathway 
Minor  pathway 


Figure  1 

Pathways  of  microRNA  (miRNA)  biogenesis  and  action.  The  canonical  pathway  {red  arrows)  and  minor  pathways  {light  blue  arrows)  are 
depicted.  Abbreviations:  miRNP,  ribonucleoprotein  complex  containing  miRNA;  ORF,  open  reading  frame;  P-body,  mRNA- 
processing  body;  PACT,  interferon-inducible  double-stranded  RNA-dependent  activator;  RISC,  RNA-induced  silencing  complex; 
TRBP,  human  immunodeficiency  virus-transactivating  RNA-binding  protein;  UTR,  untranslated  region  of  mRNA. 
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degradation  by  a  miRNA,  which  is  distin¬ 
guished  from  siRNA-mediated  mRNA  cleav¬ 
age,  can  be  explained  by  mRNA-processing 
bodies  (P-bodies),  which  are  sites  for  RNA  de¬ 
cay  (39,  40).  Plausibly,  miRNAs  inhibit  trans¬ 
lation  of  target  mRNAs,  which  are  then  se¬ 
questered  to  P-bodies  and  become  subject  to 
degradation.  This  model  fits  well  for  both 
mechanisms,  although  miRNAs  have  also  been 
proposed  to  lead  to  the  degradation  of  the  tar¬ 
get  mRNAs  without  sequestration  to  P-bodies 
(41,  42).  In  some  cases  where  mRNA  transla¬ 
tion  is  inhibited  by  a  miRNA  but  where  mRNA 
levels  remain  similar,  sequestration  to  P-bodies 
might  be  inefficient.  Nonetheless,  it  is  possible 
that  different  mechanisms  apply  to  individual 
miRNA-mRNA  interactions. 

Target  Prediction  and  Identification 

Prediction  of  miRNA  targets  (reviewed  in 
Reference  46)  is  important,  given  that  miRNAs 
exert  their  function  by  regulating  target 
mRNAs.  The  specificity  of  miRNA-mRNA  in¬ 
teraction  is  mainly  conferred  by  the  first  eight 
nucleotides  of  a  miRNA  (known  as  a  seed  se¬ 
quence)  (43).  The  likelihood  that  a  predicted 
target  is  a  bona  fide  target  is  influenced  not 
only  by  seed  pairing  but  also  by  other  factors 
such  as  the  number  of  target  sites,  the  context 
of  surrounding  sequence  in  mRNA  (44),  and  the 
occlusion  of  target  sites  by  RNA-binding  pro¬ 
teins  (45).  Currently,  several  computational  al¬ 
gorithms  (reviewed  in  Reference  46)  can  predict 
the  target  mRNA(s),  but  they  are  far  from  per¬ 
fect.  The  gold  standard  is  experimental  demon¬ 
stration  that  (a)  a  luciferase  reporter  fused  to  the 
y  UTR  of  the  predicted  target  is  repressed  by 
overexpression  of  the  miRNA  and  (b)  this  re¬ 
pression  is  abrogated  by  point  mutation  in  the 
target  sequence(s)  in  the  V  UTR.  Many  targets 
are  predicted  by  in  silico  analyses,  but  not  all  of 
them  are  confirmed  as  real  targets  in  this  bio¬ 
logical  assay. 

The  in  silico  predictions  have  been  com¬ 
plemented  by  experimental  screenings  for  tar¬ 
gets  of  miRNAs.  Because  miRNAs  destabilize 


mRNAs,  miRNA  depletion  by  a  knockdown 
of  miRNA-processing  machinery  is  expected 
to  result  in  upregulation  of  target  mRNA, 
which  can  be  assayed  by  microarray  analy¬ 
sis.  This  approach  identified  HMGA2  (high- 
mobility  group  AT-hook  2)  as  a  target  of  let- 7 
(38).  Microarrays  can  also  identify  mRNAs  that 
are  decreased  upon  transfection  of  a  miRNA, 
and  screening  these  mRNAs  for  computation¬ 
ally  predicted  sites  targeted  by  the  miRNA 
yields  some  bona  fide  targets  (37).  Aproteomics 
approach  to  identify  proteins  decreased  by 
miRNAs  led  to  the  identification  of  tropomyosin 
1  as  a  miR-21  target  (47).  Another  screen  em¬ 
ployed  a  library  of  miRNA-expressing  vectors 
and  a  sensor  plasmid  containing  green  fluores¬ 
cent  protein  fused  with  the  3'  UTR  of  the  gene 
of  interest.  p27(Kipl)  was  found  to  be  targeted 
by  miR-22 1  and  -222  by  this  approach  (48).  An¬ 
other  approach  is  to  isolate  the  miRNA-target 
mRNA  complex  by  immunoprecipitation  of  a 
component  of  RNA-induced  silencing  complex 
(RISC)  followed  by  micro  array  hybridization  of 
the  precipitated  mRNAs  (49). 

Although  a  growing  number  of  miRNA- 
target  pairs  are  being  identified  (Table  2),  the 
fraction  of  validated  pairs  is  still  small,  given 
that  30%  of  mRNAs  have  been  predicted  to 
be  miRNA  targets  (50,  51).  New  target  iden¬ 
tification  will  be  facilitated  by  the  accumula¬ 
tion  of  validated  miRNA-target  pairs,  which 
can  be  utilized  to  develop  a  better  prediction 
algorithm.  In  addition,  experimental  tools  for 
high-throughput  target  screening  need  to  be 
improved. 

Biological  Roles 

miRNAs  play  important  roles  in  cell  fate  de¬ 
termination,  proliferation,  and  cell  death.  In 
addition  to  these  vital  processes,  miRNAs  are 
implicated  in  diverse  cellular  activities,  such  as 
immune  response  (reviewed  in  References  52 
and  53),  insulin  secretion  (54),  neurotransmit¬ 
ter  synthesis  (55),  circadian  rhythm  (56),  and 
viral  replication  (57).  This  list  will  undoubtedly 
expand  as  experimental  data  accumulate. 


HMGA2:  high- 
mobility  group 
AT-hook  2 
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Table  2  Oncogenic  or  tumor-suppressive  microRNAs  (miRNAs)  and  their  direct  target  genes 


miRNA 

Target  gene 

Note 

Reference 

let-7 

RAS 

Lung  cancer 

98 

let-7 

Ccnd2,  Cdk6,  Cdc25A 

Cancer  cell  lines  (A549,  HepG2,  HeLa) 

59 

let-7 

HMGA2 

Ovarian  cancer 

132 

let-7 c 

c-Myc 

Liver  tumors 

83 

let-7 

HMGA2 

Lung  cancer  cell  lines 

38 

kt-7g 

c-Myc,  k-RAS 

70 

let-7 

HMGA2 

106 

let-7 

NF2 

Cholangiocarcinoma  cell  lines 

206 

miR-9,  -125a, -125b 

trkC 

Neuroblastoma 

107 

miR- 1  Ob 

HOXDIO 

Breast  cancer 

146 

miR-16-l,-15a 

Bcl2 

Chronic  lymphocytic  leukemia 

140 

miR-17-5p 

AIB1 

Breast  cancer 

162 

miR-17-5p, -20 

T|3RII 

74 

miR- 1 8 

CTGF 

Colon  cancer  model  of  angiogenesis 

147 

miR- 19 

THBS1 

Colon  cancer  model  of  angiogenesis 

147 

miR-20a 

E2F1,  2,  3 

126 

miR-2 1 

Pdcd4 

Colorectal  cancer 

145 

miR-21 

PTEN 

HCC 

134 

miR-2 1 

TPM1 

47 

miR-21 

PTEN 

Cholangiocarcinoma 

176 

miR-2  7b 

CYP1B1 

Breast  cancer 

207 

miR-2 9a,  b,  c 

DNMT-3A,  DNMT-3B 

Non-small  cell  lung  cancer 

82 

miR-2  9s 

Mcl-1 

Cholangiocarcinoma  cell  line 

84 

miR-2  9b,  -181b 

TCL1 

Chronic  lymphocytic  leukemia 

202 

miR- 3  4a 

E2F3 

Neuroblastoma 

76 

miR- 3  4a 

Ccne2,  MET 

62 

miR-34a,  -34b,  -34c 

Bcl2 

Non-small  cell  lung  cancer 

60 

miR- 106a  cluster 

Mylip,  Hipk3,  Rbpl-like 

T  cell  leukemia 

116 

miR- 122a 

Cyclin  G1 

HCC  versus  liver  cirrhosis 

193 

miR- 124a 

Cdk6 

Colon  cancer,  lung  cancer 

80 

miR-125a,  -125b 

ERBB2  and  ERBB3 

Breast  cancer 

144 

miR- 12 5b 

Bak-1 

Prostate  cancer 

205 

miR- 127 

Bcl6 

Bladder  (and  prostate)  cancer 

208 

miR-2  06 

ERa 

Breast  cancer 

105 

miR-22 1,-222 

p27  (Kip) 

Glioblastoma 

135 

miR-22 1,-222 

p27  (Kip) 

Glioblastoma 

48 

miR-22 1,-222 

P27  (Kip) 

Prostate  cancer 

136 

miR-372,  -373 

LATS2 

Testicular  germ  cell  tumor 

209 

miR-378 

Sufu,  Fus-1 

U87  glioblastoma  cell  line 

148 

(' Continued ) 
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Table  2  {Continued) 


miRNA 

Target  gene 

Note 

Reference 

BART  16,  17-5p,  l-5p 
(EBV-encoded  miRNAs) 

LMP1  (EBV-encoded  oncogene) 

Nasopharyngeal  carcinoma 

111 

KSHV  miRNAs  miR-K12-l, 
-K12-3-3p,  -K12-6-3p,  -K12-11 

THBS1 

112 

Listed  for  each  miRNA  is  its  target  gene(s).  A  directly  regulatory  connection  was  demonstrated  between  each  miRNA  and  target  a  by  a  luciferase  assay 
using  the  3'  untranslated  region  (UTR)  of  the  mRNA. 

Abbreviations:  AIB1,  amplified  in  breast  cancer  1;  BART,  BamW  1  A  rightward  transcript;  Ccnd2,  cyclin  D2;  Ccne2,  cyclin  E2;  CTGF,  connective  tissue 
growth  factor;  DNMT,  DNA  methyltransferase;  EBV,  Epstein-Barr  virus;  ERBB,  v-erb-b2  erythroblastic  leukemia  viral  oncogene  homolog;  E2F,  E2F 
transcription  factor;  HCC,  hepatocellular  cancer;  Hipk3,  homeodomain-interacting  protein  kinase  3;  HMGA,  high-mobility  group  AT-hook  2;  KSHV, 
Kaposi  sarcoma-associated  herpesvirus;  LATS2,  large  tumor  suppressor,  homolog  2;  MET,  hepatocyte  growth  factor  receptor;  Mylip,  myosin 
regulatory  light  chain-interacting  protein;  PTEN,  phosphatase  and  tensin  homolog;  Rbpl-like,  retinoblastoma-binding  protein  1-like;  THBS1, 
thrombospondin  1;  TPM1,  tropomyosin  1;  trkC,  tropomyosin-related  kinase  C;  T|3RII,  transforming  growth  factor  beta-receptor  type  II. 


miRNA  transfection  results  in  up-  and 
downregulation  of  a  number  of  mRNAs.  In 
several  cases,  a  set  of  genes  belonging  to  a 
particular  cellular  pathway  are  enriched  in  the 
genes  whose  expression  was  altered  (35,  58- 
62).  The  length  of  the  miRNA  seed  sequence  is 
comparable  to  that  of  the  consensus  sequence 
of  transcription  factor-binding  elements.  Thus, 
although  miRNAs  are  expected  to  directly  reg¬ 
ulate  a  large  set  of  genes  simultaneously,  ran¬ 
dom  chance  is  unlikely  to  explain  the  enrich¬ 
ment  of  genes  in  a  particular  pathway  in  the  list 
of  putative  targets. 

MicroRNA  AND  CANCER 

Profiling  of  MicroRNA  Expression 
in  Cancers 

miRNAs  are  aberrantly  expressed  in  a  variety 
of  cancers.  This  was  first  observed  in  miR-15a 
and  -16-1,  which  are  clustered  at  chromosome 
13ql4,  a  frequently  deleted  region  in  B  cell 
chronic  lymphocytic  leukemia  (CLL)  and  other 
cancers  (63).  Concomitantly,  reduction  of  these 
two  miRNAs  was  observed  in  the  cancer  sam¬ 
ples  relative  to  the  normal  tissues. 

miRNome  analyses  have  become  easier 
with  the  introduction  of  microarray  techniques 
(Table  1).  Many  miRNAs  are  found  to  be 
up-  or  downregulated  in  the  cancer  samples 


relative  to  the  normal  tissue  counterparts 
(Supplemental  Table  1;  follow  the 
Supplemental  Material  link  from  the 
Annual  Reviews  home  page  at  http://www. 
annualreviews.org).  It  is  beyond  the  capacity 
of  this  review  to  describe  all  the  miRNAs 
that  are  changed  in  cancers  from  profiling 
data.  Instead,  in  Table  1  we  summarize  the 
literature  reporting  miRNomes  in  tumors. 
If  a  change  in  miRNA  expression  has  been 
corroborated  by  a  conventional  method  (e.g., 
Northern  blot,  RNase  protection  assay,  or 
qRT-PCR),  it  is  indicated  in  Supplemental 
Table  1. 

In  addition  to  the  distinction  of  tumors  from 
normal  tissue,  miRNA  expression  is  character¬ 
istic  for  a  cancer  type,  stage,  and  other  clin¬ 
ical  variables.  The  first  systematic  analyses  of 
hundreds  of  cancer  samples  and  normal  tissues 
(64)  successfully  classified  various  cancers  based 
on  the  miRNome.  Surprisingly,  the  miRNome 
was  better  at  predicting  cancer  type  and  stage 
than  the  mRNA  expression  profile;  therefore, 
the  miRNome  has  been  proposed  as  a  useful 
tool  for  cancer  diagnosis  and  prognosis.  For 
example,  clustering  of  CLL  samples  according 
to  the  miRNome  revealed  miRNA  signatures 
that  correlate  with  overexpression  of  ZAP-70,  a 
predictor  of  early  disease  progression  (65).  The 
utility  of  miRNA  in  diagnosis  is  discussed  fur¬ 
ther  below. 
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Although  some  miRNAs  are  increased,  most 
are  repressed  in  cancers  relative  to  normal  tis¬ 
sue  counterparts  (64, 66-69).  In  agreement  with 
these  observations,  global  depletion  of  miRNAs 
by  knockdown  of  the  miRNA-processing  ma¬ 
chinery  stimulates  cell  transformation  and  tu- 
morigenesis  in  vivo  (70).  This  implies  that  the 
miRNA  alteration  is  not  simply  an  end  result 
of  tumorigenesis  but  that  it  actively  contributes 
to  cancer  development.  Despite  the  general  re¬ 
duction  of  miRNAs  in  cancers,  several  miRNAs 
are  upregulated,  some  of  which  undoubtedly 
play  oncogenic  roles. 


Mechanism  for  MicroRNA 
Deregulation  in  Cancers 

As  described  in  the  previous  section,  many 
miRNAs  are  deregulated  in  cancers.  Aberrant 
expression  of  miRNAs  can  arise  through  a  num¬ 
ber  of  different  mechanisms  (Figure  2). 

Genomic  abnormalities.  As  exemplified  by 
miR-15a  and  -16-1,  chromosomal  abnormal¬ 
ity  is  one  reason  for  miRNA  deregulation  in 
cancers.  Tumorigenesis  is  often  accompanied 
by  chromosomal  aberrations  such  as  deletion, 
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Figure  2 

Various  mechanisms  for  deregulation  of  microRNA  (miRNA)  expression.  Abbreviations:  miRNP, 
ribonucleoprotein  complex  containing  miRNA;  ORF,  open  reading  frame;  SNP,  single  nucleotide 
polymorphism;  UTR,  untranslated  region. 
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amplification,  and  translocation.  In  silico  anal¬ 
ysis  has  revealed  that  a  significant  fraction  of 
miRNAs  are  mapped  to  these  cancer-associated 
genomic  regions  or  fragile  sites  in  human  (71) 
and  in  mouse  (72).  In  many  cases,  miRNA  lev¬ 
els  are  correlated  with  changes  in  copy  num¬ 
ber  of  the  genomic  loci,  as  demonstrated  by 
comparison  between  array  comparative  ge¬ 
nomic  hybridization  data  and  miRNA  expres¬ 
sion  data  (66,  73-78). 

Epigenetic  factors.  Epigenetic  factors  can 
also  affect  miRNA  expression.  In  many  cancers, 
hypermethylation  of  CpG  islands  in  promoter 
regions  results  in  heritable  transcriptional 
silencing  of  tumor-suppressor  genes.  Gene 
silencing  by  DNA  methylation  is  closely  re¬ 
lated  to  histone  modification.  In  silico  anal¬ 
yses  have  indicated  CpG  islands  near  dozens 
of  miRNAs  (79).  In  addition,  some  miRNAs 
are  upregulated  ( a )  upon  the  exposure  of 
cells  to  the  demethylating  agent  5-aza-2'- 
deoxycytidine  (79),  ( b )  upon  mutation  of  DNA 
methyltransferases  (DNMTs)  (80),  and  (c)  upon 
histone  deacetylase  inhibitor  treatment  (81). 
These  studies  have  identified  some  miRNAs 
that  are  repressed  by  CpG  hypermethylation  in 
cancers  relative  to  normal  tissue.  Representa¬ 
tive  examples  include  miR-9-1  in  breast  cancer 


(79)  and  miR-124a  in  colorectal  tumors  (80).  In 
the  case  of  miR-124a,  hypermethylation  is  tu¬ 
mor  type  specific,  as  no  methylation  is  seen  in 
neuroblastoma.  Also,  epigenetic  silencing  of  a 
miRNA  may  be  a  reflection  of  tissue  specificity. 
For  example,  miR-124a  is  normally  highly  ex¬ 
pressed  in  neuronal  tissues,  so  its  epigenetic  re¬ 
pression  in  colorectal  tumors  is  not  surprising. 

miRNAs  may  counteract  CpG  methylation. 
For  example,  miR-29  directly  targets  DNMT- 
3 A  and  -3B.  In  agreement  with  this  observa¬ 
tion,  ectopic  expression  of  miR-29  results  in 
a  global  reduction  of  DNA  methylation,  sub¬ 
sequently  leading  to  a  depression  of  tumor- 
suppressor  genes  that  had  been  silenced  by  pro¬ 
moter  methylation  in  cancer  cells  (82). 

Transcriptional  regulation.  Transcription 
factors  may  induce  miRNAs  by  activating  the 
transcription  of  pri-miRNAs.  This  mechanism 
has  been  well  documented  in  several  cases 
where  tissue-specific  miRNAs  are  turned  on 
by  transcription  factors  during  differentiation. 
Given  the  wide  impact  of  transcription  factors 
on  fundamental  cellular  processes,  it  is  not 
surprising  that  many  oncogenes  or  tumor 
suppressors  are  transcription  factors.  Many 
miRNA-transcription  factor  relationships  have 
been  discovered  in  cancers  (listed  in  Table  3). 


Table  3  MicroRNAs  (miRNAs)  regulated  by  transcription  factors 


Transcription  factor 

Target  miRNA 

Note 

Reference 

AR 

miR-125b 

Prostate  cancer 

205 

MYCN 

miR-17-5p, -92,-320 

Neuroblastoma 

124 

Twist 

miR-lOb 

Breast  cancer 

146 

p53 

miR-34a,  miR-34b-34c  cluster 

Cell  lines 

60 

P53 

miR-34a 

Cell  lines 

137 

p53 

miR-34a,  miR-34b-34c  cluster 

Cell  lines  and  in  vivo 

62 

Stat3 

miR-21 

Myeloma  cell  lines 

91 

HIF 

miR-26,  -210 

150 

E2F3 

miR- 17-92  cluster 

125 

E2F1,  2,  3 

miR- 17-92  cluster 

126 

Myc 

miR- 17-92  cluster 

Angiogenesis  in  colon  cancer  model 

147 

Listed  are  miRNAs  whose  promoter  regions  were  analyzed  by  a  luciferase  reporter  assay  and  were  shown  to  bind  to  the  transcription  factor  by  chromatin 
immunoprecipitation  assay. 

Abbreviations:  AR,  androgen  receptor;  E2F,  E2F  transcription  factor;  HIF,  hypoxia-inducing  factor;  MYCN,  Myc  myelocytomatosis  viral-related 
oncogene,  neuroblastoma  derived;  Stat3,  signal  transducer  and  activator  of  transcription  3. 


www.annualreviews.org  •  MicroRNAs  in  Cancer  209 


Annu.  Rev.  Pathol.  Mech.  Dis.  2009.4:199-227.  Downloaded  from  arjournals.annualreviews.org 
by  University  of  Virginia  Libraries  on  02/05/09.  For  personal  use  only. 


We  discuss  three  of  these — p53,  c-Myc,  and 
E2F — in  detail  below. 

Regulation  at  microRNA  processing  steps. 

The  steady-state  level  of  a  mature  miRNA  is 
determined  not  only  by  the  transcription  rate 
of  pri-miRNA  but  also  by  the  processing  ef¬ 
ficiency  of  its  precursors  and  by  its  stability. 
miRNAs  often  exhibit  a  discrepancy  in  expres¬ 
sion  of  the  mature  form  relative  to  that  of 
the  precursor  (83-88).  Although  miRNAs  in  a 
genomic  cluster  are  usually  expressed  from  a 
common  pri-miRNA,  the  levels  of  individual 
miRNAs  in  the  cluster  are  not  necessarily  co¬ 
ordinated  (89,  90).  A  time-course  experiment 
after  induction  of  pri-miR-2 1  revealed  delayed 
kinetics  in  accumulation  of  mature  miR-2 1  (9 1). 
Collectively,  these  observations  indicate  that 
miRNA  processing  and  stability  are  impor¬ 
tant  factors  that  determine  miRNA  expression 
level. 

This  mechanism  has  been  confirmed  in  can¬ 
cers  by  a  comprehensive  analysis  of  expression 
data  (92).  As  stated  above,  miRNAs  are  gen¬ 
erally  reduced  in  cancers  relative  to  the  nor¬ 
mal  tissues.  When  a  miRNA  resides  within  a 
gene,  the  host  gene  can  be  regarded  as  the  pri- 
miRNA.  Comparison  of  microarray  data  be¬ 
tween  mRNA  and  miRNA  revealed  that  the 
miRNA  reduction  in  cancers  is  poorly  corre¬ 
lated  with  a  reduction  of  the  host  gene  expres¬ 
sion.  Aberrant  expression  of  miRNAs  during 
tumorigenesis  is  presumably  often  due  to  alter¬ 
ations  at  posttranscriptional  steps. 

The  expression  levels  of  Dicer  or  Drosha 
are  altered  in  a  number  of  cancers  (70,  93- 
96).  Drosha  upregulation  is  seen  in  more  than 
half  of  cervical  squamous  cell  carcinoma  (SCC) 
specimens  and  is  likely  due  to  the  copy  num¬ 
ber  gain  at  chromosome  5p,  where  the  Drosha 
gene  is  located  (95).  Hierarchical  clustering  of 
miRNA  expression  data  successfully  classified 
cervical  SCC  samples  into  two  groups  accord¬ 
ing  to  Drosha  overexpression.  Notably,  some 
miRNAs  were  reduced  upon  Drosha  overex¬ 
pression,  indicating  that  individual  miRNAs 
respond  differentially  to  an  elevation  of  the 
miRNA  processing  machinery.  Interestingly, 


Drosha  was  reported  to  interact  with  an  onco¬ 
genic  fusion  protein  derived  from  a  chromo¬ 
somal  translocation  in  some  leukemias  (97). 
This  interaction  affects  pri-miRNA  selection 
of  Drosha  and,  as  a  result,  influences  miRNA 
expression  patterns. 

Consequence  of  Aberrant  MicroRNA 
Expression  in  Cancers 

miRNAs  regulate  the  expression  of  their  tar¬ 
get  mRNAs,  so  over-  or  underexpression  of 
miRNAs  is  expected  to  result  in  down-  or  up¬ 
regulation,  respectively,  of  the  protein  prod¬ 
uct  of  the  target  mRNAs.  It  is  easy  to  impli¬ 
cate  a  miRNA  in  a  cancer  if  a  direct  target  of  a 
miRNA  is  an  oncogene  or  a  tumor  suppressor. 
Since  the  publication  of  a  study  showing  that 
let-  7  miRNA  directly  regulates  RAS  oncogene 
(98),  a  number  of  other  miRNA-target  pairs 
have  been  studied.  However,  the  number  of  ex¬ 
perimentally  validated  pairs  is  still  small  rela¬ 
tive  to  (a)  the  number  of  miRNAs  shown  to  be 
aberrantly  expressed  in  tumors  and  (b)  the  num¬ 
ber  of  in  silico  predicted  pairs,  mostly  because 
current  target  prediction  algorithms  are  not 
very  accurate.  A  current  list  of  miRNA-target 
pairs  implicated  in  various  cancers  is  given  in 
Table  2,  and  important  examples  are  discussed 
below. 

If  we  compare  global  gene  expression  pro¬ 
files  in  cancers  and  normal  tissues,  we  find  that 
many  miRNAs  and  mRNAs  are  deregulated. 
Therefore,  it  is  plausible  that  tumorigenesis  (or 
at  least  the  progression  of  a  cancer)  results  from 
changes  in  the  entire  miRNome,  rather  than 
from  the  change  of  a  single  miRNA  that  regu¬ 
lates  an  oncogenic  (or  tumor-suppressive)  tar¬ 
get  gene. 

Some  miRNAs  appear  to  be  deregulated 
in  cancers  much  more  frequently  than  oth¬ 
ers  (Supplemental  Table  1).  These  miRNAs 
may  play  key  roles  during  tumorigenesis.  For 
example,  the  miR- 17-92  cluster  and  miR-155 
have  been  shown  experimentally  to  be  bona  fide 
oncogenes,  as  their  ectopic  expression  acceler¬ 
ates  tumor  development.  These  two  miRNAs 
are  discussed  in  detail  in  separate  sections. 
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Aberrant  Action  of  MicroRNAs 
Without  Alteration  of  MicroRNA 
Expression  Level 

The  function  of  protein-coding  genes  is 
rendered  abnormal  by  point  mutations,  which 
either  transform  proto-oncogenes  to  onco¬ 
genes  or  abrogate  functions  of  tumor- 
suppressor  genes.  In  theory,  the  same  mech¬ 
anism  of  activation/inactivation  may  apply 
to  miRNAs.  However,  mutation  in  mature 
miRNA  sequence  seems  to  be  a  rare  event  (72, 
99-101),  presumably  due  to  the  small  size  of 
the  mature  miRNA.  In  addition,  miRNAs  are 
thought  to  be  more  tolerant  of  point  mutations 
or  of  single  nucleotide  insertion/deletion 
mutations  than  of  protein-coding  genes. 

In  contrast,  sequence  variation  in  miRNA 
target  sites  may  occur  and  may  play  a  role  in 
cancer.  In  silico  analyses  of  expressed  sequence 
tag  and  single  nucleotide  polymorphism  (SNP) 
databases  indicate  different  allele  frequencies  of 
miRNA-binding  sites  in  cancers  versus  normal 
tissues  (102).  Several  experiments  have  shown 
that  sequence  polymorphisms  in  miRNA  tar¬ 
get  sites  affect  miRNA-mRNA  interaction  and 
are  implicated  in  disease  phenotypes  (100, 103— 
105)  (Figure  2). 

An  interesting  illustration  of  this  mechanism 
can  be  found  in  let-  7  and  its  direct  target  onco¬ 
gene,  HMGA2  (38,  106).  Chromosomal  rear¬ 
rangements  at  the  HMGA2  locus  in  several  tu¬ 
mors  separate  the  open  reading  frame  (ORF) 
from  the  3'  UTR  that  contains  let- 7  target  sites. 
As  a  result,  HMGA2  escapes  from  suppression 
by  let-7,  is  overexpressed,  and  promotes  tumori- 
genesis  (Figure  2). 

An  alternative  splicing  event  may  result  in  a 
different  3'  UTR  that  displays  different  miRNA 
target  sites,  as  exemplified  in  the  targeting  of 
tropomyosin-related  kinase  C  by  miR-9,  -125a, 
and  -125b.  One  mRNA isoform  encodes  a  trun¬ 
cated  ORF  that  is  functionally  dominant  neg¬ 
ative  to  the  intact  protein.  In  this  isoform, 
the  y  UTR  contains  the  target  sites  of  these 
miRNAs.  In  contrast,  the  target  sites  are  absent 
in  another  isoform  encoding  the  intact  ORF; 
only  the  former  isoform  was  repressed  by  the 


miRNAs  (107).  Although  the  stop  codon  is  usu¬ 
ally  located  in  the  last  exon,  generation  of  differ¬ 
ent  y  UTRs  by  multiple  polyadenylation  sites 
or  alternative  splicing  has  been  known  to  oc¬ 
cur  in  a  small  but  significant  fraction  of  genes 
(108).  Thus,  variation  of  3'  UTR  and  of  atten¬ 
dant  miRNA  target  sites  is  expected  to  be  a 
mechanism  for  oncogene  activation  or  tumor- 
suppressor  inactivation  (Figure  2). 

Recently,  Steitz  and  colleagues  (109)  re¬ 
ported  that  miRNAs  activate  the  translation  of 
the  target  mRNA  in  cells  arrested  at  the  Go/Gi 
stage.  In  addition  to  aberrant  miRNA  expres¬ 
sion,  the  switch  from  repression  to  activation 
should  be  considered  in  studying  the  role  of 
miRNAs  in  differentiation  or  tumorigenesis,  as 
the  same  miRNA  may  exert  opposite  effects  in 
resting  cells  and  proliferating  cancer  cells  in  a 
given  tissue. 

Viruses  and  MicroRNAs 
in  Tumorigenesis 

Several  viruses  regulate  tumorigenesis  by 
expressing  viral  oncogenes  or  by  activating 
cellular  oncogenes  through  integration  of  vi¬ 
ral  DNA  into  genomic  loci.  Both  mecha¬ 
nisms  are  applicable  to  oncogenic  miRNAs.  For 
example,  Epstein-Barr  virus  (EBV)-encoded 
miRNAs  (110)  directly  target  a  viral  onco¬ 
gene,  LMP1,  whose  overexpression  is  dele¬ 
terious  to  host  cells.  EBV-encoded  miRNAs 
are  thought  to  enhance  EBV-mediated  cellu¬ 
lar  transformation  by  adjusting  LMP1  level 
to  a  sublethal  dose  (111).  Another  study 
showed  that  a  viral  miRNA  repressed  a  host 
tumor-suppressor  gene.  In  this  study,  miRNAs 
from  Kaposi  sarcoma-associated  herpesvirus 
(KSHV)  directly  targeted  an  antiangiogenic 
factor,  thrombospondin- 1  ( Tsp-1 )  (112).  Thus  vi¬ 
ral  miRNAs  act  as  viral  oncogenes. 

Viral  integration  near  miRNA  loci  may  lead 
to  aberrant  expression  of  miRNAs.  Indeed,  the 
BIC  locus  harboring  miR-155  was  originally 
described  as  a  frequent  integration  site  in  vi- 
rally  induced  lymphomas  (113).  Viral  integra¬ 
tion  sites  are  often  mapped  to  miRNA  loci 
(114)  such  as  the  miR- 17-92  cluster  (115),  the 


KSHV:  Kaposi 

sarcoma-associated 
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miR-106a-363  cluster  (116),  the  miR-29a-29b- 
1  cluster  (117),  and  the  miR-106a  cluster  (118). 

Various  Cancer  Types 

Comprehensive  miRNA  profiles  have  been 
reported  in  clinical  specimens  from  various 
cancers  (Table  1).  Consistent  with  the  notion 
that  miRNAs  are  expressed  in  a  tissue-specific 
manner,  miRNAs  differ  among  cancers  from 
various  tissue  origins.  Also,  a  few  miRNAs  ap¬ 
pear  to  be  frequently  deregulated  in  many  can¬ 
cers  (Supplemental  Table  1),  suggesting  that 
these  miRNAs  regulate  fundamental  processes 
such  as  cell  proliferation  and  apoptosis. 

Oncogenic  or  Tumor-Suppressive 
MicroRNAs 

miR- 17-92  cluster.  The  miR- 17-92  cluster, 
located  at  chromosome  1 3q3 1.3  in  humans,  is 
composed  of  six  miRNAs  (mir-17,  -18a,  -19a, 
-20a,  -19b-l,  and  -92a-l).  There  exists  a  similar 
cluster  at  chromosome  10  known  as  miR- 
106a-363,  which  also  contains  six  miRNAs 
(mir-106a,  -18b,  -20b,  -19b-2,  -92a-2,  and 
-363).  The  miR- 17-92  cluster  contains  the 


first  miRNAs  demonstrated  to  be  oncogenic. 
The  development  of  B  cell  lymphoma  is  sig¬ 
nificantly  accelerated  by  the  forced  expression 
of  the  miR- 17-92  cluster  in  transgenic  mice 
overexpressing  the  c-myc  oncogene  (119). 

Consistent  with  its  oncogenic  role,  the  miR- 
17-92  cluster  is  upregulated  in  a  variety  of  can¬ 
cers  including  lymphomas  (74,  75,  119,  120), 
lung  cancers  (77,  121),  and  others  (122).  There 
appear  to  be  two  mechanisms  for  upregulation 
of  this  cluster  in  cancers  (Figure  3).  One  is  the 
amplification  of  the  chromosome  13q31  locus 
in  several  lymphomas  and  other  cancers  (74,  77, 
78,  120).  The  other  is  transcriptional  activation 
of  the  pri-miRNA,  wherein  an  oncogenic  tran¬ 
scription  factor  c-Myc  binds  the  genomic  lo¬ 
cus  upstream  of  the  miR- 17-92  cluster  and  ac¬ 
tivates  its  expression  (123).  In  neuroblastoma 
cells,  MYCN,  a  protein  highly  homologous  to 
c-Myc,  appears  to  activate  this  cluster  instead 
of  c-Myc  (124).  E2Fs  also  activate  this  cluster 
(125,  126)  (Table  3). 

Given  that  the  E2Fs  are  direct  targets  of 
miR- 1 7  and  -2  0a,  the  miR- 1 7 -92  cluster  consti¬ 
tutes  a  complex  regulatory  network  with  c-Myc 
and  E2Fs.  E2F1  and  c-Myc  are  known  to  acti¬ 
vate  each  other  to  form  a  positive  feedback  loop 


chrl  3q31 .3  amplification 


Figure  3 

The  miR- 17-92  pathway. 
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(Figure  3).  All  of  the  E2Fs  (E2F1,  -2,  and  -3) 
are  capable  of  activating  the  miR- 17-92  clus¬ 
ter  and  are  subject  to  repression  by  miR- 17 
and  -20a  from  this  cluster  (Table  2).  However, 
degrees  of  activation  and  repression  vary  be¬ 
tween  individual  pairs  of  miRNA  and  E2F  fam¬ 
ily  members.  In  general,  E2Fs  are  thought  to  be 
proproliferative,  but  E2F1  is  also  proapoptotic. 
Although  this  network  is  too  complicated  to  de¬ 
scribe  in  a  straightforward  manner,  the  miR- 17- 
92  cluster  is  clearly  responsible  for  fine-tuning 
the  regulatory  network.  Upon  a  proliferative 
cue,  c-Myc  and  E2Fs  turn  on  the  miR- 17-92 
cluster.  This  cluster  represses  E2Fs,  thereby 
preventing  the  uncontrolled  amplification  of 
the  positive  feedback  loop  between  E2Fs  and 
c-Myc.  In  addition,  repression  of  E2F1  by  the 
miRNAs  may  help  minimize  the  proapoptotic 
potential  of  E2F1.  This  network  is  depicted  in 
Figure  3. 

miR-21.  miR-21  is  upregulated  in  almost  all 
kinds  of  cancers  (Supplemental  Table  1).  miR- 
21  is  transcriptionally  activated  by  signal  trans¬ 
ducer  and  activator  of  transcription  3  (Stat3) 
in  the  interleukin-6  signaling  pathway  (91) 
(Table  3).  miR-2 1  ’s  role  in  invasion  and  metas¬ 
tasis  is  well  characterized  and  is  discussed  fur¬ 
ther  below. 

miR- 155.  As  with  the  miR- 17-92  cluster, 
the  oncogenic  potential  of  miR- 155  has  been 
demonstrated  by  transgenic  mice  expressing 
miR- 155  in  B  cells  (127).  The  primary  tran¬ 
script  for  miR- 155  is  the  BIC  gene,  which  was 
originally  known  as  a  common  viral  integration 
site  in  lymphomas  in  chicken.  High  expression 
of  miR- 155  (BIC)  is  reported  in  various  B 
cell  malignancies  (128-130)  (Supplemental 
Table  1).  miR- 155  seems  to  be  regulated 
at  multiple  levels.  miR- 155  (BIC)  appears 
to  be  regulated  by  nuclear  factor-kappa  B 
(NF-kB)  through  the  B  cell  receptor-mediated 
signaling  pathway  (85)  or  through  the  Toll-like 
receptor-activated  signaling  pathway  (131). 
However,  a  discrepancy  between  the  levels 
of  miR- 155  and  BIC  suggests  that  miR- 155 
may  also  be  regulated  during  processing  steps 


(129).  For  example,  the  induction  of  BIC  does 
not  lead  to  upregulation  of  miR- 155  in  several 
Burkitt  lymphoma  cell  lines.  In  contrast, 
miR- 155  is  concordantly  increased  upon  the 
induction  of  BIC  in  Raji,  an  EBV-latency  type 
Ill-positive  Burkitt  lymphoma  cell  line  (85). 

let-  7.  let-  7  is  one  of  the  earliest-discovered 
miRNAs.  In  human,  there  are  12  paralogous 
let- 7s.  let- 7  is  almost  absent  during  embryonic 
stages  or  tissues,  whereas  high  expression  of 
let- 7  is  seen  in  most  differentiated  tissues.  The 
reduction  of  let-  7  in  cancers  is  reminiscent  of 
let-  7  expression  during  development  in  that 
it  is  most  decreased  in  less-differentiated,  ad¬ 
vanced  stages  of  cancer  cells  with  mesenchymal 
characters  (132). 

let-  7  is  probably  more  abundant  than  any 
other  miRNA,  so  reduction  of  let-  7  may  have 
a  prominent  effect  on  cell  physiology  com¬ 
pared  to  other  miRNAs.  To  support  this  hy¬ 
pothesis,  enhanced  cellular  transformation  by 
global  miRNA  depletion  was  largely  recapit¬ 
ulated  by  inhibition  of  let-  7  alone  (70).  Well- 
known  oncogenes  such  as  RAS ,  c-Myc ,  and 
HMGA2  are  validated  as  direct  targets  of  let- 
7  (Table  2).  Hence,  let- 7  is  an  important  tumor 
suppressor. 

miR-34s  (miR-34a  and  miR-34b/-34c  clus¬ 
ters).  The  miR- 3  4s  recently  acquired  notori¬ 
ety  because  they  are  induced  by  p5  3  (Figure  4). 
Consistent  with  p53’s  role  as  a  tumor  sup¬ 
pressor,  miR- 3  4s  are  downregulated  in  several 
tumors  such  as  non-small  cell  lung  cancers 
(60)  and  pancreatic  cancers  (61)  (Supplemental 
Table  1).  However,  reduction  of  the  miR- 3  4s  is 
not  always  correlated  with  p53  loss,  suggesting 
a  p5 3 -independent  mechanism  of  miR-34  re¬ 
duction  in  some  cancers.  In  fact,  miR-34a  is  lo¬ 
cated  in  lp36,  a  locus  that  is  frequently  deleted 
in  a  number  of  cancers  (76).  We  return  to  our 
discussion  of  p53  regulation  and  miR-34s’  roles 
below. 

Although  the  miR-34s  are  thought  to  be  tu¬ 
mor  suppressors,  they  have  been  found  to  be 
upregulated  in  several  cancers  including  renal 
cell  carcinoma  (133),  colon  cancer  (58),  and 
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Figure  4 

The  p53  and  miR-34  pathway.  Abbreviation:  MET,  hepatocyte  growth  factor  receptor. 


hepatocellular  carcinoma  (134)  (Supplemen¬ 
tal  Table  1).  The  role  of  high  miR-34  and  p53 
status  in  these  tumors  awaits  investigation. 

Aspects  of  Cancer  Biology  Regulated 
by  MicroRNAs 

Cell  cycle.  Cell-cycle  regulators  often  act  as 
oncogenes  or  as  tumor  suppressors.  The  best- 
characterized  example  is  the  cell-cycle  inhibitor 
p27(Kipl).  p27(Kipl)  is  a  tumor  suppressor,  as 
indicated  by  its  low  levels  in  some  cancers.  In 
addition,  p27(Kipl)  mutation  predisposes  cells 
to  tumorigenesis  upon  exposure  to  carcinogens. 
p27(Kipl)  binds  to  Cdk2-cyclin  E  and  prevents 
Gi-to-S  transition.  p27(Kipl)  is  a  direct  tar¬ 
get  of  miR-221  and  -222  in  glioblastomas  (48, 
135)  and  in  prostate  cancer  cells  (136).  In  these 
cancer  cells,  p27(Kipl)  is  anticorrelated  with 
miR-221  and  -222.  Targeting p27(Kip  1)  is  re¬ 
sponsible  for  the  proproliferative  role  of  these 
miRNAs,  as  artificial  knockdown  o{ p27 (Kip  1) 
mirrors  the  phenotype  of  these  miRNAs.  miR- 
221  and  -222  are  overexpressed  in  other  cancers 
(see  Supplemental  Table  1),  suggesting  that 
they  play  a  role  in  a  wide  range  of  cancers. 


In  addition  to  regulating  p27(Kipl), 
miRNAs  regulate  other  cell-cycle  proteins 
including  Cdk6;  Cdc25A;  Ccnd2  (cyclin  D2) 
(59);  Cdk4  (62),  a  Rb-family  protein  (89);  and 
pi 80  of  DNA  polymerase  a  (37).  A  number  of 
miRNAs  have  been  shown  to  perturb  normal 
cell  cycle  when  overexpressed  or  inhibited  (37, 
48,  62,  70,  137-139)  (Supplemental  Table  2). 
However,  oscillation  of  a  miRNA  during  the 
normal  cell  cycle  has  not  yet  been  reported. 

Programmed  cell  death.  Apoptosis  is  an  ac¬ 
tive  process  controlled  by  a  gene  expression 
program  that  varies  depending  on  the  biolog¬ 
ical  context.  Because  a  balance  between  pro¬ 
liferation  and  apoptosis  is  essential  for  tis¬ 
sue  homeostasis  and  proper  differentiation, 
aberrant  apoptosis  may  give  rise  to  tumors. 
miRNAs  participate  in  tumorigenesis  by 
directly  targeting  antiapoptotic  genes.  Repre¬ 
sentative  examples  include  the  repression  of  an¬ 
tiapoptotic  genes  Mcl-1  and  Bcl-2  by  miR-29b 
(84)  and  by  miR-34s  (60),  -15a,  and  -16 
(140),  respectively  (Table  2).  The  loss  of  these 
miRNAs  due  to  mutation  of  p53  or  deletion 
of  chromosome  13ql4  leads  to  an  increase  in 
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the  antiapoptotic  gene  expression  and  persis¬ 
tence  of  tumor  cells  that  would  normally  have 
been  removed  by  apoptosis.  It  is  very  likely  that 
miRNAs  target  other  genes  in  the  apoptotic 
pathway,  as  transfection  or  expression  of  a  num¬ 
ber  of  miRNAs  is  associated  with  apoptosis  (69, 
76,  121)  (Supplemental  Table  2). 

p53.  p53,  a  sequence-specific  transcription 
factor,  is  known  as  the  guardian  of  the  genome 
owing  to  its  critical  role  in  regulation  of  the  cell 
cycle  and  apoptosis  in  the  face  of  genomic  dam¬ 
age.  Genotoxic  stress  and  oncogene  activation 
activate  p53  to  modulate  the  transcription  of 
several  target  genes.  p53  is  the  most  extensively 
studied  tumor  suppressor  and  its  importance  is 
underscored  by  mutation  of  p53  in  almost  50% 
of  human  cancers. 

miRNA  profiling  after  p53  induction  indi¬ 
cated  miR-34a,  -34b,  and  -34c  (collectively,  the 
miR-34s)  as  the  most  upregulated  miRNAs  (61, 
62, 137, 141)  (Figure  4).  These  miRNAs  are  in¬ 
duced  after  genotoxic  stress  in  a  p5 3 -dependent 
manner  in  vitro  and  in  vivo  (62,  141).  miR- 
34b  and  -34c  are  clustered  at  chromosome  11, 
whereas  miR-34a  is  located  in  a  separate  ge¬ 
nomic  locus.  Both  pri-miRNAs  are  directly  ac¬ 
tivated  by  p53 .  The  miR-34s  seem  to  be  critical 
downstream  effectors  of  p53,  as  ectopic  expres¬ 
sion  of  the  miR-34s  recapitulate  the  phenotype 
of  p53  activation.  The  miR-34s  promote  cell- 
cycle  arrest,  apoptosis,  and  senescence  (58,  60- 
62,  137,  141)  (Supplemental  Table  2).  These 
effects  are  explained  by  the  repression  of  sev¬ 
eral  direct  targets  of  the  miR-34s  such  as  Bel- 2 
(60),  Cdk4 ,  and  hepatocyte  growth  factor  receptor 
(62)  (Table  2). 

In  addition  to  the  miR-34s,  other  miRNAs 
may  be  important  in  the  p53  pathway.  miR- 
30c,  -103,  -26a,  -107,  and  -182  were  induced 
clearly,  albeit  less  robustly,  upon  DNA  damage 
in  a  p5 3 -dependent  manner  (61).  miR-26a  ex¬ 
pression  was  also  shown  to  be  dependent  on 
p53  (142).  In  another  approach,  searching  for 
p5 3 -binding  elements  in  DNA  sequence  near 
miRNAs  identified  miR-129  as  a  good  candi¬ 
date  for  regulation  by  p53  (143). 


Invasion  and  metastasis.  Features  of  malig¬ 
nant  tumors,  distinct  from  benign  tumors,  in¬ 
clude  invasion  and  metastasis.  Malignant  tu¬ 
mors  are  fatal,  mostly  due  to  their  capacity 
to  invade  neighboring  tissues  and  metastasize 
through  the  bloodstream  to  distant  organs.  An 
effect  of  miRNAs  on  invasion  and  migration  has 
been  reported  (Supplemental  Table  2).  Ec¬ 
topic  expression  of  miR-125  impairs  cell  motil¬ 
ity  and  invasion  in  a  breast  cancer  cell  line  (144), 
and  reduction  of  global  miRNA  expression  en¬ 
hances  migration  of  cells  (70). 

These  sporadic  in  vitro  observations  were 
followed  by  more  detailed  studies  on  miR- 
10b  and  miR-21.  miR-21  is  one  of  the  most 
frequently  upregulated  miRNAs  in  cancers 
(Supplemental  Table  1).  It  promotes  cell 
motility  and  invasion  by  directly  targeting 
PTEN  {phosphatase  and  tensin  homolog ),  a  tu¬ 
mor  suppressor  known  to  inhibit  cell  invasion 
by  blocking  the  expression  of  several  matrix 
metalloproteases  (134).  Another  pathway  was 
recently  reported  in  colorectal  cancers,  where 
miR-21  promotes  invasion,  intravasation,  and 
metastasis  by  downregulating  Pdcd4  (145). 

miR-lOb  is  the  other  miRNA  implicated 
in  metastasis.  In  metastatic  breast  cancer  cells, 
miR-lOb  is  upregulated,  presumably  as  a  result 
of  transcriptional  activation  by  a  transcription 
factor,  Twist.  Ectopic  expression  of  miR-lOb 
promotes  invasion,  intravasation,  and  metasta¬ 
sis  in  otherwise  noninvasive  or  nonmetastatic 
breast  cancer  cell  lines.  miR- 1 0b  directly  targets 
HOXD10 ,  whose  reduction  induces  the  expres¬ 
sion  of  a  well-characterized  prometastatic  gene, 
RhoC  (146). 

Angiogenesis.  Recruiting  blood  vasculature  is 
crucial  to  the  survival  of  neoplastic  cells.  So  far, 
the  miR- 17-92  cluster  has  been  characterized 
in  this  process.  The  stimulation  of  neovascu¬ 
larization  by  c-Myc  involves  a  downregulation 
of  antiangiogenic  factor  Tsp-1.  c-Myc  represses 
Tsp-1  and  a  related  protein,  connective  tissue 
growth  factor  ( CTGF )  by  activating  the  miR- 17- 
92  cluster.  Tsp-1  and  CTGF  appear  to  be  direct 
targets  of  miR- 19  and  -18,  respectively.  Ectopic 
expression  of  the  miR- 1 7 -92  cluster  is  sufficient 


www.annualreviews.org  •  MicroRN As  in  Cancer  215 


Annu.  Rev.  Pathol.  Mech.  Dis.  2009.4:199-227.  Downloaded  from  arjournals.annualreviews.org 
by  University  of  Virginia  Libraries  on  02/05/09.  For  personal  use  only. 


for  promoting  angiogenesis  (147).  A  recent  ob¬ 
servation  indicates  that  other  miRNAs — miR- 
378  and  -27a — may  play  a  role  in  angiogene¬ 
sis  (139,  148)  (Supplemental  Table  2).  Viral 
miRNAs  may  also  play  a  role  in  angiogenesis, 
as  Tsp-1  has  been  shown  to  be  a  direct  target  of 
KSHV  miRNAs  (112). 

Others.  The  microenvironment  inside  a  solid 
tumor  is  usually  hypoxic.  Tumor  cells  under  hy¬ 
poxia  tend  to  be  resistant  to  therapies  and  to 
have  a  poor  prognosis.  Upregulation  of  certain 
miRNAs  under  hypoxia  (149,  150)  (Table  3) 
suggests  that  these  miRNAs  may  influence  the 
phenotype  of  hypoxic  tumor  cells.  Cancers  are 
predisposed  by  other  external  factors,  such  as 
genotoxic  stress,  a  folate-deficient  diet,  and  ex¬ 
posure  to  arsenic.  There  are  a  few  reports  of 
miRNA  profiles  under  these  conditions  (151, 
152),  but  whether  miRNA  changes  are  respon¬ 
sible  for  the  cancer  predisposition  remains  to 
be  investigated. 

Clinical  Applications 

Diagnostic  tools.  miRNAs  may  be  used  as 
diagnostic  or  prognostic  tools,  as  miRNA  ex¬ 
pression  profiles  reflect  tumor  origin,  stage, 
and  other  pathological  variables.  miRNAs  can 
function  as  accurate  molecular  markers  be¬ 
cause  they  are  relatively  stable  and  resistant 
to  RNase  degradation — probably  because  of 
their  small  size  (35,153,154).  Researchers  have 
shown  that  miRNAs  can  be  isolated  and  quan¬ 
titated  from  formalin-fixed  paraffin-embedded 
(FFPE)  specimens.  qRT-PCR  and  microarray 
data  were  reliably  and  reproducibly  obtained 
from  FFPE  samples  that  had  been  routinely 
processed  and  stored  for  10  years.  The  data 
from  FFPE  samples  are  consistent  with  those 
from  frozen  samples  (155,  156). 

The  development  of  qRT-PCR  methods  has 
improved  the  sensitivity  of  miRNA  detection 
down  to  a  few  nanograms  of  total  RNA  (68, 
157,  158).  This  amount  can  easily  be  obtained 
by  fine-needle  aspiration  biopsies  (FNABs); 
in  fact,  there  has  been  a  report  of  successful 


miRNA  measurement  by  qRT-PCR  on  FNAB 
samples  (159). 

miRNA  markers  that  could  be  used  for  can¬ 
cer  diagnosis  are  becoming  available.  For  ex¬ 
ample,  miR-196a  is  high  in  pancreatic  ductal 
adenocarcinoma  (PDAC)  but  low  in  normal  tis¬ 
sues  and  chronic  pancreatitis.  miR-2 17  exhibits 
the  opposite  expression  pattern  (Supplemental 
Table  1).  Thus,  the  ratio  of  miR-196a  to  miR- 
217,  calculated  by  qRT-PCR  measurement  of 
the  two  miRNAs  from  a  tiny  amount  of  total 
RNA,  indicates  whether  the  sample  contains 
PDAC  (160).  Once  reliable  indicator  miRNAs 
are  chosen,  they  will  likely  yield  easy  and  accu¬ 
rate  tools  for  cancer  diagnosis. 

Cancer  therapeutic  tool.  A  number  of 
miRNAs  affect  the  growth  of  cancer  cells  in 
vitro  and  in  vivo  when  overexpressed  or  inhib¬ 
ited.  Therefore,  cancer  cell  growth  can  be  con¬ 
trolled  by  manipulating  miRNAs.  Overexpres¬ 
sion  or  inhibition  of  miRNAs  can  be  achieved  in 
several  ways.  Synthetic  miRNA  mimics  include 
siRNA-like  oligoribonucleotide  duplex  (161) 
and  chemically  modified  oligoribonucleotide 
(162).  Conversely,  miRNAs  can  be  inhibited  by 
variously  modified  antisense  oligonucleotides 
such  as  2'-0-methyl  antisense  oligonucleotide, 
antagomirs,  and  so  on  (reviewed  in  Reference 
163).  As  the  first  successful  tool  for  knock¬ 
down  of  a  miRNA  in  vivo,  antagomirs  are  of 
special  interest  (164).  Antagomirs  appear  to 
be  delivered  to  all  tissues  (except  brain)  af¬ 
ter  tail-vein  injections  into  mice.  The  ther¬ 
apeutic  value  of  an  antagomir  would  be 
greatly  enhanced  by  technical  improvements 
for  selective  tumor-specific  or  tissue-specific 
delivery. 

Synthetic  oligonucleotides  are  effective  in 
vivo  for  at  most  a  couple  of  weeks,  as  has 
been  demonstrated  by  experiments  involving 
cancer  cells  engrafted  in  mice  (165)  and  tail- 
vein  injection  to  mice  (166).  To  circumvent 
this  limitation,  miRNAs  can  be  stably  expressed 
through  transcription  of  hairpin  RNA  from 
plasmid  vector  (reviewed  in  Reference  167). 
Recently,  artificial  overexpression  of  a  miRNA 
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target  sequence  was  shown  to  inhibit  the 
miRNA  function,  presumably  by  titrating  the 
miRNA  away  from  endogenous  targets  (168, 
169).  Thus  it  should  be  possible  to  apply  such 
competitive  inhibitors  for  long-term  sequestra¬ 
tion  of  a  miRNA. 

Nonspecific  side  effects  are  as  important  as 
effectiveness  and  duration  of  miRNA  expres¬ 
sion  (or  inhibition).  Some  investigators  have 
argued  that  miRNA  mimics  or  inhibitors  are 
specific  enough  to  distinguish  between  similar 
miRNAs  (170,  171).  However,  cross-reactivity 
between  miRNAs  of  similar  sequence  is  likely 
to  be  unavoidable  at  high  doses  of  antagonists 
or  agonists.  Another  possible  side  effect  is  that 
high  expression  of  miRNA  mimics  may  inter¬ 
fere  with  endogenous  miRNA  action  by  satu¬ 
rating  the  cellular  machinery  for  miRNA  pro¬ 
cessing  or  action.  This  may  result  in  a  change  in 
expression  of  other  miRNAs,  leading  to  a  dele¬ 
terious  effect  in  the  cells.  Indeed,  a  fatal  side 
effect  as  a  result  of  saturation  of  miRNA  path¬ 
way  has  been  reported  (172).  To  minimize  un¬ 
desirable  side  effects,  the  expression  or  knock¬ 
down  of  a  miRNA  should  be  improved  so  that 
it  is  more  accurate  and  controllable.  An  alter¬ 
nate  approach  to  improving  specificity  is  to  tar¬ 
get  the  pre-miRNAs  with  antisense  or  siRNA 
strategies  (173,  174). 

Some  chemical  compounds  alter  expression 
of  a  group  of  miRNAs  (81);  therefore,  it  may 
be  possible  to  screen  for  drugs  that  could  shift 
the  miRNome  in  a  cancer  cell  toward  that  of 
the  normal  tissue.  NCI-60  cancer  cell  lines 
may  serve  as  a  platform  for  this  screening, 
as  their  miRNA  expression  profiles  have  been 
published  (68,  175).  By  modulating  multiple 
miRNAs  simultaneously,  such  a  miRNome- 
modifying  approach  may  be  much  more  effec¬ 
tive  for  therapy  than  strategies  that  aim  to  reg¬ 
ulate  a  single  miRNA. 

miRNAs  affect  the  drug  sensitivity  of  a  cell 
(165,  176)  (Supplemental  Table  2).  Expres¬ 
sion  or  inhibition  of  a  miRNA  can  therefore 
be  combined  with  treatment  of  a  drug  or  other 
cytotoxic  therapy.  One  example  is  miR-2 1  inhi¬ 


bition  together  with  a  secreted  form  of  tumor 
necrosis  factor-related  apoptosis-inducing  lig¬ 
and,  which  results  in  a  complete  eradication  of 
glioblastoma  cells  (177). 

Collectively,  preliminary  results  suggest  that 
miRNAs  could  be  useful  for  cancer  therapy. 
However,  there  is  still  a  significant  gap  between 
basic  research  on  miRNAs  and  clinical  applica¬ 
tion.  Extensive  preclinical  and  translational  re¬ 
search  is  necessary  to  increase  the  efficacy  and 
decrease  the  side  effects  of  miRNAs  in  vivo.  In 
parallel,  we  need  to  expand  our  knowledge  of 
the  interactions  between  miRNAs  and  gene  ex¬ 
pression  programs  and  how  these  interactions 
are  altered  in  tumorigenesis. 

CONCLUDING  REMARKS 

The  discovery  of  miRNAs  has  added  another 
dimension  to  the  study  of  the  regulation  of 
gene  expression.  In  addition  to  transcriptional 
regulation,  posttranscriptional  repression  by 
miRNAs  may  act  like  rheostats  for  fine-tuning 
of  gene  expression.  miRNAs  at  end  stages  of 
differentiation  may  contribute  to  keeping  cells 
differentiated  by  suppressing  a  large  number  of 
genes  simultaneously. 

In  recent  years,  there  has  been  an  explosion 
of  publications  on  miRNAs.  Many  of  these  have 
emphasized  the  role  of  miRNAs  in  cancer  biol¬ 
ogy.  miRNAs  play  a  role  in  almost  all  aspects  of 
cancer  biology,  including  proliferation,  apopto¬ 
sis,  invasion/metastasis,  and  angiogenesis.  New 
technical  developments  are  enabling  investiga¬ 
tors  to  describe  a  variety  of  cancer  miRNomes. 
Therefore,  more  miRNAs  are  expected  to  be 
identified  as  oncogenes  or  as  tumor  suppressors. 
Characterization  of  individual  miRNA  path¬ 
ways  will  be  greatly  facilitated  by  an  improve¬ 
ment  of  target  prediction  algorithms. 

miRNA  expression  is  associated  with  clinical 
variables  of  cancers,  so  miRNAs  can  be  readily 
used  as  tools  for  cancer  diagnosis  and  progno¬ 
sis.  For  therapeutic  benefit,  the  methods  to  ma¬ 
nipulate  miRNAs  in  vivo  should  become  more 
robust. 
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SUMMARY  POINTS 


1 .  miRNAs  exhibit  an  expression  pattern  characteristic  of  tumor  type,  stage,  or  other  clinical 
variables.  Thus  miRNAs  can  be  used  for  cancer  diagnosis  and  prognosis. 

2.  miRNAs  exert  an  active  role  in  the  etiology  or  progression  of  cancers  by  regulating  their 
direct-target  oncogenes  or  tumor-suppressor  genes. 

3 .  miRNAs  play  roles  in  diverse  aspects  of  cancer  such  as  proliferation,  apoptosis,  inva¬ 
sion/metastasis,  and  angiogenesis. 

4.  Artificial  expression/inhibition  of  a  miRNA  can  be  therapeutic. 

APPENDIX 

Other  Small  RNAs? 

In  addition  to  siRNAs  and  miRNAs,  other  types  of  small  RNA  have  been  reported  in  yeast, 
Drosophila ,  C.  elegans ,  and  plants.  These  small  RNAs  include  tiny  noncoding  RNAs,  repeat  associ¬ 
ated  small  interfering  RNAs,  scan  RNAs,  and  others  (reviewed  in  Reference  178).  These  RNAs  are 
thought  to  be  related  to  miRNA  and  siRNA,  but  are  distinct  from  them.  Recently,  Piwi-interacting 
RNA,  another  small  RNA  that  is  slightly  larger  than  miRNA,  was  discovered  from  mammalian 
testes.  Other  classes  of  small  RNA  may  exist  but  have  yet  to  be  discovered,  given  that  a  signif¬ 
icant  portion  of  transcripts  are  noncoding  RNA.  Whereas  protein-coding  genes  comprise  only 
1-2%  of  the  human  genome,  a  significant  fraction  (~15%)  of  the  human  genome  is  transcribed; 
these  transcripts  remain  the  “dark  matter”  of  the  genome  (179,  180).  The  discovery  of  noncod¬ 
ing  RNAs  is  expected  to  be  greatly  accelerated  by  the  widespread  adaptation  of  high-throughput 
deep-sequencing  techniques. 


Other  Diseases 

miRNAs  (or  miRNA  machinery)  have  been  implicated  in  diseases  other  than  cancer  (182),  for 
example  Tourette’s  syndrome,  fragile  X  syndrome  (181),  DiGeorge  syndrome  (23),  myotonic 
dystrophy,  spinocerebellar  ataxia  type  3  (182),  and  schizophrenia  and  schizoaffective  disorder 
(183).  Of  particular  interest  is  myotonic  dystrophy  type  1,  which  is  characterized  by  expansion  of 
CTG  repeat  in  the  3'  UTR  of  the  dmpk  gene.  Tandem  CAG  sequences  are  found  in  some  miRNA 
seed  sequences  and  thus  are  able  to  interact  with  the  CTG  repeats,  raising  the  possibility  that 
these  miRNAs  may  be  involved  in  DM1  pathogenesis  (184). 
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Abstract 


A  systematic  comparison  of  the  short  RNA  universe  of  two  prostate  cancer  cell  lines  was 
performed  by  cloning  and  deep  sequencing  to  identify  microRNAs  whose  expression  may 
correlate  with  androgen  dependence/independence  or  microRNAs  that  are  regulated  by 
androgens.  The  results  were  confirmed  by  two  independent  methods  of  quantitating  microRNAs: 
locked  nucleic  acid  microarrays,  and  microRNA-specific  qPCR/RNase  protection/Northem 
assays.  Changes  that  were  consistently  detected  by  all  three  methods  are  the  decrease  of  five 
microRNAs  and  increase  of  three  microRNAs  in  the  androgen-independent  C4-2  cell  line  relative 
to  its  androgen-dependent  predecessor,  LNCaP.  Ectopic  expression  of  the  down-regulated 
microRNAs  decreased  the  growth  of  C4-2  in  the  absence  of  androgen  but  not  in  its  presence, 
suggesting  that  microRNA  changes  contribute  to  the  androgen-independence.  All  five 
downregulated  microRNAs  have  been  reported  to  be  downregulated  in  clinically  advanced 
prostate  cancer,  suggesting  that  they  constitute  a  microRNA  signature  of  prostate  cancer 
progression.  Discrepancy  in  the  profile  obtained  by  cloning  and  sequencing  versus  microarray 
hybridization  suggest  that  nine  microRNAs  suffer  an  as  yet  unknown  3’  end  modification  in 
androgen-depleted  cells  that  decreases  3’  adaptor  ligation.  Finally,  cloning  and  sequencing 
reveals  a  large  number  of  non-micro-small  RNAs  of  unknown  function  that  are  expressed 
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abundantly  in  the  prostate  epithelial  cells  and  may  be  as  useful  for  profiling  cancers  as 
microRNAs  or  mRNAs. 
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Introduction 


Small  RNAs  provide  a  new  paradigm  in  understanding  regulation  of  gene  expression. 
microRNA  (miRNA)  and  small  interfering  RNA  (siRNA)  are  two  representative  members  of 
small  non-coding  regulatory  RNAs.  Other  types  of  small  RNA  molecules  have  also  been  reported, 
including  tiny  noncoding  RNAs  (tncRNAs),  repeat  associated  small  interfering  RNAs 
(rasiRNAs),  scan  RNAs  (scnRNAs)  and  Piwi-interacting  RNA  (piRNA)  (reviewed  in  (Tolia  and 
Joshua-Tor  2007)).  These  RNAs  are  thought  to  have  functions  distinct  from  miRNAs  and  siRNAs. 

Among  these  small  RNAs,  miRNAs  have  been  studied  most  extensively.  Each  miRNA 
exhibits  a  characteristic  expression  profile,  with  levels  of  expression  of  a  few  specific  miRNA 
being  dramatically  altered  by  tissue  type  and  developmental  stage.  miRNAs  bind  to  the 
3 ’untranslated  region  (3’UTR)  of  a  target  mRNA  and  down-regulate  its  expression.  Many  studies 
demonstrated  that  miRNAs  play  critical  roles  in  cancer  (reviewed  in  (Lee  and  Dutta  2006)). 

Prostate  cancer  is  the  most  frequently  diagnosed  cancer  and  the  second  leading  cause  of 
cancer-related  deaths  in  the  male  population  of  the  United  States.  Initially,  prostate  cancers 
depend  upon  androgens  for  their  growth.  Therefore,  the  primary  treatment  for  metastatic  prostate 
cancer  is  androgen  deprivation  therapy  (ADT),  achieved  by  orchiectomy  or  anti-androgens. 
However,  prostate  cancer  often  progresses  into  an  androgen-independent  stage  and  becomes 
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resistant  to  ADT  (reviewed  in  (Suzuki  et  al.  2003)  and  (Agoulnik  and  Weigel  2006)).  Thus,  it  is 


of  clinical  importance  to  understand  the  mechanism  of  the  transition  from  androgen-dependence 
to  androgen-independence. 

Several  models  have  been  proposed  to  explain  the  molecular  basis  of  how  prostate  cancer 
cells  bypass  the  requirement  of  androgen  for  their  growth  (reviewed  in  (Pienta  and  Bradley 
2006)).  The  androgen  receptor  (AR),  a  member  of  the  superfamily  of  nuclear  hormone  receptors, 
is  a  key  molecule  in  the  androgen  pathway.  Upon  binding  to  androgen,  AR  undergoes  a 
conformational  change,  translocates  to  the  nucleus,  and  modulates  the  transcription  of  many 
target  genes  by  binding  to  androgen  response  elements  (ARE)  in  their  promoter  regions.  The 
coordinated  gene  expression  by  androgen  through  the  AR  pathway  is  vital  for  the  survival  and 
proliferation  of  prostate  cells.  Therefore,  an  activating  mutation  in  AR,  overexpression  of  AR  or 
aberrant  expression  of  AR  coregulators  might  make  the  AR  active  under  low  androgen  conditions, 
thereby  allowing  prostate  cancer  cells  to  proliferate  even  after  androgen  depletion  (reviewed  in 
(Chmelar  et  al.  2007)).  In  another  model,  prostate  cancer  cells  might  activate  a  bypass  signal 
transduction  pathway  so  as  to  survive  independent  of  androgen  and  AR. 

As  described  above,  reprogramming  of  gene  expression  is  the  basis  of  many  of  the 
proposed  mechanisms  by  which  prostate  cancer  cells  become  androgen-independent  and  by 
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which  androgens  regulate  prostate  epithelial  proliferation.  miRNAs  are  potent  regulators  that  can 


inhibit  many  genes  simultaneously  due  to  their  small  size  and  relaxed  specificity  in  target 
recognition.  Introduction  of  a  single  miRNA  into  cells  can  lead  to  a  reprogramming  of  global 
gene  expression  (Lim  et  al.  2005).  Since  miRNAs  have  been  implicated  in  a  variety  of  tumors, 
we  reasoned  that  miRNAs  might  play  a  role  during  prostate  cancer  progression  into  androgen- 
refractory  stage  and  in  the  androgen  regulation  of  androgen-dependent  cells.  Although  there  are 
two  widely  accepted  methods  for  global  profiling  miRNAs  (miRNomes),  a  head-to-head 
comparison  of  the  methods  has  not  been  published.  We,  therefore,  decided  to  profile  the  miRNAs 
in  prostate  cancer  cells  by  cloning  and  deep  sequencing  (Margulies  et  al.  2005)  and  follow  up  by 
hybridization  to  microarrays  made  of  locked  nucleic  acid  probes  (LNA-microarray)  (Castoldi  et 
al.  2006)  and  to  confirm  a  subset  of  the  changes  by  individual  miRNA-specific  biochemical 
assays.  The  results  identify  miRNAs  whose  downregulation  contributes  to  androgen- 
independence  and  serves  as  a  signature  of  advanced  prostate  cancer,  suggest  that  constitutive 
activity  of  AR  is  not  the  explanation  of  androgen-independence  of  C4-2,  identify  a  potential 
modification  at  the  3’  ends  of  specific  miRNAs  in  androgen-depleted  cells  and  uncover  a  large 
assortment  of  novel  and  abundant  small  RNAs  that  are  not  miRNAs  (nmsRNAs). 
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Results 


Small  RNA-ome  in  prostate  cancer  cells.  Small  RNAs  of  17-26  nts  were  purified  from  two 
prostate  cancer  cell  lines,  LNCaP  and  C4-2.  The  growth  of  LNCaP  is  dependent  on  androgen 
while  C4-2  is  a  cell  line  derived  from  LNCaP  that  is  independent  of  androgen  (Wu  et  al.  1994) 
and  is  more  invasive  and  metastatic  (Thalmann  et  al.  1994;  Wu  et  al.  1994;  Thalmann  et  al.  2000). 
In  addition,  to  identify  androgen-regulated  miRNAs,  short  RNAs  were  purified  from  LNCaP  cells 
grown  in  charcoal  stripped  serum  or  the  same  supplemented  with  DHT  (4,5-dihydrotestosterone), 
an  active  ligand  of  AR  (Anderson  and  Liao  1968;  Bruchovsky  and  Wilson  1968).  Charcoal 
stripped  serum  (termed  “-DHT”,  hereafter)  effectively  depleted  the  cells  of  androgen,  as 
indicated  by  RT-PCR  detection  of  PSA  (prostate  specific  antigen)  mRNA,  a  representative  target 
gene  of  the  AR  (Fig  SI).  The  level  of  PSA  is  restored  to  that  of  whole  serum  conditions  when 
DHT  is  added  to  -DHT  condition,  showing  AR  activity  is  rescued  by  DHT  (termed  “+DHT”). 

From  the  four  samples  (C4-2,  LNCaP,  -DHT  and  +DHT),  small  RNAs  of  17-26 
nucleotides  were  isolated  and  ligated  sequentially  to  the  3’  and  the  5’  adaptor  (see  Fig  3  for  a 
brief  workflow).  cDNA  from  the  ligated  RNA  was  concatamerized  before  454  deep  sequencing, 
in  order  to  increase  the  yield  of  the  amount  of  sequence  reads  per  run.  A  flowchart  for  the 
analysis  of  454  deep  sequence  data  (“454  data”  hereafter)  is  outlined  in  Fig  S2.  More  than  half 
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million  (658,068)  sequences  were  obtained  from  the  four  samples  (Table  1).  The  inserts  are 
compared  with  one  another,  to  yield  unique  sequences  and  their  cloning  frequencies.  More  than  a 
thousand  unique  sequences  from  each  sample  were  cloned  five  times  or  more,  with  a  few 
hundred  being  cloned  more  than  50  times.  The  number  of  unique  sequences  relative  to  total 
sequences  in  -DHT  (13.81%)  is  significantly  higher  than  the  other  samples,  suggesting  the 
existence  of  more  diverse  small  RNAs  in  cells  growing  in  charcoal  stripped  serum.  Global 
degradation  of  the  RNA  from  the  androgen-depleted  cells  was  ruled  out  by  ethidium  bromide 
staining  of  the  RNA  on  a  denaturing  agarose  or  acrylamide  gel  (Fig  S3)  and  so  cannot  be  the 
explanation  for  this  increase  in  diversity. 

miRNAs  were  extracted  from  the  unique  sequences  by  searching  against  miRNA 
database  (miRbase  release  10.0;  http://microrna.sanger.ac.uk/sequences/).  37%  of  total  small 
RNAs  match  miRNAs  (Table  1).  Out  of  533  human  miRNAs  in  the  database,  roughly  half  (293 
miRNAs  as  a  whole  in  the  union  of  the  four  sets,  with  221  to  240  in  each  individual  sample)  were 
cloned  from  the  prostate  cancer  cells.  The  cloning  frequency  of  each  miRNA  was  postulated  to  be 
proportional  to  the  quantity  of  the  miRNA.  Table  SI  contains  the  raw  data  for  all  miRNAs.  The 
most  frequently  cloned  miRNAs  were  let-7  family  members,  miR-125b,  -99a,  -200c,  -17,  and  -21, 
suggesting  that  these  miRNAs  are  abundant  in  these  cell  lines  (Table  SI). 


Lee,  Yong  Sun 


To  corroborate  these  profiles,  the  RNAs  from  all  four  samples  were  independently 


hybridized  to  LNA-microarray  of  miRNA probes.  The  raw  data  is  tabulated  in  Table  SI.  Based  on 
the  454  and  microarray  data,  miRNAs  with  interesting  changes  in  expression  pattern  were 
collected  and  summarized  in  Table  S2  and  S3.  These  miRNAs  were  independently  measured  by 
individual  miRNA-specific  assays  (Fig  1 ;  see  also  Fig  S4-6)  to  confirm  changes  in  their 
expression  between  samples. 

miRNAs  differentially  expressed  in  LNCaP  and  C4-2  There  were  reproducible  changes  in 
several  miRNAs  between  these  two  cell  lines  in  multiple  assays  (Table  2  and  S2;  Fig  1A  and  left 
two  lanes  of  Fig  IB;  see  also  lanes  9-10  of  Fig  S4A-B).  In  general  we  give  greater  weight  to  the 
changes  that  are  confirmed  and  accurately  quantitated  by  the  miRNA-specific  biochemical  assays 
(qPCR,  RNase  protection  or  Northern).  miR-100,  -125b,  -19b,  and  -99a  were  most  down- 
regulated  in  C4-2  relative  to  LNCaP.  The  miRNAs  in  the  miR-17-92  cluster  (miR-17-5p,  -18a,  - 
20a  and  -92a)  were  also  decreased  in  C4-2  relative  to  LNCaP.  miR-106a,  -99b,  -21,-16  were 
decreased,  but  not  to  as  great  an  extent  as  the  others. 

In  contrast,  miR-9  and  -557  are  significantly  up-regulated  in  C4-2  relative  to  LNCaP  in 
microarray  and  miRNA-specific  qPCR.  Also,  a  moderate  increase  of  of  miR-196b  in  microarray 
was  confirmed  by  miRNA-specific  qPCR  assays.  The  up-regulation  of  miR-196b  is  also 


Lee,  Yong  Sun 


9 


supported  by  the  5 -fold  higher  cloning  frequency  in  C4-2  than  LNCaP. 


Altering  the  miRNA  profile  of  C4-2  changes  its  growth  property  in  charcoal-stripped  serum. 
Having  observed  a  significant  decrease  of  miR-100,  -125b,  -19b,  and  -99a  in  C4-2,  we  tested 
whether  ectopic  expression  of  these  miRNAs  confers  androgen-dependence  on  C4-2.  miR-99b 
was  also  included  in  this  experiment,  because  it  is  similar  to  miR-99a  in  sequence  and  is 
decreased  in  C4-2  (Fig  1A).  Upon  transfection  of  a  mixture  of  these  five  miRNAs,  the  growth  of 
C4-2  was  not  affected  in  untreated  serum  but  was  significantly  reduced  in  charcoal  stripped 
serum  (Fig  2A).  Transfection  of  an  individual  miRNA,  at  the  same  concentration  that  it  was  in 
the  mixture,  showed  that  the  selective  growth  inhibition  in  charcoal-stripped  serum  was  recreated 
by  transfection  of  any  of  four  miRNAs  in  the  mix,  except  for  miR-125b  (Fig  2B).  Thus,  the 
reduction  of  miR-100,  -19b,  -99a,  and  -99b,  seen  during  the  progression  from  LNCaP  to  C4-2, 
could  contribute  to  the  transition  from  androgen-dependence  to  -independence. 
miRNAs  differentially  expressed  upon  addition  of  androgen.  The  addition  of  DHT  to  LNCaP 
cells  growing  in  charcoal-stripped-serum  produced  only  a  few  miRNA  changes  that  were 
reproduced  in  two  out  of  the  three  assays  (Table  2  and  S3;  Fig  1C  and  right  two  lanes  of  Fig  IB; 
see  also  lanes  7-8  in  Fig  S4A-B  and  Fig  S6).  Again,  greater  weight  was  given  if  a  change  was 
reproduced  and  accurately  quantitated  by  the  miRNA-specific  assays.  These  changes  were  a  2.5- 
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fold  reduction  of  miR-221  and  smaller  reductions  in  miR-196b  and  -125b.  miR-125b  was  also 
measured  by  Northern  hybridization  (lanes  7-8  in  Fig  S4A-B)  and  shown  to  be  reduced  similarly 
to  qPCR  assays.  Therefore,  the  reduction  of  miR-125b  in  the  presence  of  DHT  is  moderate  but 
reliably  reproduced  by  multiple  assays.  miR-99a  is  reduced  in  RNase  protection  assays  (lane  3-4 
in  Fig  IB),  although  this  reduction  is  not  marked  in  qPCR  assays  (Fig  1C).  The  repression  of 
miR-221,  -196b,  and  -125b  was  reproduced  by  R1881,  an  androgen  analog  (Fig.  ID;  lanes  1-6  in 
Fig  S4A-B;  see  Fig  S4C  for  PSA  level  upon  R1881  treatment).  The  slight  suppression  of  miR- 
99a  by  DHT  (Fig  1C)  was  amplified  by  R1881  (Fig  ID).  Upon  finding  that  two  of  the  androgen- 
depleted  miRNAs  (miR-125b  and  -99a)  were  also  depleted  in  the  androgen-independent  cell  line 
C4-2  (Fig  1  and  Fig  S5),  we  wondered  whether  the  miRNome  changes  in  C4-2  could  be 
explained  simply  by  the  constitutive  activation  of  AR.  However,  the  decrease  of  miR-99b  and  - 
19b  or  the  increase  of  miR-557  and  -9  in  C4-2  was  not  recapitulated  upon  addition  of  androgen  to 
LNCaP  cells  (Fig  S5).  Conversely,  the  repression  of  miR-221  and  -196b  by  androgen  in  LNCaP 
cells  was  not  paralleled  by  a  decrease  of  these  miRNAs  in  C4-2.  Thus  we  believe  that  the  miRNA 
changes  in  C4-2  are  not  solely  the  result  of  constitutive  activation  of  AR. 

In  addition,  the  slight  induction  of  let-7  in  -DHT  compared  to  LNCaP  (growing  in 
serum)  was  more  or  less  abrogated  in  +DHT  (Fig  IB;  Table  2),  suggesting  that  these  miRNAs  are 
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modestly  repressed  by  androgen. 


In  summary,  miR-221,  -196b,  -125b  and  let-7  were  repressed  by  androgens  in  LNCaP 
cells.  We  were  unable  to  find  any  miRNAs  significantly  induced  by  androgen. 

Corroboration  between  methods.  Although  the  454  data,  the  LNA-microarray  data,  and  the 
results  from  other  detection  methods  are  largely  consistent  among  one  another  for  over  several 
hundred  miRNAs,  a  couple  of  intriguing  discrepancies  were  noticed  (Table  S2).  miR-100  and  - 
557  were  not  cloned  (Table  SI)  but  were  detected  by  the  array  hybridization  and  qPCR  assay.  In 
general,  all  the  hybridization-based  techniques  (microarray,  qPCR,  RNase  protection  assay, 
Northern  hybridization)  are  very  consistent  among  themselves.  One  exception  is  miR-130b, 
which  gave  comparable  signals  between  LNCaP  and  C4-2  on  the  arrays,  was  5-fold  higher  in 
LNCaP  in  the  qPCR  assays  but  was  less  frequently  cloned  in  LNCaP  than  in  C4-2. 

A  very  interesting  pattern  of  under-reporting  of  specific  miRNAs  by  the  cloning  method 
was  seen  in  androgen-depleted  cells  (Table  S3  and  Fig  S6).  Comparison  between  “-DHT”  and 
“+DHT”  suggests  that  some  miRNAs  are  oddly  underrepresented  in  miRNA  clones  obtained 
from  androgen-depleted  LNCaP,  even  though  they  are  easily  detected  by  the  hybridization  based 
methods  (bottom  nine  rows  in  Table  S3  and  Fig  S6).  This  observation  may  be  explained  by 
several  hypotheses,  one  of  which  is  that  modification  of  5’  or  3’  end  of  a  specific  miRNA 
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decreases  its  efficiency  of  cloning  in  androgen-depleted  cells  without  affecting  its  detection  by 


hybridization  based  methods. 

Using  miR-21  as  a  test  case,  we  first  measured  the  miRNA  by  other  assays  to  ascertain 
which  method  gave  the  accurate  estimate  of  miR-21  in  “-DHT”  cells.  miR-21  was  cloned  2053 
times  in  “+DHT”  but  only  490  times  in  “-DHT”  (Table  SI),  giving  the  relative  cloning  ratio  5.86 
(Table  SI  and  S3;  Fig  S6)  after  normalization  to  each  corresponding  total  number  of  sequence.  In 
contrast,  comparable  levels  of  miR-21  were  measured  in  “-DHT”  and  “+DHT”  by  Northern 
hybridization  (Fig  IB  and  S7C)  and  qPCR  of  total  RNA  and  small  RNA  tailed  by  poly-A  (Fig  3;  i 
and  ii).  In  a  functional  assay,  where  we  measured  the  repression  of  a  luciferase  reporter 
containing  a  perfect  match  to  miR-21  in  its  3’UTR,  the  repressive  activity  of  endogenous  miR-21 
was  also  comparable  between  “-DHT”  and  “+DHT”  (Fig  S7A). 

We  next  checked  the  yield  of  intermediate  products  after  each  of  the  cloning  steps. 
Although  the  amounts  of  starting  miR-21  were  equal  in  the  two  samples,  3’  adaptor  ligation  to 
miR-21  from  “-DHT”  cells  is  less  efficient  than  that  from  “+DHT”.  This  is  revealed  by  qPCR 
amplification  of  cDNA  from  the  3’  ligated  RNA  using  the  3’  adaptor  primer  and  miR-21 -specific 
primer  (Fig  3;  iii).  This  discrepancy  persisted  (and  is  slightly  amplified)  when  the  ligation 
products  were  input  into  the  5’  adaptor  ligation  reaction  and  concatamerized  (Fig  3;  iv  and  v). 
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This  extra  decrease  is  not  due  to  an  inefficiency  of  5’  adaptor  ligation,  as  demonstrated  by  direct 


5’  adaptor  ligation  to  the  small  RNA  (Fig  3;  vii).  Our  results  suggest  that  the  lower  efficiency  of 
3’  adaptor  ligation  to  miR-21  from  “-DHT”  cells  was  the  primary  reason  why  miR-21  clones 
were  under-represented  in  this  sample. 

Novel  small  RNAs  other  than  miRNAs.  Approximately  63%  of  the  total  sequence  yield  and 
90%  of  unique  clones  (Tables  1  and  S4)  of  the  small  RNAs  do  not  match  with  miRNA  sequence 
(non-micro- short  RNA,  nmsRNA).  A  significant  portion  of  these  were  cloned  multiple  times  and 
so  could  not  be  attributed  to  mistakes  in  sequencing.  Tens  of  these  sequences  were  cloned  >100 
times  in  at  least  one  of  the  four  samples,  indicating  that  these  novel  small  RNAs  are  present  in 
abundance.  In  Table  3  are  summarized  a  representative  list  of  nmsRNAs  that  mapped  to  genomic 
sites  and  were  cloned  at  a  high  frequency.  20-22  %  of  the  nmsRNAs  are  mapped  within  known 
RNAs  such  as  rRNA,  tRNA,  snoRNA,  piRNA,  or  Refseq  genes,  suggesting  that  the  nmsRNAs 
are  derived  from  these  RNAs  (Table  S4).  Northern  hybridization  confirmed  the  existence  of  two 
nmsRNAs  derived  from  snoRNA  U3  and  a  nmsRNA  from  a  tRNA  (Table  3  and  Fig  S8).  20%  of 
the  nmsRNAs  map  to  the  human  genome  sequence  if  we  allow  two  mismatches  at  the  end  due  to 
sequencing  errors,  suggesting  that  they  are  derived  from  as  yet  unidentified  genes. 

Some  nmsRNAs  were  perfectly  aligned  to  a  part  of  human  piRNAs  or  miRNAs  from 
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chicken  (see  “note”  column  in  Table  3).  However,  it  is  very  unlikely  that  they  are  novel  miRNAs 


in  human,  as  no  secondary  structure  characteristic  of  hairpin-shaped  miRNA  precursors  was 
observed  from  the  sequence  as  embedded  in  its  flanking  genomic  or  EST  sequence  (data  not 
shown).  Some  nmsRNAs  were  encoded  within  a  repeat  element  (Table  3),  suggesting  that 
rasiRNA-like  RNA  may  exist  in  the  prostate  epithelial  cells.  In  a  few  cases,  a  nmsRNA  did  not 
match  against  the  human  genome,  but  matched  to  human  ESTs  (for  example,  ID000653_3313  in 
Table  3),  suggesting  the  possibility  of  base  modification  at  the  RNA  level. 

The  nmsRNAs  that  do  not  map  back  to  the  genome  may  be  explained  by  errors  during 
library  construction  (i.e.  during  PCR  amplification)  or  454  deep  sequencing  but  is  most  likely  due 
to  sequence  variations  introduced  by  post-transcriptional  modifications.  Recent  genome-wide 
transcript  mapping  studies  have  suggested  that  nearly  67%  of  the  genome  is  transcribed  at  low 
levels  with  ample  evidence  of  splicing  and  trans-splicing  of  transcripts  ((Kapranov  et  al.  2007) 
and  T.  Gingeras,  personal  communication).  Thus,  21  base  junctional  fragments  from  such  splices 
will  be  difficult  to  map  back  to  the  genome  by  BLAST.  RNA  editing  may  also  alter  the  small 
RNA  sequence  relative  to  the  genomic  sequence. 
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Discussion 


Recent  publications  have  reported  miRNomes  from  prostate  cancer  cell  lines  and  clinical 
samples  obtained  by  microarrays  (Mattie  et  al.  2006;  Volinia  et  al.  2006;  Porkka  et  al.  2007;  Lin 
et  al.  2008;  Ozen  et  al.  2008).  Observing  that  there  were  many  discrepancies  between  these 
reports  as  to  which  miRNAs  were  altered  and  by  how  much,  we  decided  to  first  use  cloning  and 
deep  sequencing  to  profile  the  universe  of  short  RNAs  (microRNAs  and  non-micro- short  RNAs) 
of  two  well  defined  cell  lines  under  controlled  conditions.  The  results  were  then  validated  by 
multiple  approaches  to  ascertain  the  reliability  of  the  miRNA  changes  noted.  We  thus  arrived  at  a 
validated  list  of  miRNAs  altered  in  expression  level  between  the  androgen-dependent  cell  line, 
LNCaP,  and  its  androgen-independent  derivative,  C4-2,  and  upon  addition  of  androgen  to  the 
androgen-dependent  LNCaP  cells. 

One  reason  for  the  discrepancies  between  various  profiling  approaches  is  the  small 
magnitude  of  changes  reported.  Consistent  with  this,  most  of  the  miRNA  changes  reported  in  this 
study  are  less  than  three-fold  and  there  are  significant  differences  in  the  ratios  determined  by  the 
different  methods.  However,  because  of  our  controls,  we  could  identify  five  miRNAs  (miR-99a, 
-99b,  -100,  -125b  and  -19b)  that  are  decreased  and  three  miRNAs  (miR-9,  -196b  and  -557)  that 
are  increased  in  C4-2  relative  to  LNCaP.  C4-2,  cells  that  are  normally  not  affected  by  serum 
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deprivation,  suffered  a  decrease  in  proliferation  in  charcoal-stripped  serum  upon  transfection  of 


the  miRNAs  that  are  decreased  in  these  cells.  That  proliferation  was  not  decreased  in  complete 
medium  suggests  that  this  approach  will  be  useful  to  find  miRNAs  that  contribute  to  the 
androgen-independence  of  C4-2. 

These  miRNAs  will  be  tested  in  the  future  as  molecular  markers  for  prostate  cancer 
progression.  The  utility  of  miRNA  in  cancer  diagnosis  was  discussed  in  several  reports.  An 
expression  profile  of  miRNAs  was  shown  to  be  better  at  predicting  and  classifying  cancer  types 
and  stages  than  the  mRNA  expression  profile  (Lu  et  al.  2005).  Among  the  down-regulated 
miRNAs  in  C4-2,  ectopic  expression  of  miR-19b,  -100,  -99a  and  -99b  resulted  in  a  growth 
inhibitory  effect  on  C4-2  selectively  in  the  absence  of  androgen.  This  opens  a  possibility  that 
delivery  of  these  miRNAs  combined  with  anti-androgen  treatment  could  be  considered  for  a 
therapy  for  an  advanced,  androgen-independent  prostate  cancer. 

In  agreement  with  our  observation  that  miR-100,  -125b,  -19b,  -99a  and  -99b  are  down- 
regulated  in  C4-2  than  LNCaP,  all  of  them  were  shown  to  be  markedly  decreased  in  more 
aggressive  (Gleason  score  8)  and  metastatic  prostate  tumor  (Lu  et  al.  2005;  Mattie  et  al.  2006). 
Given  that  C4-2  is  readily  metastatic  and  more  invasive  than  LNCaP  (Thalmann  et  al.  1994;  Wu 
et  al.  1994;  Thalmann  et  al.  2000),  these  miRNAs  may  also  play  a  role  in  metastasis  and  invasion. 
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Our  results  are  also  supported  by  another  profiling  study  where  miR-99a,  -125b,  -19b  and  -100 


were  shown  to  be  reduced  in  prostate  carcinomas,  with  miR-19b  and  -100  being  reduced 
especially  in  hormone-refractory  carcinomas  (Porkka  et  al.  2007).  The  corroboration  of  these 
independent  reports  of  miRNA  changes  in  clinical  prostate  cancer  by  our  results  emphasize  that 
decrease  of  these  five  miRNAs  is  a  reliable  signature  for  progression  of  prostate  cancer.  In 
addition,  profiling  of  miRNAs  in  clinical  samples  suffers  from  the  contamination  of  the  cancer 
cells  by  stromal  and  non-malignant  epithelial  cells,  so  that  any  changes  noted  could  be  due  to 
changes  in  the  non-cancer  portion  of  the  specimen.  Our  results,  however,  clearly  indicate  that  the 
decrease  of  these  miRNAs  in  the  prostate  epithelial  cells  themselves  contributes  to  the  cancer 
progression. 

A  recent  report  also  attempted  to  identify  miRNAs  implicated  in  androgen 
dependence/independence  in  comparison  of  cell  lines  (Lin  et  al.  2008).  They  observed  several 
deregulated  miRNAs  between  the  pairs  of  cell  lines,  but  comparison  between  the  published  data 
and  our  data  revealed  that  only  miR-19b  was  consistently  down-regulated  in  androgen- 
independent  cell  lines  in  the  two  studies.  This  inconsistency  may  be  due  to  the  diversity  of 
mechanisms  by  which  prostate  cancer  cells  become  androgen-independent  (Pienta  and  Bradley 
2006).  Nonetheless,  the  discrepancy  was  still  seen  when  we  compared  our  data  on  the 
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LNCaP/C4-2  pair  with  the  published  data  on  the  LNCaP/C4-2B  pair.  A  possible  explanation 
could  be  in  the  difference  between  C4-2  and  its  bone  metastatic  derivative  C4-2B  (Thalmann  et  al. 
1994),  or  experimental  variables  such  as  passage  numbers  of  LNCaP  cells  (Lin  et  al.  2003).  It  is 
worth  noting  that  in  (Lin  et  al.  2008),  after  the  initial  array  (not  on  LNA  microarrays),  the 
induction  of  miR-184  and  repression  of  miR-146a  during  prostate  cancer  progression  was 
confirmed  not  by  miRNA-specific  biochemical  assays  on  the  cell  lines  but  by  in-situ- 
hybridization  to  cancer  tissues. 

miR-221,  -196b,  -125b  and  let-7  were  moderately  repressed  by  androgens  in  LNCaP 
cells,  though  we  do  not  know  whether  this  repression  is  a  primary  effect  of  androgen  or  due  to 
secondary  changes  after  the  activation  of  the  androgen-driven  gene  expression  program.  Contrary 
to  our  expectation,  we  did  not  find  any  miRNA  that  was  robustly  induced  by  androgens.  A  recent 
report  indicated  that  miR-125b  is  induced  by  androgens  and  that  miR-125b  is  up-regulated  in 
androgen  independent  CDS  cell  lines  relative  to  LNCaP  cells  (Shi  et  al.  2007).  Our  data  from 
multiple  independent  assays  unequivocally  showed  that  miR-125b  is  repressed  by  androgens  (Fig 
1C-D  and  Fig  S4)  and  significantly  down-regulated  in  the  androgen-independent  cell  lines  C4-2 
(Fig  1A  and  S4),  PC3,  and  DU145  (data  not  shown).  Of  note,  all  these  cell  lines  express  very  low 
level  of  AR  relative  to  CDS  cells  (Wu  et  al.  1994;  Mitchell  et  al.  2000;  Alimirah  et  al.  2006),  so 
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that  miR-125b  induction  noted  in  CDS  cells  could  be  related  to  level  of  expression  of  AR. 

The  differences  in  levels  of  specific  miRNAs  in  response  to  environmental  changes  (e.g. 
addition  of  androgen)  or  during  cancer  progression  appear  to  be  smaller  in  magnitude  than  the 
>100  fold  differences  seen  when  comparing  different  tissues.  However,  our  careful  comparison  of 
miRNA  profiles  by  all  the  currently  available  techniques  reveals  that  the  smaller  differences  seen 
in  response  to  environmental  changes  or  during  cancer  progression  can  be  identified  and  must  be 
checked  by  multiple  methods  before  acceptance. 

The  significant  under-representation  of  specific  miRNAs  from  androgen-depleted  cells  in 
the  sequencing  approach  compared  to  the  hybridization  based  methods  was  cautionary.  In  case  of 
miR-21,  the  explanation  lies  in  an  inefficiency  of  3 ’adaptor  ligation  in  the  androgen-depleted 
cells,  suggesting  a  modification  at  or  near  the  3 ’end  of  miR-21.  In  plants,  it  has  been  reported  that 
methylation  on  the  ribose  of  the  3’  nucleotide  contributes  to  stabilization  of  miRNA  (Yu  et  al. 
2005).  However,  such  a  modification  in  animal  miRNAs  has  not  been  reported  yet. 

One  advantage  of  the  deep  sequencing  technique  is  that  it  allowed  us  to  discover  novel 
non-micro-small  RNAs  (nmsRNAs).  The  sheer  diversity  and  abundance  of  the  nmsRNAs 
suggests  that  they  may  become  very  important  if  one  can  implicate  them  in  specific  biological 
functions.  The  two-fold  increase  in  the  abundance  of  nmsRNAs  upon  androgen-depletion  of 
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LNCaP  is  interesting  and  unexplained.  Some  of  the  nmsRNAs  were  mapped  to  genomic  repeats, 


suggesting  that  these  are  human  counterpart  of  rasiRNA  in  Drosophila.  However,  these  novel 
small  RNAs  are  18-20  nucleotides  long  and  are  smaller  than  Drosophila  rasiRNAs  that  range 
from  23-26  nucleotides  (Aravin  et  al.  2003).  A  subset  of  the  novel  small  RNAs  are  likely  to  be 
processed  from  abundant  noncoding  RNAs,  such  as  snoRNA  or  tRNA.  These  may  not  simply  be 
a  byproduct  of  the  degradation  process  because  they  were  of  a  definite  size  and  derived  from  a 
specific  small  portion  of  the  parental  RNA.  Taken  together  with  the  signal  intensity  in  Northern 
hybridization,  their  cloning  frequency  indicates  that  some  of  these  novel  nmsRNAs  are  as 
abundant  as  miRNAs. 

In  summary,  we  obtained  a  list  of  miRNAs  that  differ  in  expression  level  between  an 
androgen-independent  and  an  androgen-dependent  prostate  cancer  cell  line  and  upon  addition  of 
androgens  to  an  androgen-dependent  cell  line.  The  decrease  of  five  miRNAs  in  the  cancer  cells 
appears  to  be  a  definitive  signature  of  prostate  cancer  progression  that  can  be  used  for  prognosis. 
The  roles  of  these  miRNAs  in  prostate  cancer  biology  and  the  mode  of  their  regulation  remain  to 
be  determined.  Since  relatively  little  is  known  about  miRNAs  in  prostate  cancer  biology,  we 
believe  that  our  list  of  miRNAs  with  well  verified  changes  in  expression  level  will  be  a  good 
starting  point  to  investigate  the  roles  of  individual  miRNAs  in  prostate  epithelial  biology.  We  also 


Lee,  Yong  Sun 


21 


report  for  the  first  time  many  novel  small  RNAs,  some  of  which  are  more  abundant  than  most 


miRNAs.  The  biogenesis,  regulation  of  expression,  and  biological  function  of  this  huge  universe 
of  non-micro-short  RNAs  will  be  an  interesting  area  of  inquiry.  Finally,  the  study  underlines  the 
extensive  variability  that  could  arise  from  methods  of  miRNome  profiling  and  emphasizes  the 
need  to  corroborate  the  small  changes  in  miRNA  expression  often  reported  in  the  Literature  by 
independent  methods  of  miRNA  quantitation. 
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Materials  and  Methods 


Cell  culture 

LNCaP  and  C4-2  cell  lines  were  maintained  in  RPMI-1640  medium  supplemented  with 
10%  fetal  bovine  serum.  For  androgen  depletion,  LNCaP  cells  were  cultured  for  four  days  in  a 
phenol  red  free  RPMI-1640  with  10%  charcoal/dextran  treated  FBS  (Hyclone).  4,5  - 

dihydrotestosterone  (DHT)  was  purchased  from  Sigma-Aldrich  and  was  added  at  10  nM.  R1881 
was  from  Perkin-Elmer. 

RNA  isolation  and  miRNA  measurements 

Total  RNA  isolation  and  miRNA  analyses  (RNase  protection  assay  and  Northern 
hybridization)  were  performed  as  previously  described  (Kim  et  al.  2006;  Lee  and  Dutta  2007). 
Microarray  profiling  of  miRNAs  was  performed  with  miRCURY™  LNA  array  (v.9.2;  Exiqon). 

qPCR  measurement  of  miRNAs  was  performed  with  Ncode  SYBR  GreenER  miRNA 
qRT-PCR  kit  (Invitrogen)  according  to  the  manufacturer’s  instructions  and  with  ABI 7300  real 
time  PCR  system.  Each  value  is  an  average  of  triplicate  samples.  miRNA  signals  were 
normalized  to  those  of  U6  small  nuclear  RNA  (snU6)  as  a  control.  snU6  primer  sequence  is  5’- 
ctgcgcaaggatgacacgca-3’ . 
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Cloning  of  small  RNA  and  454  deep  sequencing 


Isolation  of  small  RNA,  ligation  of  adaptors,  cDNA  synthesis,  PCR  amplification,  and 
concatamerization  were  performed  as  described  in  (Lau  et  al.  2001)  with  minor  modifications. 
Concatamerized  DNAs  were  subject  to  454  deep  sequencing  (VBI  Core  Lab  at  Virginia 
Bioinformatics  Institute  in  Virginia  Tech.). 

Transfection  of  miRN A  duplex  and  cell  proliferation  assay 

siRNA-like  duplexes  containing  miRNA  sequence  at  one  strand  (Hutvagner  and  Zamore 
2002)  were  synthesized  Invitrogen  Life  Technologies  and  used  as  a  miRNA  mimic,  with  GL2 
duplex  as  a  negative  control  (Elbashir  et  al.  2001).  Transfection  into  C4-2  cells  was  performed 
with  Lipofectamine  2000  reagent  (Invitrogen  Life  Technologies)  at  a  final  concentration  of  20 
nM  of  each  siRNA  duplex.  After  72  hrs,  BrdU  incorporation  was  measured  as  previously 
described  (Machida  et  al.  2006)  and  was  normalized  to  cell  density  measured  by  MTT  assay 
(CellTiter  96  non-radioactive  cell  proliferation  assay  kit  from  Promega). 
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Supplemental  Data 


Supplemental  Tables  and  Figures  are  included  as  a  single  pdf  file  (Lee  et  al_RNA_supple.pdf). 
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Figure  legends 


Figure  1.  Expression  of  miRNAs  measured  by  Northern,  RNase  protection  assay,  and  qPCR. 
A  and  C.  qPCR  was  performed  as  indicated  in  Experimental  Procedures.  A  miRNA  value  was 
normalized  to  that  of  snU6.  The  average  and  standard  deviation  from  triplicate  samples  are 
indicated.  miRNAs  are  arranged  from  top  to  bottom  according  to  the  fold  change.  The  qPCR 
values  of  miRNAs,  together  with  their  other  results,  are  tabulated  in  Tables  S2  (panel  A;  C4-2 
versus  LNCaP)  and  Table  S3  (panel  C;  +DHT  versus  -DHT).  In  C4-2/LNCaP  comparison,  the 
higher  value  was  set  as  1  (panel  A);  while  -DHT  value  was  set  as  1  in  +DHT/-DHT  comparison 
(panel  C).  miRNAs  with  an  asterisk  are  also  measured  in  panel  B. 

B.  RNase  protection  assays  (top  panels)  and  Northern  hybridization  (bottom  panels)  of  selected 
miRNAs.  In  case  of  Northern  hybridization  of  let-7,  we  used  a  mixture  of  probes  against  let-7b 
and  -7e,  the  two  most  divergent  let-7  variants,  to  cover  as  many  kinds  of  paralogous  let-7s  as 
possible  (Lee  and  Dutta  2007).  miRNAs  with  an  asterisk  are  also  measured  in  panel  A  or  C. 

D.  qPCR  assays  of  selected  miRNAs  from  panel  C  after  treating  androgen  analog  R1881  at 
indicated  concentrations.  The  value  of  no  R1881  (-DHT  condition;  plain  bar)  is  set  as  1. 
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Figure  2.  The  proliferation  of  C4-2  after  transfection  of  miR-19b,  -99a,  -99b,  -100,  and  - 


125b  in  the  presence  or  absence  of  androgen. 

A.  Upon  transfection  of  indicated  miRNA(s)  or  GL2  (Experimental  Procedures),  C4-2  cells  were 
cultured  in  untreated  serum  or  charcoal  stripped  serum.  After  72  hrs,  BrdU  incorporation  was 
measured  and  normalized  to  cell  density  from  MTT  assay  (y-axis).  The  average  and  standard 
deviation  from  triplicate  samples  are  shown.  “miRs”  indicates  equimolar  mixture  of  the  five 
miRNAs  in  panel  B.  The  value  of  GL2  in  each  condition  is  set  as  1  (panel  A). 

B.  Transfection  of  individual  miRNA  mimics.  Each  value  in  charcoal- stripped  serum  is  again 
normalized  to  that  of  untreated  serum  and  the  relative  value  is  shown.  GL2  =  1 . 

Figure  3.  3’  adaptor  ligation  to  miR-21  from  -DHT  is  impaired 

Measurement  of  various  cloning  intermediates.  Steps  i  and  ii;  qPCR  with  total  RNA  (step  i)  or 
gel-isolated  17-26  mer  small  RNA  (sRNA;  step  ii).  Steps  iii  and  iv  are  qPCR  of  cloning 
intermediates.  qPCR  assays  after  control  ligation  of  5 ’adaptor  directly  to  sRNA  are  shown  in  a 
separate  step  vii.  The  PCR  primer  pairs  are;  a  miRNA  sense  primer  and  the  3’  adaptor  antisense 
primer  (step  iii),  or  a  miRNA  antisense  primer  and  the  5’  adaptor  primer  (step  iv  and  vii).  miR-21 
values  are  normalized  to  that  of  miR-200b  (y-axis).  All  values  are  an  average  of  triplicate 
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samples.  For  step  v,  40  nanogram  of  each  concatamer  (the  last  intermediate  in  the  cloning 
procedure)  was  subjected  to  Southern  hybridization  with  a  probe  specific  to  miR-21  (left  panel). 
A  mixture  of  probes  against  let-7b  and  -7e  was  used  as  a  control  (right  panel).  In  step  vi,  the 
cloning  frequency  of  miR-21  relative  to  let-7  is  shown. 


Table  legends 

Table  1.  Summary  of  454  deep  sequencing  data 

Numbers  of  total  and  unique  sequences  from  each  sample  (C4-2,  LNCaP,  +DHT  and  -DHT;  see 
text  for  details)  are  shown.  The  unique  sequences  are  classified  according  to  cloning  frequency. 
miRNA  sequences  are  extracted  from  the  total  sequences,  yielding  numbers  of  total  and  unique 
miRNA  sequences. 

Table  2.  List  of  selected  miRNAs  in  prostate  cancer  cells 

miRNAs  of  interesting  expression  pattern  were  selected  according  to  the  relative  expression 
levels  estimated  from  qPCR,  RNase  protection  assays,  and  Northern  hybridization  (Fig  1  and  S4). 
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The  approximate  relative  quantities  are  indicated  by  the  number  of  from  -/+  (the  lowest 
expression)  to  five  +’s  (the  highest  expression) 

Table  3.  Novel  small  RNAs 

A  list  of  non-micro-short  RNA  (nmsRNA)  that  were  cloned  abundantly  (>100  times  from  at  least 
one  sample).  Some  sequences  map  to  a  known  non-coding  RNA  or  genomic  repeat.  If  BLAST 
maps  a  sequence  to  more  than  three  loci  in  the  human  genome  or  ESTs,  one  representative 
coordinate  is  shown.  Another  BLAST  search  was  done  against  non-coding  RNA  database 
(http://www.noncode.org/)  and  the  results  are  shown  in  the  “note”  column. 
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Table  SI.  Sequencing  frequency  and  LNA  -microarray  data  of  all  miRNAs. 

Comprehensive  list  of  454  deep  sequencing  data  and  LNA-microarray  data.  Cloning 
frequencies  (raw  value)  for  each  miRNA  are  indicated.  Each  cloning  frequency  of  a 
miRNA  was  divided  by  the  total  number  of  sequence  of  the  corresponding  sample  (see 
Table  1;  for  example,  192,025  for  LNCaP  and  191,970  for  C4-2)  and  then  the 
normalized  values  were  used  to  calculate  ratio  of  C4-2/LNCaP  or  +DHT/-DHT  (for 
nomenclature,  see  the  text).  Each  value  in  microarray  data  is  ratio  of  median  values 
from  LN  CaP/C 4-2  or  +DHT/-DHT  pair.  "NA”  indicates  insignificant  or  poor  detection 
of  a  miRNA  in  microarray,  probably  due  to  low  expression  of  the  miRNA.  In  cases  of 
some  miRNAs  that  are  cloned  in  our  experiment  but  not  present  in  miRCURY  ™  LNA 
array  (v.9.2;  Exiqon),  they  are  indicated  as  "N  P". 

An  asterisk  was  given  to  miRNAs  of  interesting  expression  profile  which  include;  1) 
miRNAs  with  a  two-fold  or  higher  change  in  LNA-microarray,  or  2)  miRNAs  with  a 
change  by  more  than  33%  in  cloning  ratio  among  the  abundantly  (>30  in  either  one 
sample)  cloned  ones.  A  fraction  of  these  miRNAs  were  measured  by  miRNA  -specific 


assays  and  the  results  were  summarized  in  Table  S2  and  S3. 


Table  S2.  miRNAs  that  change  with  androgen  independence. 


miRNAs  were  selected  from  LNCaP  and  C4-2  based  on  changes  in  expression  in 
microarray  data  and  454  deep  sequencing  data  (T able  SI).  The  expression  of  each 
miRNA  was  measured  by  mi R N A -specific  assays  (Fig  1) 

Table  S3.  miRNAs  that  change  with  androgen  addition 

miRNAs  were  selected  from  +DHT  and  -DHT.  All  other  descriptions  are  same  as  Table 
S2. 

T able  S4.  Analyses  of  unique  sequences  that  are  not  miR  N  As. 

Out  of  unique  sequences,  those  that  do  not  match  miRNA  sequences  were  further 
analyzed  as  depicted  in  Fig  S2.  The  unique  sequences  are  subject  to  BLAST  search 
against  collection  of  rRNA,  tRNA,  snoRNA,  piRNA,  or  mRNA  (refseq)  sequences. 

Figure  SI.  RT-PCR  of  PSA  (prostate  specific  antigen)  and  p-actin 

A.RT  (reverse  transcri  ption)-PC  R  was  performed  with  the  total  RNAsfrom  LNCaP,  - 

DHT,  and  +DHT  to  measure  PSA  and  p-actin  control.  Each  amplified  fragment  with  its 


expected  size  in  bp  (base  pairs)  is  indicated.  The  primers  for  PSA  amplification  are  5'- 


attcccaaccctggcaggtg-3'  and  5' -tcgtagagcgggtgtgggaa-3' .  p-actin  primers  were 
previously  described  in  (Lee,  Y .  S.  &  Dutta,  A.  (2007)  Genes  Dev  21, 1025-1030). 

B.  PS  A /p-actin  values  from  real  time-qPCR  are  shown. 

Figure  S2.  Pipeline  of  analyses  of  454  deep  sequencing  data 

A  schematic  of  the  454  sequencing  data  analyses  is  depicted.  See  the  text  for  details. 

Figure  S3.  Quality  control  of  total  RNA  byethidium  bromide  staining 

Total  RNAs  are  subject  to  gel  electrophoresis  and  visualized  by  ethidium  bromide 
staining.  Left  panel',  1%  denaturing  agarose  gel  with  28S  and  18S  ribosomal  RNA  bands 
indicated .  Right  panel,  15%  denaturing  polyacrylamide  gel  with  the  positions  of  two 
dyes  (bromophenol  blue  and  xylene  cyanol  FF)  marked. 

Figure  S4.  M  easurement  of  mi R -125b  in  the  four  samples  and  after  R1881 
treatment. 

A.  N  orthern  hybridization  was  performed  to  measure  miR-125b  after  treatment  of 
indicated  concentrations  of  androgen  analog  R1881  (lanes  1-6)  and  in  the  four  samples 
(lanes  7-10).  5S  rRN  A  is  a  loading  control.  The  blot  was  exposed  to  phosphorimager 


Storm  860  (M  olecular  Dynamics)  and  each  band  was  quantified  with  ImageQuant  5.2 


software.  The  intensity  of  miR-125b  was  normalized  to  that  of  5S  rRNA.  "-DHT" 
samples  in  each  set  (lane  1  and  7)  were  set  as  1. 

B. qPCR  assays  to  detect  miR-125b.  Results  from  Fig  ID  (lane  1-6)  and  Fig  1A  and  C 
(lanes  7-10)  were  re-drawn  for  direct  comparison  with  Northern  data  in  panel  A. 

C.  PSA/p-actin  values  from  real  time-qPCR  after  R1881  treatment. 

D.  miR-98/snU6  measured  by  qPCR,  to  rule  out  global  miRNA  reduction  after  R 1881 
treatment. 

Figure  S5.  qPCR  assays  to  compare  C4-2/LNCaP  with  +DHT/-DHT 

Selected  miRNAsin  Fig  1  were  measured  by  qPCR  assays  for  side-by-side  comparison 
between  C4-2/LNCaP  (left  panel)  and  +DHT/-DFIT  (right  panel).  Some  of  the  data  was 
reclaimed  from  Fig  1A  and  C  for  ease  of  comparison. 

Figure  S6.  qPCR  assays  to  measure  miRN  As  underrepresented  in-DHT  in  cloning 
deep  sequencing 

Selected  miRN  As  in  Table  S3  were  measured  by  qPCR  assays.  M  icroarray  and  deep 
sequencing  values  of  each  miRNA  are  shown  in  the  left.  All  other  descriptions  are  same 


as  Fig  1C. 


Figure S7.  M easurement of  miR-21  in  +DHT  and  -DHT 


A.  Luciferase  assays  were  performed  as  described  (Lee,  Y.  S.  &  Dutta,  A .  (2007)  Genes 
Dev  21, 1025-1030).  pRLCM  V-miR21as  is  a  plasmid  containing  a  sequence  perfectly 
complementary  to  miR-21  downstream  from  the  Renilla  luciferase  ORF  in  the  parental 
vector,  pRL-CM  V.  pRLCM  V-miR21as  (or  pRL-CM  V  vector  control),  together  with 

pG L 3-C ontrol  vector  (Promega  Corp.),  were  transfected  into -DHT  and  +DHT  using 
Lipofectamine  2000  reagent  (I  nvitrogen).  Renilla  luciferase  activity  was  normalized  to 
the  f i  refly  ( Photinus  pyralis)  I  uciferase  activity  from  the  co-transfected  pG  L  3-C  ontrol 
vector,  y-axis  indicates  ratio  of  pRLCM  V-miR21as  to  pRL-CM  V.  The  average  and 
standard  deviation  of  triplicate  samples  are  indicated. 

B.  miR-21  sequences  were  cloned  with  variable  ends,  probably  due  to  a  flexibility  of 
Dicer  or  Drosha  processing.  Among  these  sequences  (490  times  in  -  DHT  and  2,053 
times  in  +DHT  as  a  whole)  are  shown  the  two  most  abundantly  cloned  miR-21 
sequences  (italicized)  and  their  cloning  frequencies  in  the  454  data.  At  the  bottom, 
mature  miR-21  sequence  (capital  letters)  and  its  partial  flanking  sequence  from  the 
precursor  (lower  case  letters)  are  shown  from  the  database  (miRBase; 


http://microrna.sanqer.ac.uk/). 


C.  In  an  attempt  to  detect  3'end  modification  in  the  miR-21  that  may  interfere  with  the 


3'  adaptor  ligation,  another  N orthern  hybridization  was  performed  on  an  acrylamide  gel 
of  higher  resolution  (18%  denaturing  acrylamide  sequencing  gel).  For  comparison,  a  26- 
meroligodeoxyribonucleotide  containing  miR-21  sequence  was  phosphorylated  at  the 
5'end  by  polynucleotide  kinase  (New  England  Biolabs),  was  run  in  parallel  with  the 
unphosphorylated  form,  and  was  subject  to  hybridization  to  a  probe  against  miR-21.  As 
seen  in  right  panel,  a  single  phophate  group  was  resolved  in  this  gel.  In  the  miR-21 
lanes  (left  two  lanes),  we  observed  two  bands  that  are  presumably  22-  and  23-mer 
which  represent  the  two  most  abundantly  cloned  forms  of  miR-21  (panel  B).  We  failed 
to  detect  any  aberrant  migration  of  miR-21  in  -DHT  compared  to  +DHT,  suggesting 
that  any  modification  at  the  3'end  is  too  small  to  be  resolved  by  this  technique. 

Figure  S8.  Northern  hybridization  detection  of  novel  small  RNAs 

Northern  hybridization  was  performed  to  detect  selected  non-micro-small  RNAs 
(nmsRNAs).  Two  fragments  from  U  3  snoRN  A  (ID  000134_1499  and  ID  004265_1575) 
and  a  fragment  from  tRN  A  (ID  001257  0208)  were  chosen  from  Table  3.  A  discrete 
band  of  expected  size  (marked  by  arrowheads)  confirmed  the  existence  of  the  small 


RN  A  at  a  measurable  steady-state  level.  The  very  strong  bands  on  top  of  the  gel  are 


most  likely  due  to  full  length  snoRNA  U3  ortRNA  molecules.  We  detected  additional 


bands  of  intermediate  size  whose  intensity  is  roughly  proportional  to  that  of  the  small 
RNA,  raising  the  possiblity  that  the  novel  smalIRNAs  were  derived  from  snoRNA  U3 
ortRNA  through  a  stepwise  processing,  with  the  intermediate  sized  bands  being 
processing  intermediates. 
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Abstract 


New  types  of  small  RNAs  distinct  from  microRNAs  (miRNAs)  are 
progressively  being  discovered  in  various  organisms.  In  order  to  discover  such  novel 
small  RNAs,  a  library  of  17-26  base  long  RNAs  was  created  from  prostate  cancer  cell 
lines  and  sequenced  by  ultra-highthroughput  sequencing.  A  significant  number  of  the 
sequences  are  derived  from  precise  processing  at  the  5'-  or  3'-end  of  mature  or  precursor 
tRNAs,  to  form  three  series  of  tRFs,  the  -5,  -3  and  -1  series.  These  sequences, 
collectively  named  as  tRFs  (tRNA-derived  RNA  Fragments),  constitute  a  class  of  short 
RNAs  that  are  second-most  abundant  to  miRNAs.  Northern  hybridization  and  qRT-PCR 
assays  independently  measured  the  levels  of  at  least  16  tRFs.  tRF-1001,  derived  from 
the  3'-end  of  a  Ser-TGA  tRNA  precursor  transcript  that  is  not  retained  in  the  mature 
tRNA,  was  further  investigated.  tRF-1001  is  expressed  highly  in  a  wide  range  of  cancer 
cell  lines  but  much  less  in  tissues  and  its  expression  in  cell  lines  was  tightly  correlated 
with  cell  proliferation.  siRNA-mediated  knockdown  of  tRF-1001  impaired  cell 
proliferation  with  the  specific  accumulation  of  cells  in  G2,  phenotypes  that  were 
reversed  by  co-introducing  a  synthetic  2'-0-methyl  tRF-1001  oligoribonucleotide 
resistant  to  the  siRNA.  Although  tRF-1001  is  localized  in  the  cytoplasm  and  similar  in 
size  to  miRNAs  and  siRNAs,  it  did  not  repress  a  target  mRNA,  suggesting  that  it  does 
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not  enter  into  RISC  or  miRNP  complexes.  Our  data  suggest  that  tRFs  are  not  random 
by-products  of  tRNA  degradation  or  biogenesis  but  an  abundant  and  novel  class  of  short 
RNAs  with  precise  sequence  structure  that  have  specific  expression  patterns  and 
specific  biological  roles. 
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Introduction 


Small  non-coding  regulatory  RNAs  have  emerged  as  important  players  in 
diverse  aspects  of  biology.  The  most  extensively  studied  among  them  are  microRNAs 
(miRNAs)  that  regulate  target  mRNA  post-transcriptionally  in  animals  and  plants 
(reviewed  in  (Bartel  2004;  Du  and  Zamore  2005;  Lee  and  Dutta  2008)).  Another  class 
of  small  RNA  closely  related  to  miRNAs  is  small  interfering  RNA  (siRNA),  a  key 
intermediate  in  the  RNA  interference  pathway  (Elbashir  et  al.  2001b).  Although  siRNAs 
are  a  popular  tool  to  knockdown  genes  in  mammalian  cells,  the  existence  of  endogenous 
siRNAs  is  still  questionable  in  mammals.  This  is  in  contrast  to  plants  where  naturally 
occurring  siRNAs,  tasiRNA  (trans-acting  siRNA)  and  natsiRNA  (natural  antisense 
transcripts  siRNA),  have  been  discovered  (reviewed  in  (Chapman  and  Carrington 
2007)). 

A  number  of  other  small  RNAs  have  also  been  identified  in  various  organisms 
(reviewed  in  (Farazi  et  al.  2008)),  including  tncRNA  (tiny  noncoding  RNAs),  21U-RNA, 
rasiRNA  (repeat  associated  siRNA),  hcRNA  (heterochromatic  small  RNA),  scnRNA 
(scan  RNA),  and  piRNA  (piwi-interacting  RNA).  Although  some  of  these  have  been 
shown  to  play  biological  roles  in  processes  such  as  heterochromatin  formation  and 
transposon  silencing,  most  of  them  are  recently  discovered  and  not  fully  understood. 
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We  believe  that  additional  small  RNAs  await  discovery  expecially  through  the 
application  of  high-throughput  sequencing  technology  (reviewed  in  (Morozova  and 
Marra  2008)). 

In  this  paper,  we  analyzed  the  global  expression  profile  of  small  RNA  (small 
RNome)  in  human  prostate  cancer  cell  lines  by  454  deep  sequencing.  Next  to  miRNAs, 
the  most  abundant  group  of  small  RNA  was  fragments  derived  from  tRNA  (termed  here 
as  “tRNA-derived  RNA  fragments;  tRFs”).  Besides  their  well-known  role  in  translation, 
tRNAs  have  been  shown  to  function  in  other  cellular  processes  such  as  reverse 
transcription,  porphyrin  biosynthesis,  etc  (reviewed  in  (Giege  2008)).  Here,  we  identify 
a  new  aspect  of  tRNA  biology.  The  characteristic  and  abundant  expression  of  individual 
tRFs  reported  here  as  well  as  their  precise  sequence  composition,  indicates  that  tRFs  are 
not  random  degradation  intermediates  during  tRNA  biogenesis  and  turnover,  but 
constitute  a  novel  class  of  small  RNA  that  deserve  further  scrutiny.  One  member  of  tRFs 
is  shown  to  be  required  for  cell  viability. 


5 


Results 


Cloning  and  identification  of  tRF-5,  -3  and  -1  series.  Small  RNAs  of  17-26 
nucleotides  were  cloned  from  prostate  cancer  cells  and  sequenced  by  454  deep 
sequencing  technology  (Lee  et  al,  in  preparation).  More  than  half  of  the  sequence  reads 
did  not  match  miRNAs.  Among  these  non-micro-small  RNAs  (nmsRNAs),  we  focused 
on  30  sequences  that  were  abundantly  cloned  (>300  times;  Table  1).  In  comparison,  a 
well-known  miRNA  with  a  demonstrated  biological  function,  miR-18a,  was  cloned  473 
times  while  a  highly  abundant  miRNA,  let- 7a,  was  cloned  17,856  times.  Of  695 
miRNAs  in  the  database  (miRBase;  http://microma.sanger.ac.uk/),  a  majority  of  them 
(635  miRNAs)  were  cloned  <200  times. 

Each  of  these  abundant  nmsRNA  sequences  were  subjected  to  BLAST  search 
of  the  human  genome  sequence  and  was  further  identified  by  examining  its  genomic 
locus  on  the  UCSC  genome  browser.  28  sequences  were  mappable  to  the  human 
genome,  in  many  cases  at  multiple  loci  (Fig  1  A).  Six  of  these  are  located  within  known 
transcripts  and  snoRNA  (small  nucleolar  RNA),  from  which  the  small  RNAs  are  likely 
to  be  derived.  Of  note,  more  than  half  (17  out  of  30)  of  the  nmsRNA  sequences  mapped 
to  tRNA  loci  (Fig  1A  and  Table  1).  Given  that  a  cloning  frequency  approximately 
reflects  the  intracellular  level  of  a  small  RNA,  these  tRNA-related  small  RNAs  are 


expected  to  be  more  abundant  than  most  miRNAs. 


Primary  transcripts  (pre-tRNA)  for  mature  tRNA  are  transcribed  by  RNA 
polymerase  III  and  contain  a  5'  leader  sequence  and  a  3'  trailer  sequence  that  are 
trimmed  during  tRNA  maturation  (reviewed  in  (Mori  and  Marchf elder  2001;  Hopper 
and  Phizicky  2003)).  For  most  eukaryotic  tRNAs,  a  CCA  sequence  is  enzymatically 
added  to  the  3'-end  of  the  trimmed  tRNA  intermediate.  Among  the  17  tRNA-related 
small  RNAs,  five  and  eight  sequences  are  precisely  aligned  to  the  5'-  and  3'-ends  of 
mature  tRNA,  respectively  (Fig  1  and  SI,  tRF-5  and  tRF-3  series).  The  3'-end  sequence 
of  all  the  eight  tRF-3  RNAs  includes  the  CCA  that  is  post-transcriptionally  added  to  the 
3'-end  of  tRNAs  during  tRNA  maturation.  We  postulate  that  the  tRF-5  and  tRF-3  series 
of  small  RNAs  are  derived  from  cleavage  of  mature  tRNAs  after  they  have  undergone 
5'-/3'-end  trimming  and  CCA  addition.  The  precise  sequence  conservation  of  the  ends  of 
all  the  clones  strongly  suggests  that  a  precise  cleavage  event  is  necessary  for  the 
generation  of  these  tRFs. 

Four  small  RNAs  of  the  tRF-1  series  are  located  within  pre-tRNA,  3'  to  the 
mature  tRNA  sequence.  In  all  cases,  their  5'-ends  start  precisely  after  the  3'-ends  of 
mature  tRNA  (Fig  1  and  SI).  The  3'-ends  fall  within  DNA  sequences  of  5-6  consecutive 
Ts  (Fig  SI),  a  canonical  termination  site  of  RNA  polymerase  III  transcription 
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(Hagenbuchle  et  al.  1979;  Koski  and  Clarkson  1982;  Cozzarelli  et  al.  1983).  Most  likely, 


the  tRF-1  series  of  small  RNAs  have  been  released  from  pre-tRNAby  an 
endonucleolytic  cleavage  during  the  3'-end  maturation  of  tRNA. 

To  obtain  a  more  comprehensive  list  of  tRFs,  we  expanded  our  analysis  to  all 
small  RNA  sequences  cloned  more  than  5  times  in  the  library  (see  Table  SI  and  S2). 
Out  of  1541  nmsRNA  sequences  that  mapped  to  the  human  genome,  approximately 
40%  mapped  to  genomic  loci  containing  a  tRNA  and  flanking  25  nucleotides  (Table  S 1 
and  Fig  S2).  77%  of  these  were  identified  as  tRFs  by  their  precise  start  and  end  sites, 
while  the  remaining  -23%  match  tRNA  at  random  locations  and  were  not  further 
investigated.  The  fraction  of  tRF-1  series  (14.2%)  is  lower  than  that  of  tRF-3  (44.8%) 
or  -5  series  (41.0%). 

The  lower  frequency  of  tRF-1  clones  is  explained  by  the  hypothesis  that  tRF-1 
are  released  by  3'  cleavage  of  the  tRNA  from  the  pre-tRNA.  Of  the  622  tRNA  loci  in 
the  human  genome  (http://lowelab.ucsc.edu/GtRNAdb/)  the  distance  between  the  3'- 
end  of  the  mature  tRNA  and  the  termination  signal  of  the  RNA  polymerase  III  is 
variable  (Fig  S3).  Only  14.6%  of  the  tRNA  loci  could  give  rise  to  tRF-1  of  the  16-27 
base  size  range  that  was  included  in  our  cloning,  accounting  for  the  lower  yield  of  tRF- 


1  series. 


Thus,  tRF-5,  -3  and  -1  (Table  S2)  constitute  the  most  abundant  class  of  short 


RNAs  in  human  prostate  epithelial  cells,  second  only  to  miRNAs.  In  addition,  their 
precise  start  and  end  sites  at  or  near  the  tRNA  ends  strongly  suggest  that  the  tRFs  are 
derived  from  tRNA  cleavage  in  a  specific  manner. 

Expression  of  tRFs  validated  by  other  techniques.  Northern  hybridization 
detected  six  tRFs  tested  (Fig  S4A).  Not  only  did  we  detect  a  band  of  the  correct  size,  but 
other  bands  in  the  range  of  80-100  nucleotides  were  apparent  and  could  be  precursors  to 
the  final  tRFs.  For  tRF-1001,  -1002,  and  -1003,  the  larger  bands  are  likely  to  be  the  pre- 
tRNA  from  which  the  tRF  is  derived.  For  tRF-5  and  -3  series,  the  probe  also  detects  the 
corresponding  mature  tRNA,  which  is  present  at  a  much  higher  abundance  than  the  tRF. 
Below  the  mature  tRNA  bands,  there  are  additional  bands  of  intermediate  size  (35-60 
nucleotides  long)  that  could  correspond  to  the  tRNA  cleavage  products  that  were 
recently  reported  (Li  et  al.  2008;  Thompson  et  al.  2008). 

The  expression  level  of  some  tRFs  were  too  low  to  be  detected  by  Northern 
hybridization,  particularly  for  tRF-3  and  tRF-5  series,  where  a  strong  signal  from  the 
mature  tRNA  and  intermediate  bands  interferes  with  the  tRF  signal.  To  circumvent  this 
problem,  we  employed  qRT-PCR  detection  in  a  method  similar  to  the  detection  of 
miRNAs  (reviewed  in  (Schmittgen  et  al.  2008)).  We  gel-purified  the  RNA  of  17-26 
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nucleotides  prior  to  cDNA  synthesis  to  exclude  contaminating  signals  from  larger 
mature  tRNAs  or  other  RNAs.  qRT-PCR  for  12  tRFs  successfully  detected  all  of  them 
(Fig  S4B).  Interestingly,  each  tRF  exhibited  a  characteristic  expression  pattern  across 
the  four  RNA  samples  (see  Fig  S4  and  legend),  arguing  against  the  possibility  that  the 
tRFs  are  random  by-products  from  nonspecific  degradation  of  tRNA.  Such  non-specific 
degradation  would  have  resulted  in  a  similar  pattern  of  expression  of  all  tRFs  amongst 
the  four  RNA  preparations. 

Expression  of  tRF-1001.  tRF-1001  was  chosen  for  further  study,  as  it  is  the 
most  abundantly  cloned  tRF  and  was  clearly  detectable  by  Northern  hybridization.  tRF- 
1001  is  expressed  more  abundantly  in  cell  lines  than  in  tissues  (Fig  2 A;  upper  panel). 
Despite  the  greater  degradation  of  RNA  from  tissues,  as  indicated  by  the  ethidium 
bromide  (EtBr)  staining  (Fig  2A;  bottom  panel),  tRF-1001  levels  are  not  increased  in 
the  tissues,  reinforcing  the  notion  that  tRF-1001  is  not  a  byproduct  of  non-specific  RNA 
degradation.  The  higher  expression  of  tRF-1001  in  cancer  cell  lines  of  many  different 
lineages  suggests  that  tRF-1001  expression  may  be  related  to  cell  proliferation. 

Consistent  with  this,  tRF-1001  is  decreased  upon  serum  starvation  of  DU145 
prostate  cancer,  LNCaP  prostate  cancer  and  HCT116  colon  cancer  cell  lines  and  is 
restored  upon  serum  re-addition  (Fig  2B,  S5A-B).  tRF-1001  levels  are  also  reduced 


10 


when  cell-density  is  very  high  (Fig  S5C).  Therefore,  conditions  of  poor  cell 
proliferation,  as  with  serum  depletion  or  increased  cell  density,  are  associated  with  a 
down-regulation  of  tRF-1001.  In  contrast,  the  corresponding  mature  tRNA  level  stayed 
unchanged. 

tRF-1001  is  required  for  cell  proliferation.  The  best  argument  against  the 
tRFs  being  a  non-specific  byproduct  of  tRNA  metabolism  is  to  demonstrate  a  functional 
role  for  at  least  one  of  the  tRFs.  To  investigate  this,  we  knocked  down  tRF-1001 
through  an  siRNA  strategy.  We  designed  an  siRNA  duplex  that  targets  the  pre-tRNA 
such  that  it  anneals  to  part  of  the  tRF-1001  and  the  mature  tRNA  (without  the  CCA). 
The  siRNA  is  centered  at  +2/+3  position  of  tRF-1001  (Fig  3A).  Thus,  siRNA-mediated 
cleavage  of  the  pre-tRNA  at  the  center  of  siRNA  target  sequence  (Elbashir  et  al.  2001b), 
is  expected  to  cut  the  pre-tRNA  at  a  site  (arrowheads  in  Fig  3A)  incompatible  with  tRF- 
1001  generation.  Alternatively,  the  siRNA  could  sterically  interfere  with  the  cleavage 
event  necessary  for  tRF-1001  generation. 

In  this  experiment,  Northern  hybridization  was  not  suitable  for  measuring  tRF- 
1001  knock-down,  because  the  passenger  strand  of  si-tRFlOOl  shares  12  nucleotides 
with  tRF-1001  and  obscured  the  latter  (data  not  shown).  Instead,  we  took  advantage  of 
the  fact  that  qRT-PCR  assay  of  tRF-1001  will  fail  to  detect  the  DNA  passenger  strand  of 
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si-tRF  1 00 1  for  the  following  reasons.  First,  the  first-step  in  qRT-PCR  is  polyadenylation 


of  the  small  RNA  by  poly  (A)  polymerase  that  requires  3'  ribonucleotides  for  efficient 
poly(A)  synthesis  (Zhelkovsky  et  al.  1998).  Second,  the  3'-end  of  the  specific  primer 
used  in  the  second  step  of  qRT-PCR  anneals  to  the  tRF  but  not  to  the  passenger  strand  of 
si-tRF  1001.  The  reliability  of  qRT-PCR  in  measuring  tRF- 1001  was  demonstrated  by 
side-by-side  comparison  with  Northern  hybridization  (Fig  S6). 

qRT-PCR  assay  confirmed  the  reduction  of  tRF-1001  by  the  siRNA  (Fig  3B).  In 
addition  to  standard  qRT-PCR  procedure  for  tRFs  using  gel-isolated  small  RNA,  qRT- 
PCR  with  total  RNA  yielded  a  similar  result  (Fig  3B;  right  panel). 

Upon  tRF-1001  depletion,  cell  proliferation  was  impaired  as  indicated  by  a 
decrease  in  the  number  of  viable  cells  (MTT  assay  in  Fig  3C),  which  was  accompanied 
by  a  reduction  of  DNA  synthesis  (BrdU  incorporation  assay  in  Fig  3C)  and  an 
accumulation  of  cells  in  the  G2  phase  of  the  cell  cycle  (FACS  profile  in  Fig  3D).  A 
closer  observation  of  individual  cells  revealed  an  increase  in  the  population  of  large 
cells  with  big  nuclei  (Fig  S8).  Although  the  si-tRFlOOl  is  expected  to  decrease  the 
maturation  of  the  tRNA  associated  with  tRF-1001,  the  same  tRNA  (Ser-TGA)  is 
expressed  from  several  alternate  genomic  loci  different  from  the  tRF-1001  locus.  Indeed, 
Ser-TGA  tRNA  level  is  unchanged  after  si-tRFlOOl  treatment  (Fig  S7).  Thus,  the  cell- 
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proliferation  defects  induced  by  si-tRFlOOl  cannot  be  attributed  to  a  decrease  in  the 


corresponding  mature  tRNA . 

Most  importantly,  the  decrease  in  cell  proliferation  and  G2  accumulation 
induced  by  si-tRFlOOl  were  rescued  by  co-transfection  of  a  synthetic  tRF-1001  in  the 
form  of  a  2'-0-methyl  oligoribonucleotide  (Me-tRF-1001,  Fig  3C  and  D),  but  not  by  a 
similar  sequence  from  an  irrelevant  gene,  luciferase  (Me-GL2).  We  used  2'-0-methyl 
derivatives  because  this  modification  is  known  to  make  the  oligoribonucleotide 
chemically  stable  and  resistant  to  cellular  nucleases  (reviewed  in  (Lamond  and  Sproat 
1993)).  Because  the  sequence  of  Me-tRF-1001  overlaps  si-tRFlOOl  by  12  nucleotides 
(Fig  3A  and  S9A),  we  wanted  to  eliminate  the  possibility  that  the  rescue  was  because 
the  Me-  tRF-1001  anneals  to  the  siRNA  and  titrates  it  away  from  a  non-specific  cellular 
target.  The  Me-control  oligoribonucleotide  (Fig  S9A)  overlaps  with  another  12 
nucleotides  of  si-tRFlOOl,  but  is  different  from  the  tRF-1001  sequence.  Me-control  did 
not  rescue  the  growth  impairment  induced  by  si-tRFlOOl  (Fig  S9B),  strongly 
suggesting  that  the  rescue  of  cell  proliferation  by  Me-tRF-1001  was  because  it  acted  as 
a  tRF-1001  mimic. 

tRF-1001  does  not  act  as  an  si/miRNA.  tRF-1001  is  exclusively  localized  in 
the  cytoplasm  (Fig  4A).  Its  size  (19-mer)  and  intracellular  localization  raise  the 
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possibility  that  tRF-1001  may  act  like  an  siRNA  or  miRNA.  To  address  this  question  we 


designed  a  reporter  plasmid  containing  a  perfect  antisense  to  tRF-1001  downstream  of 
the  luciferase  ORF  (pRL-anti-tRFlOOl  in  Fig  4B).  If  endogenous  tRF-1001  is 
incorporated  into  the  siRNA/miRNA  pathway,  it  is  expected  to  suppress  the  expression 
of  luciferase  that  has  anti-tRF-1001  in  its  3'UTR. 

Reduction  of  endogenous  tRF-1001,  either  by  si-tRFlOOl  or  by  serum 
starvation,  did  not  increase  the  luciferase  activity  from  pRL-anti-tRFlOOl  (Fig  4C  and 
D).  Transfection  of  single-stranded  Me-tRF-1001,  a  functional  form  of  tRF-1001,  also 
did  not  inhibit  the  luciferase  expression  (Fig  4E).  As  a  technical  control,  tRF-1001 
synthesized  and  transfected  as  an  siRNA  duplex  significanly  inhibited  the  luciferase 
expression  of  pRL-anti-tRFlOOl  (Fig  4E),  Thus,  tRF-1001  will  act  as  an  siRNA  or 
miRNA  if  it  is  artificially  delivered  into  a  RISC  complex.  Endogenous  tRF-1001, 
however,  appears  to  not  be  in  a  RISC,  consistent  with  the  lack  of  any  evidence  that  it  is 
ever  part  of  an  RNA  duplex  :  1)  tRF-1001  did  not  form  a  stem- loop  secondary  structure 
with  its  flanking  genomic  sequence  and  2)  a  sequence  complementary  to  tRF-1001  was 
not  cloned  in  the  454  high-throughput  sequencing.  Collectively,  these  results  suggest 
that  tRF-1001  is  a  functionally  distinct  small  RNA  from  siRNA/miRNA. 
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Discussion 


We  report  that  tRFs  are  a  novel  class  of  small  RNAs  that  are  only  second  in 
abundance  next  to  miRNAs.  tRF-5  and  -3  are  derived  from  mature  tRNAs.  Although  it 
is  conceivable  that  as  many  of  these  tRFs  will  exist  as  there  are  mature  tRNAs,  we  did 
not  clone  all  the  possible  tRFs-5  and  -3.  Whether  this  is  due  to  our  limited  cloning  and 
sequencing  effort  or  due  to  biological  differences  in  the  processing  of  tRNAs  will  be 
determined  as  we  accumulate  more  reports  of  cloning  and  sequencing  of  short  RNAs. 

Of  the  622  tRNA  loci,  an  analysis  of  the  distance  between  the  3 '-end  of  the  mature 
tRNA  and  3'-termination  site  of  the  pre-tRNA  (Fig  S3)  indicates  that  fewer  tRNA  loci 
can  produce  tRF-1  series  of  17-26  nucleotides.  However,  it  is  entirely  possible  that 
future  cloning  of  longer  or  shorter  RNAs  will  yield  more  tRF-1  clones.  Like  miRNAs, 
we  expect  different  sets  of  tRFs  to  be  expressed  in  different  cell  lines  and  tissues  under 
different  growth  conditions  and  future  efforts  will  be  directed  at  expanding  the 
repertoire  of  tRFs. 

The  mechanism  of  biogenesis  of  the  tRFs  is  not  clear  yet,  even  though  the 
juxtaposition  of  the  start  and  end  points  of  tRFs  with  that  of  mature  tRNAs  suggests  that 
they  are  generated  from  cleavage  of  tRNA  or  pre-tRNA  by  specific  nucleases.  Two 
endonucleases  for  3'-trimming  of  pre-tRNA,  human  ELAC1  and  ELAC2  (Takaku  et  al. 
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2003),  are  the  best  candidates  for  generation  of  tRF-1  series.  It  is  intriguing  that  ELAC2 


was  originally  identified  as  a  candidate  prostate  cancer  susceptibility  gene  (Tavtigian  et 
al.  2001),  given  that  tRF-1001  is  highly  expressed  in  cancers  and  required  for  cancer 
cell  growth. 

Ranpirnase  or  Onconase  is  an  endonuclease  for  tRNAs  that  is  also  known  for  its 
cytotoxic  and  antitumor  activity  (Ardelt  et  al.  1991;  Saxena  et  al.  2002).  Ranpirnase 
belongs  to  the  pancreatic  RNase  family  that  has  several  members  in  humans  (Rybak  and 
Newton  1999).  Ranpirnase,  however,  cleaves  preferentially  between  the  two  G  residues 
of  UGG  sequence  in  vitro,  in  the  variable  arm  or  the  D-arm  regions  (Suhasini  and 
Sirdeshmukh  2006).  UGG  does  not  define  any  of  the  cleavage  sites  utilized  to  generate 
tRF-3  or  -5,  so  that  either  Ranpirnases  are  not  relevant  for  the  generation  of  these  tRFs 
or  one  of  the  Ranpirnases  has  a  different  cleavage  site  in  cells  than  what  has  been 
defined  in  vitro. 

The  correlation  of  tRF-1001  expression  with  cell  proliferation  raises  the 
possibility  that  tRF-1001  can  be  used  as  a  molecular  marker  for  proliferating  cancer 
cells.  The  utility  of  miRNAs  as  a  cancer  diagnostic  tool  has  been  well  described 
(reviewed  in  (Lee  and  Dutta  2008)).  miRNAs  exhibit  characteristic  expression  profile 
according  to  a  cancer  type,  stage,  and  other  clinical  variables  (Lu  et  al.  2005).  The 
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miRNAs  in  serum  or  plasma  are  in  a  special  form  that  renders  them  more  resistant  to 


RNase  degradation,  making  them  reliable  bio-markers  in  tissues  and  body -fluids 
(Mitchell  et  al.  2008).  The  sensitive  qRT-PCR  assays  for  miRNA  allow  their 
measurement  in  even  a  few  nanograms  of  total  RNA  from  formalin-fixed  paraffin- 
embedded  (FFPE)  specimens  (Xi  et  al.  2007).  All  the  advantages  of  miRNAs  as  a 
molecular  marker  are  applicable  to  tRF-1001,  especially  since  it  can  be  reliably 
measured  from  total  RNA  by  qRT-PCR  assays  (Fig  3B).  Since  the  expression  level  of 
tRF-1001  seems  to  respond  quickly  to  cell  proliferation  (Fig  2B  and  Fig  S5),  tRF-1001 
can  be  used  as  a  proliferation  marker  similar  to  the  Ki-67  protein  (Scholzen  and  Gerdes 
2000). 

The  tRF-1  series  are  of  particular  interest,  not  only  because  we  show  that  one  of 
them  is  required  for  cell  proliferation,  but  also  because  the  sequence  downstream  of  the 
tRNA  that  is  encoded  in  the  pre-tRNA  is  mostly  unique  among  different  tRNA  genes 
(Table  S2).  Mature  tRNAs  of  identical  or  highly  similar  sequences  are  repeated  at 
multiple  loci  in  the  genome  as  expected  from  the  fact  that  622  tRNA  loci  code  for 
tRNAs  containing  at  most  64  anticodons.  We  speculate  that  some  tRNA  genes  have 
additional  gene-specific  functions  by  generating  unique  tRFs  from  the  3'-flanking 
sequences,  perhaps  explaining  the  duplication  and  survival  of  so  many  tRNA  genes 
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during  evolution. 


Much  more  work  needs  to  be  done  to  investigate  how  many  tRFs  are  expressed, 
what  regulates  their  levels,  how  many  have  discernible  biological  roles  and  whether  the 
tRFs  exert  their  biological  functions  as  part  of  ribonucleoprotein  complexes.  Despite 
this,  however,  it  is  exciting  to  consider  the  obvious  parallels  between  the  current  status 
of  tRFs  with  what  was  known  about  miRNAs  soon  after  it  became  evident  that 
hundreds  of  different  miRNAs  are  present  in  mammalian  cells  (Lagos-Quintana  et  al. 
2001;  Lau  et  al.  2001;  Lee  and  Ambros  2001). 
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Materials  and  Methods 


Isolation  of  RN  A,  454  deep  sequencing,  and  measurement  of  tRFs 

Isolation  of  total,  nuclear,  and  cytoplasmic  RNA  from  cell  lines  was  performed  as 
described  in  (Lee  and  Dutta  2007).  Small  RNAs  of  17-26  nucleotides  were  cloned  and 
subjected  to  454  deep  sequencing  as  described  in  (Lee  et  al,  in  preparation).  qRT-PCR 
measurement  of  tRFs  and  miRNAs  was  performed  with  Ncode  SYBR  GreenER  miRNA 
qRT-PCR  kit  (Invitrogen  Corp.)  and  ABI  7300  real  time  PCR  system  (Applied 
Biosystems).  Unless  described  otherwise,  total  RNA  was  subjected  to  15%  denaturing 
acrylamide  gel  electrophoresis  for  isolation  of  small  RNA  (17-26  nucleotides  long)  prior 
to  qRT-PCR  steps.  The  qRT-PCR  primer  for  measuring  small  nuclear  RNA  U6  in  Fig. 

3B  is  5'-ctgcgcaaggatgacacgca-3'. 

Cell  culture  and  transfection  of  RNA  oligonucleotides 

Cell  lines  and  culture  conditions  are  described  in  Supplemental  Materials  and  Methods. 
siRN A  duplexes  were  synthesized  by  Invitrogen  Corp.  si-tRFlOOl  sequence  is  shown  in 
Fig  3A  and  siGL2  sequence  is  described  in  (Elbashir  et  al.  2001a).  2'-0-methyl 
modified  GL2  (5'-cguacgcggaauacuucga-3')  and  tRF-1001  (5'-gaagcgggugcucuuauuu- 
3')  were  synthesized  by  Dharmacon  RNA  Technologies.  Transfection  of  short  RNA  was 
performed  with  Lipofectamine™  RNAiMAX  reagent  (Invitrogen  Corp.)  using  33  nM 
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siRNA  duplex  and  67  nM  2'-0-methyl  oligonucleotide,  according  to  the  reverse 
transfection  procedure  in  the  manufacturer’s  instructions.  The  first  transfection  (day  0) 
into  HCT1 16  cells  was  followed  by  a  second  transfection  at  day  2.  Total  RNA  was 
isolated  and  various  assays  were  performed  at  day  4,  unless  otherwise  indicated.  During 
the  course  of  the  transfection  experiment,  HCT116  cells  were  kept  below  a  density  of 

n 

2X10  cells  per  10  cm  culture  dish. 

The  probe  sequences  for  Northerns  in  Fig  2B  are  5'-aaataagagcacccgcttc-3' 
(tRF-1001),  5'-tggcgcagcgagcagggttcgaa-3'  (mature  tRNA-Ser-TGA),  and  5'- 
gatcgggcgcgttcagggtggtat-3'  (5S  rRNA).  The  probe  for  snoU38b,  a  small  nucleolar 
RNA,  for  the  Northern  in  Fig  4A  is  5'-agaactggacaaagttttcatcac-3'. 

Measurement  of  cell  proliferation  and  DN  A  synthesis 

Cell  proliferation  was  measured  with  CellTiter  96  non-radioactive  cell  proliferation 
assay  kit  (MTT  assay;  Promega  Corp.).  DNA  synthesis  was  assayed  by 
bromodeoxyuridine  (BrdU)  ELISA  as  described  in  (Machida  et  al.  2006).  All 
experiments  were  performed  in  triplicates,  from  which  average  and  standard  deviation 
are  calculated  and  plotted. 

L  ucif erase  assays 


20 


Luciferase  assays  were  performed  as  described  in  (Kim  et  al.  2006;  Lee  and  Dutta  2007) 
and  Supplemental  Materials  and  Methods. 
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Table  Legends 


Table  1.  List  of  abundantly  cloned  tRFs 

17  tRFs  are  displayed  with  their  sequence,  length  (in  nucleotides),  cloning  frequency 
(number  of  sequence  reads  in  454  deep  sequencing  data),  and  relative  location  in  the 
corresponding  tRNA  (see  also  text  and  Fig  1  and  SI).  For  comparison,  we  included 
sequence  reads  of  two  miRNAs,  miR-18a  and  let-7a,  and  the  number  of  miRNAs  cloned 
more  than  or  less  than  200  times. 
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Figure  Legends 


Figure  1.  Classification  and  nomenclature  of  tRFs  revealed  by  analysis  of  highly 
abundant  nmsRNAs  (non-micro-small  RNAs). 

A.  Identification  and  classification  of  30  nmsRNA  sequences  that  were  cloned 
abundantly.  The  nmsRNAs  are  grouped  by  their  identity  in  the  pie  chart.  “snoRNA”, 
small  nucleolar  RNA;  “unmappable”:  nmsRNAs  unidentifiable  in  the  human  genome 
sequence  or  the  human  EST  database,  probably  due  to  as  yet  unreported  splicing  or 
sequence  editing  events.  tRNA-related  small  RNAs  (tRFs;  tRNA-derived  RNA 
Fragment)  are  further  classified  by  their  relative  location  (see  panel  B)  in  the  tRNA 
primary  transcript  (pre-tRNA). 

B.  Diagram  of  tRFs  aligned  at  a  tRNA  locus  (see  Fig  SI  for  actual  sequence  alignment). 
Wavy  line:  the  genomic  sequence  for  pre-tRNA  ending  with  an  oligo-dT  stretch,  the  3'- 
termination  signal  of  RNA  polymerase  III.  Grey  bar:  A  mature  tRNA  after  3 '-trimming 
of  the  pre-tRNA  and  addition  of  the  CCA  that  is  not  present  in  the  genomic  sequence. 
Black  bars:  Three  groups  of  tRFs  (tRF-5,  -3  and  -1  series)  are  aligned  with  the  pre- 
tRNA  or  mature  tRNA.  tRFs  are  numbered  with  the  first  digit  indicating  the  group  and 
the  other  digits  forming  a  serial  number  in  that  group. 
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Figure  2.  tRF-1001  is  highly  expressed  in  proliferating  cancer  cells. 

A.  Northern  hybridization  of  tRF-1001  (top  panel)  and  ethidium  bromide  staining  of  the 
gel  (bottom  panel).  Source  of  tissue  RNAs  is  described  in  (Kim  et  al.  2006).  tRF-1001 
(19-mer)  is  indicated  by  an  arrowhead  on  the  left.  The  positions  of  18-  and  24-mer 
synthetic  oligoribonucleotides  (arrowheads  on  the  right),  along  with  10-bp  DNA  ladder, 
are  shown  as  molecular  size  markers. 

B.  Measurement  of  tRF-1001  upon  serum  depletion  in  DU145  cell  line.  At  72  hrs  after 
serum  starvation  (see  Supplemental  Materials  and  Methods),  DU145  cells  were 
replenished  with  10%  FBS  medium  (“re-addition”).  Northern  hybridization  of  tRF-1001 
(bottom  panel),  its  mature  tRNA  (middle  panel),  and  5S  rRNA  as  a  loading  control  (top 
panel).  Different  panels  have  different  levels  of  exposure. 

Figure  3.  The  impaired  growth  of  HCT116  cells  upon  siRNA-mediated  knock¬ 
down  of  tRF-1001  is  rescued  by  Me-tRF-1001. 

A.  siRNA  duplex  was  designed  against  the  underlined  region  of  pre-tRNA 
encompassing  tRF-1001  and  the  mature  tRNA.  Arrowheads  indicate  expected  sites  in 
the  pre-tRNA  that  will  be  cleaved  by  the  siRNA  (Elbashir  et  al.  2001b). 

B.  qRT-PCR  measurement  of  tRF-1001  upon  siRNA  transfection.  Left  panel:  qRT-PCR 


27 


was  performed  with  gel-isolated  small  RNA  and  the  relative  expression  of  tRF-1001 


(normalized  by  tRF-5003)  is  indicated.  Right  panel:  total  RNA  was  used  for  qRT-PCR, 
with  snU6  for  normalization.  Each  value  is  an  average  of  triplicate  samples  and  the 
error  bars  indicate  standard  deviation.  The  normalized  level  of  tRF-1001  after  siGL2 

transfection  is  set  as  1. 

C.  MTT  and  BrdU  incorporation  assays  performed  after  transfection  of  siGL2  (plain 
bar),  si-tRFlOOl  (black  bar),  and  si-tRFlOOl  with  Me-GL2  (light  grey  bar)  or  Me-tRF- 
1001  (dark  grey  bar).  Left  panel:  A  scan  of  wells  after  MTT  assay  at  day  4.  Right  panel: 
BrdU  incorporation  normalized  to  MTT  value  at  days  3  and  4.  Average  and  standard 
deviation  from  triplicates  are  shown.  The  BrdU/MTT  for  cells  transfected  with  siGL2  is 
set  as  1 . 

D.  FACS  of  HCT116  cells  transfected  with  indicated  oligonucleotides  and  stained  with 
propidium  iodide  for  DNA  content  (4  days  after  first  transfection).  Cells  in  Gl,  S,  or  G2 
were  quantitated  by  ModFit  LT  for  Mac  (version  3.2.1).  Right  panels:  change  in  %  of 
cells  in  each  phase  of  the  cell-cycle  between  the  indicated  two  samples. 

Figure  4.  tRF-1001  is  localized  in  the  cytoplasm  but  is  not  incorporated  into 
miRNA/siRNA  pathway 
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A.  Total,  nuclear,  and  cytoplasmic  RNAs  were  subjected  to  Northern  hybridization  of 


tRF-1001,  snoU38  as  a  nuclear  marker,  and  5S  rRNA  as  a  loading  control. 

B.  Diagram  of  pRL-anti-tRFlOOl  which  contains  an  antisense  to  tRF-1001  downstream 
of  the  Renilla  luciferase  ORF  (open  reading  frame)  in  pRLCMV  vector  (Promega  Corp. 
and  (Kim  et  al.  2006)). 

C-E.  Luciferase  assays  were  performed  in  HCT116  cells.  Luciferase  values  were 
measured  at  24  hrs  post-transfection  of  the  plasmids.  Each  Renilla  luciferase  value  (Rr) 
was  normalized  to  that  of  firefly  (Pp)  from  the  co-transfected  pGL3-Control  vector 
(Promega  Corp.).  In  panel  C,  the  luciferase  plasmids  were  transfected  at  day  3  after  two 
rounds  of  siRNA  transfection  (at  day  0  and  2)  and  assays  were  performed  at  day  4  when 
knock-down  of  tRF-1001  was  confirmed  (see  Fig  3B).  In  panel  D,  serum  was  depleted 
at  the  time  of  transfection  so  that  assays  were  done  at  24  hrs  after  serum  depletion  when 
tRF-1001  was  significantly  reduced  (Fig  S5B).  In  panel  E,  transfection  of  indicated 
oligoribonucleotides  was  followed  by  transfection  of  the  luciferase  plasmids  at  24  hrs. 
GL2  value  in  each  pair  was  set  as  1 . 
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name 


sequence  (length) 


cloning  freq 
(seq  reads) 


tRNA  mapped 


tRF-1001 

GAAGCGGGTGCTCTTATTT 

19 

3995 

tRNA-Ser-TGA 

3'-out 

tRF-1002 

GCCGGGTACTTTCGTATTTT 

20 

1090 

tRNA-Asp-GTC 

3'-out 

tRF-1003 

GCTAAGGAAGTCCTGTGCTCAGTTT 

25 

378 

tRNA-Ser-GCT 

3'-out 

tRF-1004 

GTGTGTAGCTGCACTTTT 

18 

362 

tRNA-Asp-GTC 

3'-out 

tRF-3001 

ATCCCACCGCTGCCACCA 

18 

2852 

tRNA- Leu -TAG 

3' 

tRF-3002 

ACCCTGCTCGCTGCGCCA 

18 

672 

tRNA-Ser-TGA 

3' 

tRF-3003 

TCCCCGGCACCTCCACCA 

18 

657 

tRNA-Ala-AGC 

3' 

tRF-3004 

TCCCCGGCATCTCCACCA 

18 

775 

tRNA-Ala-TGC 

3' 

tRF-3005 

ATCCTGCCGACTACGCCA 

18 

1219 

tRNA-Ser-AGA 

3' 

tRF-3006 

TCGATTCCCGGCCCATGCACCA 

22 

358 

tRNA-Gly-GCC 

3' 

tRF-3007 

TCGATTCCCGGCCAACGCACCA 

22 

749 

tRNA-Gly-TCC 

3’ 

tRF-3008 

GTCCCACCAGAGTCGCCA 

18 

308 

tRNA-Arg-TCT 

3’ 

tRF-5001 

GACGAGGTGGCCGAGTGG 

18 

1797 

tRNA-Ser-GCT 

5’ 

tRF-5002 

TCCCTGGTGGTCTAGTGGTTA 

21 

906 

tRNA-Glu-CTC 

5’ 

tRF-5003 

GGTAGCGTGGCCGAGCGGTC 

20 

710 

tRNA- Leu- TAG 

5’ 

tRF-5004 

GCGTTGGTGGTATAGTGGT 

19 

1271 

tRNA-Gly-TCC 

5’ 

tRF-5005 

GCATGGGTGGTTCAGTGGTA 

20 

1221 

tRNA-Gly-GCC 

5’ 

microRNA 


miR-18a 


let- 7a 


635  miRNAs 


60  miRNAs 


seq  reads 


473 


17,856 
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Supplemental  Information  for  “A  novel  class  of  small  RNAs,  tRNA-derived  RNA 


fragments  (tRFs),  is  highly  expressed  in  proliferating  cells  and  required  for  cell 
proliferation” 

Supplemental  Materials  and  Methods 

Cell  allure 

Prostate  cancer  cell  lines  LNCaP  and  C4-2  (Wu  et  al.  1994)  were  maintained  in  RPMI- 
1640  media  supplemented  with  10%  FBS  (fetal  bovine  serum).  LNCaP,  an  androgen 
dependent  cell  line,  was  deprived  of  androgen  (“LNA-”  condition)  by  culturing  for  four 
days  in  a  phenol  red  free  RPMI-1640  medium  with  10%  charcoal/dextran  treated  FBS 
(Hyclone).  For  “LNDHT”  condition,  10  nM  4,5-a-dihydrotestosterone  (DHT;  from 
Sigma-Aldrich  Co.)  was  added  to  LNA-  condition. 

HCT116  colon  cancer  cell  line  and  DU  145  prostate  cancer  cell  line  were  maintained  in 
McCoy’s  5A  medium  (Iwakata  &  Grace  modification)  and  D-MEM  (Dulbecco’s 
Modification  of  Eagle’s  Medium;  4.5  g/L  glucose),  supplemented  with  10%  FBS 
respectively.  For  serum  starvation,  the  medium  was  replaced  with  D-MEM  (1.0  g/L 
glucose)  containing  0.1%  BSA  (bovine  serum  albumin;  fraction  V  from  Sigma-Aldrich 
Co.). 
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Staining  of  cells  urtfi  DAPI  endSP-DiOCisO) 


SP-DiOCig(3)  (a  sulfonated  derivative  of  DiO)  was  purchased  from  Molecular  Probes™. 
Cells  were  washed  with  PBS  (phosphate  buffered  saline)  twice  between  each  treatment. 
HCT116  cells  on  sterile  glass  coverslip  were  treated  with  SP-DiOCig(3)  for  5  min  at 
4°C  to  stain  cellular  membranes.  Cells  were  fixed  with  3.7%  (w/v)  HCHO  in  PBS  at 
37°C  for  10  min  and  were  permeabilized  with  chilled  acetone  at  -20°C  for  10  min. 

Nuclei  were  stained  with  DAPI  (H-1200)  (Vector  Laboratories,  Inc.)  for  1  min  before 
mounting.  Images  were  observed  under  Nikon  Microphot-SA  microscope,  captured  by 
Nikon  UFX-DX  camera,  and  processed  by  SPOT  (version  3.5.4  for  MacOS™)  software. 
Ludferaseassays 

An  artificial  target  site,  which  is  perfectly  complementary  to  tRF-1001,  was  cloned 
downstream  of  the  Renilla  lucif erase  ORF  (open  reading  frame)  in  pRL-CMV(MCS) 
(Kim  et  al.  2006),  to  construct  pRL-anti-tRFlOOl  (Fig  4B).  pRL-anti-tRFlOOl  or  pRL- 
CMV(MCS)  vector  control  was  co-transfected  with  pGL3-Control  vector  (Promega 
Corp.)  expressing  firefly  (Photinus  pyralis )  luciferase.  Plasmid  transfection  was 
performed  with  Lipofectamine  2000  reagent  (Invitrogen  Corp.)  according  to  the 
manufacturer’s  instructions.  Luciferase  assays  were  performed  at  24  hours  post¬ 
transfection  of  the  plasmids,  using  Dual-Luciferase  Reporter  Assay  System  (Promega 
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Corp.)  and  Pharmingen  Monolight™  3020  luminometer.  Each  value  of  Renilla 


luciferase  ( Rr )  was  normalized  to  the  firefly  luciferase  value  (Pp).  Each  value  is  an 
average  of  three  transfections  with  the  standard  deviation  indicated. 


Supplemental  Table  and  Figure  Legends 
Table  SI.  Number  of  tRNA-related  small  RNAs 

Among  non-micro- small  RNAs  (nmsRNAs)  that  were  cloned  more  than  5  times  and 
mappable  in  the  human  genome,  tRNA-related  sequences  (defined  as  sequences 
matching  anywhere  in  mature  tRNA  sequence  plus  the  flanking  25  nucleotides  at  both 
ends)  were  counted  and  tabulated.  tRNA-related  sequences  were  further  sorted  into  the 
three  tRFs  (tRNA-derived  RNA  Fragments;  defined  and  described  in  the  text  and  Fig  1 
and  SI)  or  the  rest  (“non-specific”).  The  numbers  are  plotted  in  pie  charts  in  Fig  S2. 

Table  S2.  tRFs  captured  in  the  454  high-throughput  sequencing 

All  tRFs  (cloned  >5  times)  are  tabulated  with  their  name,  sequence,  cloning  frequency 
(seq  reads),  and  corresponding  genomic  tRNA  loci  (chr  -  tRNA  #).  The  tRNA  number 
on  a  chromosome  is  based  on  the  tRNA  #  in  the  Genomic  tRNA  Database 
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(http://lowelab.ucsc.edu/GtRNAdb) 


Figure  SI.  Sequence  alignment  of  abundantly  cloned  tRFs  (tRNA-derived  RNA 
Eragment)  to  the  corresponding  tRNA 

Each  tRF  sequence  (bold  letters;  see  also  Table  1)  and  its  relative  position  (see  also  Fig 
1)  in  the  mature  tRNA  (shaded  capital  letters)  or  the  flanking  genomic  regions 
(lowercase  letters). 

Figure  S2.  Fraction  of  tRNA-related  sequences  in  nmsRNAs 

The  numbers  in  Table  SI  are  depicted  in  pie  charts.  See  Table  SI  legend  for  details. 

Figure  S3.  Length  distribution  of  tRF-1  series. 

For  622  tRNA  loci  in  the  tRNA  database,  3 '-end  of  pre-tRNAs  is  estimated  in  silico 
from  RNA  polymerase  III  termination  signal.  Pre-tRNAs  are  known  to  terminate  in 
oligo-U  stretches  (Hagenbuchle  et  al.  1979;  Koski  and  Clarkson  1982)  or  in  other  non- 
canonical  signals  such  as  UUCUU,  GUCUU  or  AUCUU  (Thomann  et  al.  1989).  3'-end 
of  pre-tRNA  is  defined  as  the  third  T  in  four  or  more  T  stretch  or  the  fourth  nucleotide 
in  TTCTT,  GTCTT  or  ATCTT.  Distribution  of  the  distances  between  3'-end  of  mature 
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tRNA  and  pre-tRNA  (=  estimated  lengths  of  tRF-1  series)  are  represented  in  a  pie  chart. 


The  portion  representing  91  tRNA  loci,  which  can  potentially  generate  tRF-1  series  of 
16-27  nucleotides  long  that  can  be  capture  by  our  cloning,  is  highlighted  by  grey  color. 

Figure  S4.  Detection  of  tRFs  by  Northern  hybridization  or  qRT-PCR. 

A.  Northern  hybridization  of  tRFs.  The  four  samples  are;  “LNCaP”  (cultured  in  a 
medium  with  untreated  serum),  “LNA-”  (LNCaP  cells  cultured  in  charcoal  stripped 
serum),  “LNDHT”  (LNA-  supplemented  with  a  synthetic  androgen  dihydrotestosterone 
(DHT)),  and  “C4-2”  (an  androgen  independent  cell  line  derived  from  LNCaP)  (see  also 
Supplemental  Materials  and  Methods).  Ethidium  bromide  (EtBr)  staining  of  total  RNA 
is  shown  for  equal  loading  (leftmost  panel).  18-  and  24-mer  oligoribonucleotides  (arrow 
heads)  and  10-bp  DNA  ladder  are  shown  as  molecular  size  markers. 

B.  tRFs  were  measured  by  qRT-PCR  as  described  in  Materials  and  Methods.  qRT-PCR 
value  of  each  tRF  was  normalized  by  that  of  miR-200b,  a  miRNA  similarly  expressed 
across  the  four  samples  (data  not  shown).  Each  value  is  an  average  of  triplicate  samples, 
with  the  standard  deviation  indicated. 

Figure  S5.  Decrease  of  tRF-1001  upon  serum  depletion. 
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A.  Northern  hybridization  of  tRF-1001  upon  serum  starvation  in  LNCaP  cells.  EtBr 
staining  of  total  RNA  is  shown  for  equal  loading. 

B  and  C.  tRF-1001  measurement  in  HCT116  cell  line  upon  serum  starvation  (panel  B) 
and  different  cell  density  (panel  C).  In  panel  B,  the  number  of  cells  at  the  time  of  serum 
starvation  was  -1X107  (in  10-cm  culture  dish).  In  panel  C,  number  of  cells  in  10-cm 
culture  dish  was  counted  at  the  time  of  RNA  isolation.  All  other  descriptions  are  same 
as  Fig  2B. 

Figure  S6.  Comparison  of  qRT-PCR  and  Northern  hybridization  in  the 
measurement  of  tRF-1001. 

A.  Northern  hybridization  of  miR-21  and  tRF-1001,  as  well  as  EtBr  staining  of  total 
RNA  as  a  loading  control.  tRF-1001  result  was  reclaimed  from  Fig  S5B  (rightmost  lane) 
and  S5C  (left  lane)  for  side-by-side  comparison. 

B.  qRT-PCR  measurement  of  tRF-1001.  The  data  was  normalized  by  miR-21  (left 
panel)  or  tRF-5003  (right  panel).  Plain  bar  is  set  as  1.  All  other  descriptions  are  same  as 
Fig  S4B. 


Figure  S7.  Level  of  the  mature  tRNA-Ser-TGA  was  unaffected  after  transfection  of 


si-tRFlOOl. 


Northern  hybridization  of  tRNA-Ser-TGA  (top  panel)  and  EtBr  staining  of  total  RNA 
(bottom  panel). 

Figure  S8.  Cell  morphology  after  tRF-1001  knock-down 

Cellular  contour  (cellular  membranes  stained  green  by  SP-DiOCig(3))  and  nuclei 
(stained  blue  by  DAPI)  are  visualized  as  described  in  Supplemental  Materials  and 
Methods. 

Figure  S9.  Reduced  growth  upon  si-tRFlOOl  treatment  was  not  rescued  by  another 
control  2'-0-methyl  oligonucleotide  whose  sequence  overlaps  si-tRFlOOl  by  12 
nucleotides 

A.  Sequence  of  Me-tRF-1001  and  Me-control.  The  portions  overlapping  si-tRFlOOl 
(shown  in  the  top)  are  underlined. 

B.  MTT  assay  is  performed  as  described  in  Materials  and  Methods.  The  value  of  siGL2 
is  set  as  1 . 
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Supplemental  Information  for  “A  novel  class  of  small  RNAs,  tRNA-derived  RNA 


fragments  (tRFs),  is  highly  expressed  in  proliferating  cells  and  required  for  cell 
proliferation” 

Supplemental  Materials  and  Methods 

Cell  allure 

Prostate  cancer  cell  lines  LNCaP  and  C4-2  (Wu  et  al.  1994)  were  maintained  in  RPMI- 
1640  media  supplemented  with  10%  FBS  (fetal  bovine  serum).  LNCaP,  an  androgen 
dependent  cell  line,  was  deprived  of  androgen  (“LNA-”  condition)  by  culturing  for  four 
days  in  a  phenol  red  free  RPMI-1640  medium  with  10%  charcoal/dextran  treated  FBS 
(Hyclone).  For  “LNDHT”  condition,  10  nM  4,5-a-dihydrotestosterone  (DHT;  from 
Sigma-Aldrich  Co.)  was  added  to  LNA-  condition. 

HCT116  colon  cancer  cell  line  and  DU  145  prostate  cancer  cell  line  were  maintained  in 
McCoy’s  5A  medium  (Iwakata  &  Grace  modification)  and  D-MEM  (Dulbecco’s 
Modification  of  Eagle’s  Medium;  4.5  g/L  glucose),  supplemented  with  10%  FBS 
respectively.  For  serum  starvation,  the  medium  was  replaced  with  D-MEM  (1.0  g/L 
glucose)  containing  0.1%  BSA  (bovine  serum  albumin;  fraction  V  from  Sigma-Aldrich 
Co.). 
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Staining  of  cells  urtfi  DAPI  endSP-DiOCisO) 


SP-DiOCig(3)  (a  sulfonated  derivative  of  DiO)  was  purchased  from  Molecular  Probes™. 
Cells  were  washed  with  PBS  (phosphate  buffered  saline)  twice  between  each  treatment. 
HCT116  cells  on  sterile  glass  coverslip  were  treated  with  SP-DiOCig(3)  for  5  min  at 
4°C  to  stain  cellular  membranes.  Cells  were  fixed  with  3.7%  (w/v)  HCHO  in  PBS  at 
37°C  for  10  min  and  were  permeabilized  with  chilled  acetone  at  -20°C  for  10  min. 

Nuclei  were  stained  with  DAPI  (H-1200)  (Vector  Laboratories,  Inc.)  for  1  min  before 
mounting.  Images  were  observed  under  Nikon  Microphot-SA  microscope,  captured  by 
Nikon  UFX-DX  camera,  and  processed  by  SPOT  (version  3.5.4  for  MacOS™)  software. 
Ludferaseassays 

An  artificial  target  site,  which  is  perfectly  complementary  to  tRF-1001,  was  cloned 
downstream  of  the  Renilla  lucif erase  ORF  (open  reading  frame)  in  pRL-CMV(MCS) 
(Kim  et  al.  2006),  to  construct  pRL-anti-tRFlOOl  (Fig  4B).  pRL-anti-tRFlOOl  or  pRL- 
CMV(MCS)  vector  control  was  co-transfected  with  pGL3-Control  vector  (Promega 
Corp.)  expressing  firefly  (Photinus  pyralis )  luciferase.  Plasmid  transfection  was 
performed  with  Lipofectamine  2000  reagent  (Invitrogen  Corp.)  according  to  the 
manufacturer’s  instructions.  Luciferase  assays  were  performed  at  24  hours  post¬ 
transfection  of  the  plasmids,  using  Dual-Luciferase  Reporter  Assay  System  (Promega 
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Corp.)  and  Pharmingen  Monolight™  3020  luminometer.  Each  value  of  Renilla 


luciferase  ( Rr )  was  normalized  to  the  firefly  luciferase  value  (Pp).  Each  value  is  an 
average  of  three  transfections  with  the  standard  deviation  indicated. 


Supplemental  Table  and  Figure  Legends 
Table  SI.  Number  of  tRNA-related  small  RNAs 

Among  non-micro- small  RNAs  (nmsRNAs)  that  were  cloned  more  than  5  times  and 
mappable  in  the  human  genome,  tRNA-related  sequences  (defined  as  sequences 
matching  anywhere  in  mature  tRNA  sequence  plus  the  flanking  25  nucleotides  at  both 
ends)  were  counted  and  tabulated.  tRNA-related  sequences  were  further  sorted  into  the 
three  tRFs  (tRNA-derived  RNA  Fragments;  defined  and  described  in  the  text  and  Fig  1 
and  SI)  or  the  rest  (“non-specific”).  The  numbers  are  plotted  in  pie  charts  in  Fig  S2. 

Table  S2.  tRFs  captured  in  the  454  high-throughput  sequencing 

All  tRFs  (cloned  >5  times)  are  tabulated  with  their  name,  sequence,  cloning  frequency 
(seq  reads),  and  corresponding  genomic  tRNA  loci  (chr  -  tRNA  #).  The  tRNA  number 
on  a  chromosome  is  based  on  the  tRNA  #  in  the  Genomic  tRNA  Database 
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(http://lowelab.ucsc.edu/GtRNAdb) 


Figure  SI.  Sequence  alignment  of  abundantly  cloned  tRFs  (tRNA-derived  RNA 
Eragment)  to  the  corresponding  tRNA 

Each  tRF  sequence  (bold  letters;  see  also  Table  1)  and  its  relative  position  (see  also  Fig 
1)  in  the  mature  tRNA  (shaded  capital  letters)  or  the  flanking  genomic  regions 
(lowercase  letters). 

Figure  S2.  Fraction  of  tRNA-related  sequences  in  nmsRNAs 

The  numbers  in  Table  SI  are  depicted  in  pie  charts.  See  Table  SI  legend  for  details. 

Figure  S3.  Length  distribution  of  tRF-1  series. 

For  622  tRNA  loci  in  the  tRNA  database,  3 '-end  of  pre-tRNAs  is  estimated  in  silico 
from  RNA  polymerase  III  termination  signal.  Pre-tRNAs  are  known  to  terminate  in 
oligo-U  stretches  (Hagenbuchle  et  al.  1979;  Koski  and  Clarkson  1982)  or  in  other  non- 
canonical  signals  such  as  UUCUU,  GUCUU  or  AUCUU  (Thomann  et  al.  1989).  3'-end 
of  pre-tRNA  is  defined  as  the  third  T  in  four  or  more  T  stretch  or  the  fourth  nucleotide 
in  TTCTT,  GTCTT  or  ATCTT.  Distribution  of  the  distances  between  3'-end  of  mature 
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tRNA  and  pre-tRNA  (=  estimated  lengths  of  tRF-1  series)  are  represented  in  a  pie  chart. 


The  portion  representing  91  tRNA  loci,  which  can  potentially  generate  tRF-1  series  of 
16-27  nucleotides  long  that  can  be  capture  by  our  cloning,  is  highlighted  by  grey  color. 

Figure  S4.  Detection  of  tRFs  by  Northern  hybridization  or  qRT-PCR. 

A.  Northern  hybridization  of  tRFs.  The  four  samples  are;  “LNCaP”  (cultured  in  a 
medium  with  untreated  serum),  “LNA-”  (LNCaP  cells  cultured  in  charcoal  stripped 
serum),  “LNDHT”  (LNA-  supplemented  with  a  synthetic  androgen  dihydrotestosterone 
(DHT)),  and  “C4-2”  (an  androgen  independent  cell  line  derived  from  LNCaP)  (see  also 
Supplemental  Materials  and  Methods).  Ethidium  bromide  (EtBr)  staining  of  total  RNA 
is  shown  for  equal  loading  (leftmost  panel).  18-  and  24-mer  oligoribonucleotides  (arrow 
heads)  and  10-bp  DNA  ladder  are  shown  as  molecular  size  markers. 

B.  tRFs  were  measured  by  qRT-PCR  as  described  in  Materials  and  Methods.  qRT-PCR 
value  of  each  tRF  was  normalized  by  that  of  miR-200b,  a  miRNA  similarly  expressed 
across  the  four  samples  (data  not  shown).  Each  value  is  an  average  of  triplicate  samples, 
with  the  standard  deviation  indicated. 

Figure  S5.  Decrease  of  tRF-1001  upon  serum  depletion. 
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A.  Northern  hybridization  of  tRF-1001  upon  serum  starvation  in  LNCaP  cells.  EtBr 
staining  of  total  RNA  is  shown  for  equal  loading. 

B  and  C.  tRF-1001  measurement  in  HCT116  cell  line  upon  serum  starvation  (panel  B) 
and  different  cell  density  (panel  C).  In  panel  B,  the  number  of  cells  at  the  time  of  serum 
starvation  was  -1X107  (in  10-cm  culture  dish).  In  panel  C,  number  of  cells  in  10-cm 
culture  dish  was  counted  at  the  time  of  RNA  isolation.  All  other  descriptions  are  same 
as  Fig  2B. 

Figure  S6.  Comparison  of  qRT-PCR  and  Northern  hybridization  in  the 
measurement  of  tRF-1001. 

A.  Northern  hybridization  of  miR-21  and  tRF-1001,  as  well  as  EtBr  staining  of  total 
RNA  as  a  loading  control.  tRF-1001  result  was  reclaimed  from  Fig  S5B  (rightmost  lane) 
and  S5C  (left  lane)  for  side-by-side  comparison. 

B.  qRT-PCR  measurement  of  tRF-1001.  The  data  was  normalized  by  miR-21  (left 
panel)  or  tRF-5003  (right  panel).  Plain  bar  is  set  as  1.  All  other  descriptions  are  same  as 
Fig  S4B. 


Figure  S7.  Level  of  the  mature  tRNA-Ser-TGA  was  unaffected  after  transfection  of 


si-tRFlOOl. 


Northern  hybridization  of  tRNA-Ser-TGA  (top  panel)  and  EtBr  staining  of  total  RNA 
(bottom  panel). 

Figure  S8.  Cell  morphology  after  tRF-1001  knock-down 

Cellular  contour  (cellular  membranes  stained  green  by  SP-DiOCig(3))  and  nuclei 
(stained  blue  by  DAPI)  are  visualized  as  described  in  Supplemental  Materials  and 
Methods. 

Figure  S9.  Reduced  growth  upon  si-tRFlOOl  treatment  was  not  rescued  by  another 
control  2'-0-methyl  oligonucleotide  whose  sequence  overlaps  si-tRFlOOl  by  12 
nucleotides 

A.  Sequence  of  Me-tRF-1001  and  Me-control.  The  portions  overlapping  si-tRFlOOl 
(shown  in  the  top)  are  underlined. 

B.  MTT  assay  is  performed  as  described  in  Materials  and  Methods.  The  value  of  siGL2 
is  set  as  1 . 


7 


Supplemental  References 

Hagenbuchle,  O.,  Larson,  D.,  Hall,  G.I.,  and  Sprague,  K.U.  1979.  The  primary 

transcription  product  of  a  silkworm  alanine  tRNA  gene:  identification  of  in  vitro 
sites  of  initiation,  termination  and  processing.  Cell  18(4):  1217-1229. 

Kim,  H.K.,  Lee,  Y.S.,  Sivaprasad,  U.,  Malhotra,  A.,  and  Dutta,  A.  2006.  Muscle-specific 
microRNA  miR-206  promotes  muscle  differentiation.  J  Cell  Biol  174(5):  677- 
687. 

Koski,  R.A.  and  Clarkson,  S.G.  1982.  Synthesis  and  maturation  of  Xenopus  laevis 

methionine  tRNA  gene  transcripts  in  homologous  cell-free  extracts.  J  Biol  Chem 
257(8):  4514-4521. 

Thomann,  H.U.,  Schmutzler,  C.,  Hudepohl,  U.,  Blow,  M.,  and  Gross,  H.J.  1989.  Genes, 
variant  genes  and  pseudogenes  of  the  human  tRNA(Val)  gene  family.  Expression 
and  pre-tRNA  maturation  in  vitro.  J  Mol  Biol  209(4):  505-523. 

Wu,  H.C.,  Hsieh,  J.T.,  Gleave,  M.E.,  Brown,  N.M.,  Pathak,  S.,  and  Chung,  L.W.  1994. 
Derivation  of  androgen-independent  human  LNCaP  prostatic  cancer  cell 
sublines:  role  of  bone  stromal  cells.  Int  J  Cancer  57(3):  406-412. 


